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TITLE: PROJECT MANAGEMENT 

Principal Investigator: Dr. Tim Anderson, Professor and Associate Dean 
Program Manager: Canan Balaban, Chemical Engineer, M.Sc. 

GOALS 

Project management activities for this grant includes internal oversight, coordination and 
facilitation among research groups (internal to UF and external); system and/or trade studies to 
assist in evaluating and prioritizing new hydrogen production processes that may become 
available for potential research investigations; a technical and programmatic outreach activity; 
and formal reporting, including periodic technical reviews with the NASA team. 

ACCOMPLISHMENTS 

Proposal 

The original proposal including the budget was revised based on the results of the NASA 
rankings after the NASA August 2003 review meeting. The revised proposal was reviewed with 
the NASA technical review team at GRC to provide in-depth briefings on the content of the 
various tasks and to assure that the programmatic content of this grant is fully understood and 
supported by the NASA review team. 

Review Meetings 

1- September 26, 2003 Review Meeting: “Laser Instrumentation for Hydrogen Leak Detection, 
Rayleigh/Mie Scattering” was reviewed at Glenn Research Center with Dr. Richard Seasholtz, 
Amy Mielke, and James Burkhart. 

2- November 14, 2003 Review Meeting: A NASA review meeting was held in Gainesville 
Florida, for the following projects: 

• Instrumentation for Hydrogen Leak Detection, Raman and Rayleigh Scattering 

• Sensors for Hydrogen Leak Detection 

• The Ammonia- Water Combined Power and Refrigeration Cycle 

• Advanced Protonic Conductors for Hydrogen Production 

3- March 4-5, 2004 Review Meeting: Held at the Florida State University, Tallahassee Florida 

4- March 30-31, 2004 Review Meeting: Held at the University of Central Florida, Orlando 
Florida. 

5- November 3-4, 2004 Review Meeting was held at the Florida Solar Energy Center, Cocoa 
Beach, Florida. 

All the review presentations are posted at http://www.mae.ufl.edu/NASAHvdrogenResearch 
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Reporting 

Five progress reports were submitted to NASA Glenn Research Center. The reporting dates are 
given below: 

1 - December 3 1 , 2003 

2- April 30, 2004 

3- September 30, 2004 

4- January 3 1 , 2005 

5- June 30, 2005 

In addition to the above progress reports a report on “NASA Grants Impact” was submitted in 
May, 2005. All the reports are posted at http://www.mae.ufl.edu/NASAHvdrogenResearch 

Fuel Cell Workshop 

A UF interdisciplinary fuel cell research planning activity was initiated in February 2004. The 
focus of the planning activity was hydrogen-based fuel cells for space, aviation, and ground 
based applications. This activity was supplemented by a UF Advanced Fuel Cell Workshop. 

The UF Advanced Fuel Cell Workshop was held on May 27-28 in Gainesville, FL. A total of 92 
people participated. 

Allan Bedwell, Deputy Secretary for Regulatory Programs for the Florida Department of 
Environmental Protection, was the invited opening speaker. The other speakers were: 

- Anita Liang, Chief, Aero-propulsion Projects Office, NASA Glenn Research Center 

- Dr. Ajay K. Misra, Chief, Ceramics Branch, NASA Glenn Research Center 

- Katrina Fritz, Interim Managing Director, Power Partnership for Ohio 

- Mr. Wayne Surdoval, SECA Coordinator, National Energy Technology Laboratory, U. S. 
Department of Energy 

- Dr. David Schiraldi, Macromolecular Science & Engineering, Case Western Reserve 
University 

- Dr. Hansung Kim, University of South Carolina 

- Dr. Eric Wachsman, Department of Materials Science & Engineering, University of Florida 

- Dr. Lisa McElwee- White, Department of Chemistry, University of Florida 

Invention Disclosures / Patent Applications 

The project management functions facilitated the applications for the invention disclosures and 
NASA new technology reports on: 

• Mass Gauging of Cryogenic Fluids in Tanks Experiencing Arbitrary Gravity or 
Acceleration 

• Fluid Distribution for In-Space Cryogenic Propulsion 

Dr. Hugh Fan filed a patent application on “Hydrogen Sensor Using Enzyme-catalyzed 
Reaction” 
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Other Funding Opportunities 

UF continues to identify and compete for other government and private funding opportunities in 
the hydrogen and fuel cell area. The project manager forwards to the faculty, relevant Requests 
for Proposals (RFP’s) from federal, state, and private funding sources. She also identifies and 
facilitates collaborative research possibilities among various faculty by meeting with faculty that 
may eventually submit competitive proposals. Several DOE funding opportunities were shared 
with the interested faculty members. 

Dr. Eric Wachsman received two DOE awards in the areas of “Mixed Ionic-electronic Conductor 
Electrodes” and “Fuel Cell Cathodes”. The total award is $1.25M. 

Outreach 

The Education/outreach effort consists of a website. The website documents and describes 
hydrogen research activities carried out at the University of Florida and other universities that 
may be collaborating with UF. All researchers involved, their research activities, and 
publications reside in the website. The hydrogen technology web course is also posted at the web 
site. The website connects all the researchers and improves communications within the program 
and with external hydrogen and energy research communities. The target audience for this web 
site is research professionals, students, the public, and government agencies. 

The web site is functional and in use. In addition to the NASA Supported Hydrogen Research, a 
web site on “UF Faculty Expertise in Energy” was developed. The links for both sites are given 
below: 

http://www.mae.ufl.edu/NASAHvdrogenResearch 

http://www.mae.ufl.edu/Energy 
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PROJECT TITLE: SYSTEM AND TRADE STUDIES 

Consultant: Dr. Dale W. Kirmse 

BACKGROUND 

The system and trade studies were performed for the Ammonia-Water Combined Power and 
Refrigeration Cycle. 

Dr. Dale Kirmse, recently retired Professor from UF’s Chemical Engineering Department, 
worked as a consultant to perform simulation, system/technology trades studies. ASPEN 
chemical process modeling software was used to perform the studies. 

ACCOMPLISHMENTS 

Since FSEC has a number of hydrogen production processes under this grant, the system and 
trade study efforts were coordinated with them. UF team visited FSEC group on Dec 08, 04 to 
establish comparison criteria for the hydrogen production processes under this grant. 

Dr. Kirmse submitted 4 reports. A summary of each report is given below. 

Report 1: Ammonia- Water Combined Cycle - Model 1 

Summary: Initial steady state ASPEN models of the Ammonia- Water Combined Power and 
Refrigeration Cycle process have been developed. Simple process modules and simplifying 
operation assumptions were used in this initial model. Heater modules were used to simulate heat 
exchangers and a stream splitter block was used to simulate the rectifier. Simplifying 
assumptions included no pressure losses through operating equipment. The recommended 
computational sequence is to break the cycle between the boiler and separation operations to 
form an equivalent open cycle, steady state model. The process flow sheet is given below. 



Figure 1: Ammonia-Water Combined Cycle Model 1.0 Flow sheet 
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Report 2: Ammonia-Water Combined Cycle - Model 2 

Summary: An open cycle, steady state ASPEN model of the Ammonia- Water Combined Power 
and Refrigeration Cycle has been developed. The simplistic modules used in the initial model 
have been replaced with rigorous modules. This model represents the demonstration facility 



Figure 2: Ammonia-Water Combined Cycle Model 2.0 Flow sheet 


Report 3: Ammonia- Water Combined Cycle - Response Exploration 

Summary: The responses to changes in parameter values in the open cycle, steady state ASPEN 
Model 2 of the Ammonia- Water Combined Power and Refrigeration Cycle have been studied. 
The sensitivity study found that four of the independent variables (the ammonia fraction of the 
strong aqueous solution, absorber pressure, fractionator pressure, and reboiler temperature) 
evoked strong responses to cycle performance. The process flow sheet and turbine exhaust 
temperature sensitivity to pressure are given below. 
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Figure 3: Ammonia-Water Combined Cycle Model 2.1 Flow sheet 



Figure 4: Turbine exhaust temperature sensitivity to high pressure section pressure 

Report 4: Ammonia- Water Combined Cycle - Module/Process Analysis 

Summary: Steady state ASPEN models of individual unit operations and processes of the 
Ammonia- Water Combined Power and Refrigeration Cycle were developed and analyzed. Each 
of the modules was inserted on a flow sheet and tested to make sure that the individual blocks 
and streams behaved as expected. 

The reports submitted by Dr. Kirmse are very detailed. Please contact Canan Balaban 
('cbalaban@ufl.edul to request the detailed reports. 
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HYDROGEN PRODUCTION PROCESSES 

Faculty Lead: Dr. Skip Ingley, Mechanical &Aerospace Engineering 

PROJECT TITLE: THE AMMONIA -WATER COMBINED POWER & 
REFRIGERATION CYCLE 

Task PI: Dr. H.A. (Skip) Ingley, P.E., Mechanical & Aerospace Engineering 
Graduate Students: Robert Reed, Gerard O’Sullivan, Mark Taurence, Numpol 
Sathonpattanakij, Allen Mathis, Sirisha Govindaraju 

PROJECT GOALS 

The project goal is to design and construct a combined cycle that creates 5kW of electricity, a 
refrigeration effect and 1 8 L/min hydrogen gas from a low temperature heat source and to house 
the system within a mobile test facility. The mobile facility will serve as an educational tool to 
demonstrate various uses for hydrogen and methods for its production around the state. 

ACCOMPLISHMENTS 

The project includes several sub-tasks: a 1st law analysis of the system, the design and 
construction of the combined cycle shown in Figure 1, including the heating and cooling water 
systems, and the design of the mobile test facility with a viewing area and a demonstration of 
how the device can produce hydrogen and cooling, with subsequent conversion to electricity 
with a fuel cell. 

Constraints have been placed on pressures and temperatures of the system to both simulate a 
low-temperature heat source such as that available from flat-plate solar collectors or waste heat 
and produce 5kW of electricity in the generator. The high and low cycle pressures are defined as 
1 lOpsig and 40 psig. The temperatures exiting the boiler/generator and absorber are fixed at 
170°F and 100°F, respectively. Results from the 1st law analysis for the thermodynamic analysis 
show a required heat input 387.6 MBTUH and heat rejection of 372.9 MBTUH to produce 4600 
BTU/h of cooling and 5.6 kW of shaft power at the turbine (0.6 kW of which is required to 
overcome the inefficiency of the generator). Assumptions for the 1st law analysis are zero 
pressure drops through the pipes and components, negligible pump work, 80% isentropic turbine 
efficiency, 90% electric generator efficiency, and pure ammonia vapor exiting the rectifier. 

The 1st law analysis was used to design each component of the combined cycle. Design is given 
below for each of the following: absorber, vapor generator/separator, rectifier, super-heater, 
cooler, turbine/generator set, electrolyzer, power conditioner, fuel cell, and process water 
systems (heating hot water and cooling water). 

• Absorber: A falling film absorber is being used due to the large liquid to vapor interface and 
large cooling area. It is constructed from flat plate heat exchanger plates, such as those 
commonly used in the dairy industry and thus readily available. The pressure vessel is 
constructed from mild steel. An alternative absorber was also designed using a packed 
column/heat exchanger arrangement. The basic cycle can accommodate either of these absorber 


9 


UNIVERSITY OF 

^FLORIDA 


packages. The study and design of the falling film absorber was the subject of a Master’s thesis 
for Sirisha Govindaraju who graduated in 2004. 

• Vapor Generator/Separator: The vapor generator is a shell and tube type heat exchanger. A 
commercially available heat exchanger has been purchased and modified to serve as the vapor 
generator. 

• Rectifier: The rectifier consists of a vertical length of pipe packed with a material such as burl 
saddles, followed by a cooling coil. The purity of the ammonia vapor is critical to the cooling 
capacity, so the rectifier stack has a modular construction for simple modification. It is mounted 
directly above the vapor generator to allow the reflux to be re-boiled. 

• Super-heater: As superheating is detrimental to cooling capacity, the superheater is used to 
reheat the rectified vapor up to the original boiler temperature. It is part of the rectifier stack, and 
identical in construction to the cooling coil. 

• Cooler: The cooler is used to cool air using the cold vapor at the turbine exit. An intercooler 
from a small diesel truck is used as an air-to-air heat exchanger. 

• Turbine: The turbine is of a scroll type because of its compact size and steady power output. A 
scroll compressor was purchased and modified to operate as an expander. Testing of this device 
has been conducted to determine its efficiency and operating characteristics. Any efficiency 
below 80% will cause the complete loss of cooling capability for the specified operating 
conditions. The scroll compressor has proven to be inadequate for the efficiencies required. A 
new scroll expander was designed to provide the efficiencies required. The subject of this 
investigation led to a Master’s thesis by Robert Reed (2005). This work was recently presented 
to the 2005 World Solar Congress held in Orlando Florida and is published in the Congress 
proceedings. A copy of this paper is attached to this report. 

• Electric Generator: The generator provides 400 Hz power, requiring the researchers to provide 
power conditioning. 

• Electrolyzer: The electrolyzer is required to produce hydrogen with a standard purity of 99.9% 
at 40 psig operating pressure. The IMET Electrolysis from Hydrogen System Inc. has been used 
as a model for calculations. However its cost was prohibitive. An alternate electrolyzer has been 
purchased and tested. The operating pressure and temperature ranges are 50-150 psig and 36- 
104°F, respectively. The unit is capable of delivering hydrogen at 99.9% purity as predicted from 
calculations. A functionl test has been successfully completed by one our graduate students, 

Mark Taurence. Mark has connected the electrolyzer to the Ballard fuel cell demonstration and 
demonstrated its ability to operate the various components of the demonstration unit. 

• Fuel Cell: The fuel cell selected to utilize the gaseous hydrogen produced by the electrolyzer is 
the Ballard Nexa. This PEM fuel cell provides 1.2 kW of DC power with a voltage range of 22- 
50 V. The fuel cell is on site and has already been used in a study on the effects of air 
contaminants on PEM fuel cell performance. The subject of this study provided a Master’s 
thesis for Gerard O’Sullivan (2005). A paper on this subject was presented and published in the 
proceedings of the 2005 World Solar Congress. A copy of the paper is attached. 

• Process water systems: The boiler for the heating hot water system was selected to satisfy the 
calculated heat input requirement, leaving extra capacity to allow for reheating of the cooling 
water by a maximum of 5°F and for an estimated 10% of losses through the piping. The cooling 
tower was also selected for the cooling water system to satisfy and exceed the heat rejection 
requirement, as well as the load and size constraints. A remote sump option for the cooling tower 
was selected to reduce the operating weight by locating the majority of the water outside of the 
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facility. The assumption used to calculate the required flow rates of the process water systems is 
a 10° AT for the water through the absorber, rectifier, vapor generator, and super-heater. 

Both items have been installed in the mobile test facility. The piping for the cooling water 
system is schedule-40 PVC for the cooling water and schedule-40 steel with fiberglass pipe 
insulation and an ASJ jacket for the heating hot water system. The layout, sizing and installation 
of the pipe and pumps is completed. The system is currently being tested. The process water 
flow diagram is given in Figure 2. 

• Mobile Test Facility: 

Overview: The mobile test facility is a modified mobile building unit to house the 5 kW 
combined cycle with the capability of being moved/towed from place to place. The trailer is able 
to support the total load of all the components of the combined cycle. 

Design and Layout: A 1-hour fire-rated mechanical room has been built at the end of the trailer 
to house the boiler as required. 

Fire Protection: The University of Florida building codes do not require fire sprinklers. Carbon 
monoxide, ammonia, and smoke detectors as well as ammonia respirators and fire extinguishers 
are being provided at different locations. At least two means of egress have been provided 
Selection: The trailer size is 14 ft x 60 ft and 8 ft tall. All interior wall sections with view panels 
will be provided by the researchers. A schematic of the trailer is given in Figure 3. The trailer 
includes a 100 panel and a through-the-wall air conditioning unit. 

• Hydrogen/Fuel cell demonstration: The portion of the system that illustrates the application of 
hydrogen technology is the fuel cell demonstration. A Ballard Nexa Fuel Cell uses the hydrogen 
provided by the electrolyzer. The fuel cell provides a net output of 1 .2 kW with a DC voltage 
range from 22-50 Volts. This output will power several components on a graphical 
representation of a vehicle. 

A large decal of an outline of a car has been prepared. This is visible from the viewing area of 
the trailer. Within this decal the fuel cell and all necessary connections are visible. The fuel cell 
is connected to three motors and two lamps. The two motors rotate the wheels of the car, a fan 
motor blows cool air into the viewing area, and two halogen lamps represent the rear and head 
lights. The overall purpose of this display is to show the usefulness of the fuel cells’ power. 

The display components have been purchased and assembled. 

Two thesis topics have resulted from this work leading to the graduation of two Masters level 
students in Mechanical and Aerospace Engineering. Another graduate student working on this 
project acquired a non-thesis Masters degree. The two thesis topics included research of the 
effects of air contaminants on the performance of a PEM fuel cell (“The Effect of Ambient Air 
Contaminants on the Performance of a Proton Exchange Membrane Fuel Cell”, Gerard Michael 
O’Sullivan, 2005) and the study of a scroll expander as applied to the combined ammonia/water 
combined cycle (“Study of the Feasibility and Energy Savings of Producing and Pre-Cooling 
Hydrogen with a 5kW Ammonia Based Combined Power/Cooling Cycle”, Robert Reed, 2005). 
Currently three other graduate students are preparing to continue these studies with focus on the 
control and optimization of the 5 kW ammonia/water combined cycle and the destructive effects 
of urban air contaminants on PEM fuel cells. 
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Figure 1. Schematic of the combined cycle showing the state point numbers, system components, 
and flow directions 
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Figure 3. Mobile test facility floor plan, 


Color Code: 

Blue - Cooling water system components 
Red — Heating hot water system components 
Green - Combined cycle components 
Violet — Fuel cell demonstration components 


Figure 4. Mobile Test Facility 











Figure 5. Scroll Expander 


Figure 5a. Scroll Expander 


Figure 5b. Scroll Expander Test Set-up 


Figure 6. Cooling Tower 
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PROJECT TITLE: ADVANCED PROTONIC CONDUCTORS 

Task PI: Dr. Eric. D. Wachsman, Department of Materials Science and Engineering 


PROJECT GOALS 

Membrane reactors can be used to produce hydrogen from hydrocarbon feedstocks, from natural 
gas to coal to a variety of hydrocarbon containing waste streams such as landfill gas and swamp 
gas. The overall objective of this project is to demonstrate the feasibility of producing hydrogen 
from hydrocarbon based fuels using advanced proton conducting membrane. The objectives of 
this phase are to develop thin film proton conducting membranes on porous supports and to 
advance our fundamental understanding of these materials. The anticipated result is 
demonstration of high hydrogen fluxes through these thin supported membranes. 


BACKGROUND 

Hydrogen is the most abundant chemical-energy resource in the world. Historically, hydrogen 
has been produced from fossil fuels and for the foreseeable future will continue to be so. The 
major source of H2 is steam reformation of natural gas. Therefore, improvements in the 
efficiency and cost of H2 production from natural gas are necessary in the near term. Gas 
separation membranes and membrane reactors based on ion conducting ceramics may provide 
the technological advance necessary to increase the efficiency and reduce the cost of H2 
production from natural gas. However, other sources of H2 must be developed for the envisioned 
hydrogen economy, and coal provides the greatest U.S. domestic resource-based option. The 
U.S. DOE is developing a FutureGen plant based on coal gasification, solid oxide fuel cells 
(SOFCs), and ion conducting membranes that will produce H2 and electricity with zero 
emissions and carbon sequestration; thereby, not contributing to global warming. 

Membrane reactor technology holds the promise to circumvent thermodynamic equilibrium 
limitations by in situ removal of product species, resulting in improved chemical yields. Mixed- 
conducting oxide-membrane technology presents the possibility for a dramatic reduction in the 
cost of converting petroleum and coal derived feed stocks to hydrogen and other value-added 
hydrocarbons. The membranes are based on metal oxides (e.g., SrCeiJvLCb, where M is a metal 
dopant, e.g. Eu) that exhibit both ionic and electronic (mixed) conductivity. Because of their 
significant electronic conductivity, these mixed ionic-electronic conductors (MIECs) have an 
internal electrical short and the ionic species selectively permeates through a dense film of the 
material under a differential partial pressure, as shown conceptually for H2 permeation in Figure 
1 . The potential permeation flux rates of these materials are extremely high. 
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Partially Oxidized 
Hydrocarbon Feed 


Figure 1 : Integration of proton transport membranes (including SEM of membrane) in 
conversion of hydrocarbon fuels to H2. 

Membrane reactors can be used to produce hydrogen from hydrocarbon feedstocks, from natural 
gas to coal to a variety of hydrocarbon containing waste streams such as landfill gas and swamp 
gas. The overall objective of this project is to demonstrate the feasibility of producing hydrogen 
from hydrocarbon based fuels using advanced proton conducting membranes. The objectives of 
this phase are to optimize electron and proton conductivity in SrCe03 membrane materials, 
develop thin film proton conducting membranes on porous supports and to advance our 
fundamental understanding of these materials. The anticipated result is demonstration of high 
hydrogen fluxes through these thin supported membranes. 

1. Fundamental Investigation of Proton Conductors 

1.1. Electrical properties of Eu-doped SrCeCh 

1.1.1. Defect Structure 

Based on internal and external equilibria of the system, the ionic and electronic defect diagram of 
Eu-doped SrCeC>3 was successfully developed as shown in Figure 2. The possibility of interstitial 
defects was neglected because they are energetically unfavorable in the perovskite structure. 
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Figure 2: Ionic and electronic defect diagram in the P 02 -PH 20 plane. 

1.1.2. Effect of Eu-doping level 

The conductivities of SrCei. x Eu x 03 (x = 0.05, 0.1, 0.15, 0.2) were measured in hydrogen using 
impedance spectroscopy. The highest conductivity was achieved with the 1 0 mol% Eu-doped 
SrCeC >3 as shown in Figure 3. 



Figure 3: Conductivity as a function of Eu-dopant level for SrCei_ x Eu x 03 (x = 0.05, 0.1, 0.15, 
0.2)inH 2 . 


19 





UNIVERSITY OF 

FLORIDA 


1.1.3. Effect of oxygen partial pressure 


The conductivity of SrCeo. 9 Euo. 1 O 3 was measured in a hydrogen-nitrogen mixture and in a 
nitrogen -oxygen mixture using impedance spectroscopy. The results are given in Figure 4. 
SrCeo. 9 Euo. 1 O 3 exhibits n-type electronic conductivity in the hydrogen-nitrogen mixture (Figure 
4(a)), while the conductivity of the same material is of p-type when measured in the nitrogen and 
oxygen mixture (Figure 4(b)). The n-type conductivity of SrCeo.9Euo.1O3 in the hydrogen- 
nitrogen mixture is relatively high at temperatures above 700 °C, indicating that this material 
should have good hydrogen permeability at temperatures above 700 °C. 


* 


Figure 4-(a). Total Conductivity of SrCeo.9Euo.1O3 as a function of Ph 2 I/4 
at different temperatures. 



p 
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Figure 4 - (b). Total Conductivity of SrCeo.9Euo.1O3 as a function of P 0 2 1/4 at different 
temperatures. 

1.2. UF proprietary proton conductor 

The transference number of the UF proprietary proton conductor was investigated to determine 
the extent of mixed conduction (protonic and electronic). The UF proprietary proton conductor 
was used as electrolyte in the following hydrogen concentration cell: 

dry H 2 ( P H2 ,I )+Ar, Pt/sample / Pt, dry H 2 ( P Hi , II , higher)+Ar 

The EMF was measured at 500, 550 and 600 °C. The measured voltage is determined by the 
following equation: 


E 


measured 


= t 


H 


RL Ln ME 

IF 


The results are shown in Figure 5. The transference number is about 0.772, 0.785 and 0.788 at 
600, 550 and 500°C, respectively. This shows that it is a proton conductor in this temperature 
range. It is thus a proton conductor that can be used as the solid electrolyte in fuel cells. But the 
electronic conduction is not high enough to make the UF proprietary proton conductor a mixed 
conductor. 
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Figure 5: Proton transference number as function of temperature. 


2. Hydrogen Permeation Membrane Development 

2.1. Membrane supporting tube: materials and fabrication 

2.1.1. NiO-GDC 

The initial concept for a supported membrane focused on co-sintering a nickel oxide and 
gadolinium-doped ceria (NiO-GDC) support and a europium-doped strontium cerate (ESC) film. 
The porous support would result from in-situ reduction of NiO to Ni in hydrogen. For this 
concept to succeed, the ESC must be thermodynamically stable with both the NiO and GDC in 
air, and the ESC and NiO-GDC must have comparable thermal expansion properties in air and 
hydrogen. 

The thermodynamic stability of ESC with both GDC and NiO was investigated in parallel with 
the thermal expansion properties of these materials. Figure 6(a) demonstrates the lack of 
reactivity between ESC and GDC at 1500 °C. Unfortunately, as shown in Figure 6(b), ESC does 
react with NiO above 1400 °C. At 1400 °C, only ESC and NiO are present. While at 1500 °C, 
the powder had completely reacted and the resultant peaks were identified as Ce02 and NiO. 
Thus, based on reactivity, a NiO-GDC support would work if the membrane is sintered < 
1400°C. 

Figure 7 demonstrates the thermal expansion behavior of both ESC and NiO-GDC in (a) air and 
(b) hydrogen. The thermal expansion behavior of the ESC was similar in both air and hydrogen 
(thermal expansion coefficient of 9.4 in air and 9.0 in hydrogen). However, thermal expansion 
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Although the difference in thermal expansion between ESC and NiO-GDC in air can be 
tolerated, the difference in hydrogen results in failure of the membrane-support as shown in 
Figure 8. Figure 8(a) illustrates the surface of a film before hydrogen treatment while Figure 
8(b) shows cracks which have propagated from the support through the film after hydrogen 
treatment. 



Figure 8: Scanning electron micrographs of ESC film (a) before hydrogen treatment and (b) 
after hydrogen treatment 

The results indicate that the NiO-GDC system is not appropriate for an ESC supported 
membrane that is sintered above 1400 °C. 


2.1.2. YSZ 

Another possible system investigated was yttria-stabilized zirconia (YSZ). Unfortunately, the 
YSZ and ESC react at 1500 °C as shown in Figure 9. 


Figure 9: 

1500°C. 
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2.1.3. Membrane supporting tube based on NiO-SC 

Undoped strontium cerate (SC) was investigated as a support material. The main complication 
using SC as the support material is the reaction between SC and NiO at 1450 °C which is below 
the sintering temperature used with the GDC and YSZ supports. Two approaches were taken to 
overcome this problem. First, SC tape was prepared using carbon rather than NiO as the pore 
former. The green strength of the SC-carbon tape was too low to prepare porous support tubes 
via the rolling method. Second, SC tape was prepared with NiO, and temperatures below 
1450°C were investigated. The slurry composition for the SC-NiO tape is shown in Table 1. 

Table 1: Slurry composition of SC containing nickel catalyst for tape casting process 


Process step Material 

Function 

Weight (g) 

SrCe0 3 


52.55 


Powder 


NiO 


22.52 

Mill for 24 hr (l'st stage) Solsperse 

Dispersant 

0.75 

Toluene 


13.11 


Solvent 


E.A 


8.18 

DBP 


3.30 

Add to above and mill for 

Plasticizer 


PFG 

24 hr (2'nd stage) 


0.66 

PVB 

Binder 

4.50 



900 1000 1100 1200 1300 1400 1500 1600 


Temperature (°C) 

Figure 10: Shrinkage rates for the un-doped and europium-doped strontium cerate, prepared by 
uniaxial pressing, and un-doped strontium cerate with nickel oxide, prepared by tape casting, 
according to the sintering temperatures. 
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The shrinkage behavior of uniaxially pressed ESC and SC and tape casted SC-NiO was studied. 
The shrinkage rates for the SC, ESC, and SC-NiO, shown in Figure 10, demonstrate that the ESC 
sinters more easily than SC below 1450 °C and that the ESC and SC-NiO have comparable 
shrinkage rates which are necessary for dense thin film preparation. The limited sinterability of 
SC should improve the porosity of the support tube. Also, the samples with higher europium 
dopant levels demonstrate similar shrinkage rates. 

The change in apparent porosity before and after hydrogen treatment for both 10 wt% and 30 
wt% NiO-SC samples is shown in Figure 1 1. At 1425 °C, the apparent porosity for the 10 wt% 
NiO sample only increased from 5% to 12.5% while for the 30 wt% NiO sample the apparent 
porosity increased from 5% to 25%. Following in situ reduction of NiO to Ni, the SC support 
containing 30 wt% NiO is satisfactorily porous. 

The thermal expansion behaviors of ESC and 30 wt% NiO-SC must be similar in both air and in 
hydrogen to avoid failure of the ESC film since the samples are sintered in air and tested in 
hydrogen. The expansion behavior of ESC is very similar to that of 30 wt% NiO-SC in air as 
shown by Figure 12(a) and the ESC films do not crack after sintering which is confirmed by the 
SEM micrograph shown in Figure 13(a). In hydrogen, the thermal expansion coefficient (TCE) 
of the 30 wt% NiO-SC is approximately 13% higher than the TCE of ESC. As shown by Figure 
13(b), this expansion difference is not significant enough to result in cracking of the ESC film. 

The film shown in Figure 13 was formed using a slurry coating process combined with a 
filtration apparatus. Earlier films that were formed by a conventional slurry coating process 
developed pin-hole defects. The new coating method eliminated the pin-hole defects. 
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Figure 1 1 : Apparent porosity for the sintered NiO-SC support before and after reduction 
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Figure 12: Thermal expansion behavior of ESC and NiO-SC in (a) air and (b) hydrogen. 



Figure 13: Scanning electron micrographs of ESC film on NiO-SC support (a) before and (b) 
after hydrogen treatment. 


The new coating method allowed for preparation of six inch NiO-SC support tubes with a dense 
ESC film. The film remains crack free after treating in hydrogen. 

2.2. Preparation of ESC film on NiO-SC support 

According to our investigations, the final sintering temperature should be below 1450 °C to 
prevent a reaction between NiO and strontium cerate (SC), and above 1400 °C to obtain dense 
films of Eu-doped strontium cerate (ESC). In addition, the pre-sintering temperature for the NiO- 
SC support should be in the range of 1 150 to 1200 °C, in order to obtain compatible shrinkage 
rates between NiO-SC and ESC film at the sintering range. The cracking of the NiO-SC support 
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was observed when the support was pre-sintered at temperatures below 1150 °C. The optimal 
pre-sintering and sintering temperature ranges are given in Figure 10. 


As summarized in Table 2, the preparation of the ESC hydrogen membrane on NiO (or Ni)-SC 
tubular support consists of four steps; and the pictures of a six-inch unit cell thus prepared at 
each preparation step are shown in Figure 14. 


Table 2: Preparation steps of the tubular hydrogen membrane cell and related processes 


Preparation step 

Related process 

Green state tube 

Slurry preparation, tape casting and rolling 
(Tubular type 30wt% NiO in SrCe03 support) 

Partially sintered 

Pre-sintering NiO-SC support at 1150 °C to 1200 °C 

support 
and membrane 

Citrate process/solid state reaction 
(Preparation 10 mol% Eu-doped SrCeOs) 

coating 

Slurry coating ESC on NiO-SC support 
(combined with a filtration apparatus) 

Sintered cell 

Sintering ESC/NiO-SC at 1400 to 1450 °C 

Reduced cell 

Reducing NiO to Ni in support by annealing in H2 



tSiCe 1|(1 Lu t| 0 M uiW qu tfcn inner fiife uf the »iipp»n) 


1 Sintered cell j 



Figure 14: Pictures of a six-inch unit cell of 10 mol% Eu-doped SrCe03 membrane coated on the 
inner side of NiO(or Ni)-SrCe03 tubular support at each processing step. 


2.3. Hydrogen permeation 

By assembling the unit cell with gas inlet and outlet systems, we obtained a hydrogen permeation 
reactor, the picture of which is shown in Figure 15. Installing the hydrogen permeation reactor in 
a furnace and connecting the reactant and sweep gases, we are now ready for testing the 
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hydrogen permeability of the six-inch unit cell. Mass flow controllers are used to control the gas 
flow rate, and hydrogen permeation rate is measured by a mass spectrometer. The actual 
experimental set-up is shown in Figure 16. 



Figure 15: The hydrogen permeation reactor with the six-inch unit cell. 



Figure 16: Experimental set-up for hydrogen permeation test. 
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Figure 17: Hydrogen flux and water vapor formation at 900 °C. Corresponding experimental 
conditions are given in the figures. 

Gas flow rates were controlled by mass flow controllers. One side of the membrane cell was 
exposed to H 2 (99.999%) diluted to the desired concentration using He (99.999%) and/or Ar 
(99.999%) at 20 cm 2 3 /min total flow rate (feed gas). The other side was swept with He at 20 
cm 3 /min (sweep gas). Hydrogen permeation was measured using a mass spectrometer. Figure 17 
shows hydrogen fluxes obtained from the hydrogen permeation reactor at 900 °C under different 
experimental conditions. It is obvious that the hydrogen flux is sensitive to the partial pressure of 
hydrogen in the feed gas. The water vapor fluxes are also shown in the figures. 

The result shows that at 900 °C a single tube of our current generation of membrane is capable of 
producing about 5 cc/min of pure H 2 . The hydrogen flux can be significantly increased by 
optimizing the membrane, and proportionally increased by increasing the number of tubes. 

3. Publications and Presentations 

3.1. Publications 

1. “Hydrogen Permeability and Microstructure Effect on Mixed Protonic-Electronic Conducting Eu- 
doped Strontium Cerate”, S.-J. Song, E.D. Wachsman, J.M. Rhodes, H.-S. Yoon, G. Zhang, K.- 
H. Lee, S.E. Dorris and U. Balachandran, Journal of Materials Science, accepted. 

2. “Hydrogen Production from Fossil Fuels with High Temperature Ion Conducting 
Ceramics”, E. D. Wachsman and M.C. Williams, Interface, 13-3, 32-37 (2004). 

3. “Hydrogen Permeability of SrCeo.gsEuo.osOa-g (x = 0.05, M = Eu, Sm)”, S. J. Song, E. D. 
Wachsman, J. Rhodes, S. E. Dorris and U. Balachandran, Solid State Ionics, 167, 99-105 
(2004). 
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4. “Numerical Modeling of Hydrogen Permeation in Chemical Potential Gradients”, S. J. 

Song, E. D. Wachsman, J. Rhodes, S. E. Dorris and U. Balachandran, Solid State Ionics, 164 , 
107-116(2003). 

5. “Defect Structure and n-Type Electrical Properties of SrCeo.gsEuo.osCW’, S. J. Song, E. D. 
Wachsman, S. E. Dorris and U. Balachandran, Journal of the Electrochemical Society, 150 , 
A1484 (2003). 

6. “Electrical Properties of p-Type Electronic Defects in the Protonic Conductor 
SrCeo. 95 Euo.o 503 -d”, S. J. Song, E. D. Wachsman, S. E. Dorris and U. Balachandran, Journal of 
the Electrochemical Society, 150 , A790 (2003). 
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470 (2003). 
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SENSORS FOR HYDROGEN LEAK DETECTION 

Faculty Lead: Dr. Mark Law and Dr. Jenshan Lin, Electrical & Computer 

Engineering 

OVERVIEW 

In this FY03 program, we have successfully demonstrated novel hydrogen sensors including 
wide energy band-gap based sensors and biological assays sensors. We also demonstrated the 
feasibility of remote detection using wireless transceiver and microcontroller. By using low 
power solid state sensors, the integrated wireless hydrogen sensor can be made compact and 
lightweight to facilitate a fast deployment of a hydrogen leak detection network. 

In the beginning of the project period, our first prototype sensor needs to operate at elevated 
temperature. Through the improvement in device structure, we achieved room-temperature 
(25°C) sensing using AlGaN/GaN MOS-HEMT sensor. Toward the end of this project period, 
we developed the novel ZnO nanorod sensor, which has much lower power consumption yet 
better sensitivity. This paves the way for the development of self-powered wireless hydrogen 
nano-sensor to be developed in the next phase. 

PROJECT TITLE: WIDE ENERGY BAND-GAP BASED GAS SENSORS 

Task PI: Dr. Fan Ren, Department of Chemical Engineering 

Co-I: Dr. Steve J. Pearton, Department of Materials Science and Engineering 

Collaborator: Dr. Jenshan Lin, Department of Electrical and Computer Engineering 

PROJECT GOALS 

Demonstrate room temperature hydrogen sensor with wide energy bandgap semiconductor 
devices. 

ACCOMPLISHMENTS 

SUMMARY 

This report describes the development of solid-state hydrogen sensors using GaN-based devices 
and ZnO-based devices. Device structures and measured performances of hydrogen detection 
using GaN Schottky diode, AlGaN/GaN Schottky diode, AlGaN/GaN HEMT, Pt-ZnO Schottky 
diode, and ZnO nanorod are compared. These devices, operating with low power consumption, 
are suitable for integration with electronic circuits to form a low power wireless hydrogen sensor 
module. A prototype of wireless sensor system is also described. 

BACKGROUND 

Sensors for hydrogen leak detection are crucial to ensure safety in many applications that use 
hydrogen as fuel, e.g., space exploration missions in NASA, hydrogen-powered vehicles in 
automobile industry. Semiconductor gas sensors that have advantages of small size and 
reusability are promising technologies for hydrogen leak detection. The solid-state sensors, 
which respond to chemical detection with an electrical signal, can also be easily integrated with 
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solid-state electronic circuits that process and transmit the sensor data. If the sensor data is 
transmitted using radio waves, a wireless sensor network can be designed to link multiple 
sensors together to a central monitoring station. Research efforts are currently carried out in 
different disciplines to achieve this goal. Nevertheless, the fundamental chemical sensing 
mechanism occurred in the sensor device plays the most important role, and the search for the 
optimum hydrogen sensor has been an interesting research subject recently. Pioneering research 
in the past few years have produced promising results. The silicon-based approach developed in 
NASA used Schottky diodes as the sensors in a hydrogen detection system for various space 
exploration missions. The approach can also integrate novel technologies such as MEMS, 
nanomaterials, and SiC [1], Through this program, we demonstrated several different types of 
hydrogen sensors using GaN-based diodes and transistors [2], In this report, the development of 
GaN-based hydrogen sensors using diodes and MOS-HEMT is reviewed and their performances 
are compared. As nanotechnology getting matured, nano-scale devices were proposed for gas 
sensing. In the end of this project period, our research effort shifted to using ZnO 
nanorods/nanowires as hydrogen sensor. Regardless of the type of the sensor, as long as the 
sense of hydrogen results in a change in its electrical characteristics, the sensor device can be 
integrated with electronics to form a complete sensor system and network. The integration with 
electronic circuits and a prototype wireless sensor system are also described in this report. 

GaN-BASED HYDROGEN SENSOR 

GaN Schottky Diode Sensor 

Figure 1 shows the microphotographs of the GaN Schottky diode hydrogen sensor and the test 
fixture where the substrate was mounted on for testing in chamber. Figure 2 shows the cross- 
sectional view of the Schottky diode structure. Front-side ohmic contacts of Ti/Al/Pt/Au were 
formed by lift-off and subsequent annealing at 600°C. The Schottky contacts were formed by 
lift-off of e-beam deposited Pd or Pt, 24nm and 1 5nm thick respectively. The contact diameter is 
80pm in both cases. The devices were wire-bonded to a test fixture using Ti/Au bond-pads and 
Au wires for contact [3]. 



Figure l.GaN Schottky diode hydrogen sensor (left) and the substrate mounted on a test fixture 

used for testing in chamber (right). 
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Figure 2. Cross-sectional view of the GaN Schottky diode structure. 

Figure 3 shows the forward current-voltage characteristics of the Pd/GaN diode at 170°C in pure 
N 2 and 1 0% H 2 in N 2 ambients. The addition of 1 0% H 2 to the ambient produces a shift of 74mV 
at a forward current of 2mA. Note that at higher voltage and current bias, the signal change is 
even more discemable, e.g., 1mA at a forward bias of IV or lOOmV at a bias current of 6mA. 
The sensing mechanism in semiconductor gas sensors is thought to be creation of a polarized 
layer in the semiconductor surface by hydrogen atoms diffusing through the metal contact [4,5]. 
The adsorbed hydrogen is then assumed to change the work function of the metal [5]. Recent 
elastic recoil detection measurements on Pt-GaN Schottky diodes confirm the presence of 
hydrogen at the metal semiconductor interface after exposure to molecular hydrogen at 400°C 
[ 6 ]. 



Figure 3. Current- voltage characteristics of the GaN Schottky diode hydrogen sensor at 170°C in 
pure N 2 and 10% H 2 in N 2 ambients. 
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Figure 4. Time response of the GaN Schottky diode hydrogen sensor at 150°C 

The time response of the GaN Schottky diode hydrogen sensor is demonstrated in Figure 4. In 
this case, Pd/GaN diodes were biased at a constant forward bias of 0.5V at 150°C and the 
resulting forward current measured as a function of time as the ambient was cycled from N 2 to 
10% H 2 in N 2 and then back to N 2 . The forward bias current increase about 0.3mA as H 2 is 
introduced and the subsequent decay back towards the initial current as N 2 is introduced into the 
chamber. The increase in forward current is very rapid after switching in the 10% H 2 in N 2 gas. 
From the observation of this rapid rise in current, we believe that diffusion of hydrogen through 
the metal contact is not the limiting factor in the time response of the diodes but rather the mass 
transport of gas into the enclosure. Similarly, the longer recovery times are at least partially 
caused by the time needed to completely flush the H 2 out of the measurement enclosure. 

Schottky Diode Sensor on AlGaN/GaN HEMT Structure 

A variation of Schottky diode sensor was fabricated on AlGaN/GaN High Electron Mobility 
Transistor (HEMT) layer structure. Figure 5 shows the microphotographs of the AlGaN/GaN 
MOS diode hydrogen sensor and the test fixture where the substrate was mounted on for testing 
in chamber. Figure 6 shows the cross-sectional view of the AlGaN/GaN Schottky diode 
structure. 

This diode structure also has the capability of detecting hydrogen. Figure 7 shows the time 
response of the AlGaN/GaN Schottky diode sensor when ambient conditions of pure N 2 and 10% 
H 2 in N 2 were cycled through the test chamber. The diode was biased at 2V with a bias current of 
about 10mA. The introduction of hydrogen into the chamber caused a current increase of more 
than 0.4mA. The measured measurement result also demonstrated the reversibility of this diode 
sensor. 
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Figure 5. AlGaN/GaN Schottky diode hydrogen sensor (left) and the substrate mounted on a test 
fixture used for testing in chamber (right). 
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Figure 6. Cross-sectional view of the AlGaN/GaN Schottky diode structure. 



Figure 7. Time response of the AlGaN/GaN Schottky diode sensor. 
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AlGaN/GaN MOS HEMT Sensor 


In addition to Schottky diodes, the AlGaN/GaN HEMT structure can also be used to make a 
sensor based on three-terminal transistor topology. In this case, drain to source current is 
measured, which is modulated by the gate that is affected by the hydrogen concentration in the 
ambient. A MOS-HEMT sensor was designed and fabricated [Error! Bookmark not defined.]. 
Figure 8 shows the micrograph of the sensor. The structure is the same as other FET-type 
transistors. SC 2 O 3 was used as the gate oxide. The DC I-V measurement shown in Figure 9 
confirms its transistor characteristics and the response to hydrogen. As can be seen in the figure, 
this sensor has largest current change in response to hydrogen when biased in the range of 1 0 - 
20mA. 



Figure 8 . AlGaN/GaN MOS HEMT sensor. 
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Figure 9. I-V curves of the AlGaN/GaN MOS-HEMT sensor. 
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Figure 10 shows the time response of the AlGaN/GaN MOS-HEMT sensor when ambient 
conditions of pure N 2 and 1% H 2 in N 2 were cycled through the test chamber. The gate was 
biased at -3V with a bias current of about 16.25mA in pure N 2 . As 1% H 2 was introduced into 
the chamber, the current increased to 16.7mA. This 0.45mA increase is about the same as 
observed in the case of AlGaN/GaN diode sensor fabricated on the same layer structure, but with 
only 1% H 2 . The factor of 10 improvement is due to the gain of the HEMT transistor. The result 
clearly shows the advantage of using AlGaN/GaN MOS-HEMT as hydrogen sensor as compared 
to Schottky diodes. The sharp response also shows the fast response time of less than one second 
and the reusability. In addition, since the sensor device uses the same HEMT layer structure that 
can be used to build radio frequency transceiver circuits, potential monolithic integration of 
sensor device and wireless transceiver is feasible. Noted that GaN-based HEMT is expected to 
be the first GaN-based electronic device to be commercialized for next generation radar and 
wireless communication systems. Since GaN-based devices can operate at elevated temperatures, 
the GaN-integrated wireless sensor can be used in harsh environment where sensor devices 
fabricated on other semiconductor technologies cannot survive. 

ZnO-BASED SENSOR 

Pt-ZnO Schottky Diode Sensor 

We have been investigating ZnO for sensor application since it is a stable and bio-safe material 
(bandgap=3.3eV). A Pt-ZnO Schottky diode was fabricated, as shown in Figure 11. 
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Figure 1 1 . Pt-ZnO Schottky diode photo (left) and structure (right). 



Figure 12. Current-voltage characteristics of Pt-ZnO Schottky diode sensor in response to 
different levels of hydrogen concentration. 

The I-V characteristic of the Pt-ZnO diode is shown in Figure 12. The I-V curve changes in 
response to different levels of hydrogen concentration. In pure N 2 , the characteristic is similar to 
other diodes with exponential increase in current at forward bias. As hydrogen concentration 
increases, the current at reverse bias increases and the forward bias turn on characteristics 
gradually disappears. This interesting behavior shows that as hydrogen concentration increases, 
the device behaves more like a resistor than a diode. At 1 % H 2 level, the device is essentially a 
resistor of 83Q. Compared to GaN Schottky diode, this Pt-ZnO is more sensitive and can detect 
ppm level hydrogen concentration change. 

ZnO Nanorod Sensor 

ZnO can also be used to fabricate nanoscale devices like nanowires and nanorods [7]. A 
depletion-mode Field-Effect Transistor (FET) was demonstrated using a ZnO nanorod [8], as 
shown in Figure 13. The device was found to be sensitive to UV with significant changes in DC 
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I-V characteristics (current increased by larger than 4 times). A comparison of I-V characteristics 
measured at room temperature is shown in Figure 14. 



— 1pm 

Figure 13. ZnO nanorod FET. 



Figure 14. I-V characteristics of ZnO nanorod FET measured at room temperature in the dark 

(left) and under UV illumination (right). 

The result indicates that ZnO has a very good potential for low power operation since the current 
is very low. By integrating with low power electronic circuits, a long lifetime battery-powered 
wireless sensor can be achieved. With energy harvesting devices, a self-powered wireless sensor 
is also possible. 

The same ZnO nanorod can be used for hydrogen detection as well. Figure 15 shows the photo of 
a ZnO nanorod hydrogen sensor. Compared to Figure 13, the gate metal is removed. The 
nanorod channel resistance is modulated by the hydrogen concentration. Figure 16 shows the 
measured I-V characteristics of ZnO nanorod sensor under pure N 2 and 10% H 2 in N 2 , at 
different temperatures. The existence of H 2 changed the resistance of ZnO nanorod, and the 
change is larger at higher temperature. Compared to ZnO Schottky diode sensor, the current 
change is smaller, which suggests that an array of multiple ZnO nanorods is needed to increase 
the sensitivity. 
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Figure 15. ZnO nanorod sensor with length=3.5pm and diameter=130nm. 



Figure 16. Measured I-V characteristics of ZnO nanorod sensor. 


WIRELESS SENSOR SYSTEM PROTOTYPE 

Integration with Wireless Transceiver 

There are strong reasons for developing distributed arrays of integrated semiconductor gas 
sensors for detecting hydrogen leaks and being able to transmit data on the location of the leaks 
back to a central monitoring station or control facility. In addition, these detectors would have 
dual-use in automobiles and aircraft, fire detectors, and exhaust diagnosis from industrial 
processes. 

By using low noise amplifier, microprocessor, and wireless transceiver, we have successfully 
demonstrated a low-cost wireless sensor system prototype [3]. Figure 17 shows the remote 
wireless transmitter module and the monitoring receiver module. The size of the prototype 
module is smaller than palm-size and can be easily integrated in a wireless sensor network. 
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Figure 17. Photos of the remote wireless transmitter module (left) and the monitoring 
receiver module (right). 

Operation of Wireless Hydrogen Sensor System Prototype 

The detailed block diagram of the modules is shown in Figure 18. The remote wireless 
transmitter module has an input terminal that can detect voltage-current change from any sensor 
with electrical signal output. In the case of hydrogen sensor, voltage sensing is used. If the sensor 
is small enough, it can be integrated with the module. A low noise amplifier amplifies the sensor 
signal for digital signal processing in the micro-controller. The micro-controller has built-in 
analog-digital converter and encoder for secured wireless transmission. After the signal 
processing, the sensor data is sent by a transceiver and an antenna operating at 916MHz. The 
antenna was designed as a low-cost low-profile printed loop antenna fabricated on the circuit 
board inside the module. 



Figure 18. Prototype wireless hydrogen sensor system. Left hand side of the figure shows the 
wireless transmitter module and its block diagram. Right hand side of the figure 
shows the monitoring receiver module and its block diagram. Small loop antennas 
are integrated inside the module 

On the receiver side, same antenna and transceiver are used as radio frequency interface. After 
the signal is received and down-converted to low frequency, the same micro-controller is used to 
process and decode the sensor data. The sensor data is then displayed on the LCD screen. The 
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voltage reading on the screen is the same voltage detected at the remote sensor, and the change 
would indicate the detection of hydrogen leakage. A more complicated system can be built upon 
this platform. For example, a programmed threshold level to trigger the alarm and shut off the 
gas valve. With various multiple access techniques, e.g., code-division multiple access (CDMA), 
a wireless hydrogen sensor network can be realized. 
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PROJECT TITLE: SENSING HYDROGEN BY USING BIOLOGICAL 
ASSAYS 

Task PI: Dr. Z. Hugh Fan, Mechanical & Aerospace Engineering 
Collaborator: Dr. Jenshan Lin, Electrical & Computer Engineering 

PROJECT GOALS 

The goal of this project is to develop a novel hydrogen sensor using an enzyme-catalyzed 
reaction and microfluidic technology. The anticipated benefits over the state-of-the-art 
palladium- or its alloy-based hydrogen sensors include the ability to operate at ambient 
temperatures, the ability to operate in background gases due to the specificity of the enzyme, fast 
response time, and no recovery time. 

ACCOMPLISHMENTS 

Concept. The principle of using a biological assay for detecting hydrogen is shown as follows: 

H 2 + 2BV^ Wo*™* , 2H . + W r 

The enzyme, hydrogenase, catalyzes the oxidation of H 2 ; BV 2+ (benzyl viologen) acts as the 
electron acceptor, producing BV + that can be detected electrochemically. The enzyme was the 
cytoplasmic hydrogenase from Ralstonia eutropha, supplied by a collaborating research group in 
Germany. 

Experimental Setup. The H 2 testing station was assembled as shown in Figure 1 . It consists of 
a H 2 tank, a flow cell, an electrochemical detector, and a data acquisition system. Varying 
concentrations of H 2 gas (mixed with N 2 ) flows through the electrochemical flow cell. The cell 
contains 1 mL of an aqueous solution with the pressure at ~1 atm and at ambient temperature. 
The electrochemical detection is achieved using a potentiostat that controls the cell voltage and 
measures the electrical current. 
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Figure 1. The H 2 testing station consists of a H 2 tank, a flow cell, an electrochemical 
detector, and a data acquisition system. An exploded view of the electrochemical flow 
cell is shown on the left. 


Electrochemical Detection. When a H 2 gas sample flowed through the detection system, H 2 
was collected by reacting with a reagent mixture containing hydrogenase. The resultant product, 
BV + , was measured using chronoamperometry technique and the electrochemical current was 
used as the detection signal. An example showing the detection of H 2 is in Figure 2. 



Time (s) 

Figure 2. Hydrogen detection using chronoamperometry. Four repeats were performed 
at each H 2 concentration (diluted in N 2 ). Background is 100% N 2 . 

Calibration Curve. After optimizing the experimental protocol to achieve repeatable detection 
at low H 2 concentration, we obtained the calibration curve shown in Figure 3. A linear response 
was observed from 1 to 100% H 2 , and electrochemical signal is normalized against enzyme 
reaction time. A calibration curve was also developed from gas mixtures containing 20% O 2 , 
which is similar to atmospheric air. 
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Figure 3. Calibration curve for current (normalized by enzyme reaction time) and 
percentage of H 2 in gas streams mixed with N 2 . The current is obtained by 
chronoamperometric detection of BV + , which is the reaction product of BV 2+ and H 2 
catalyzed by hydrogenase. 

Work in Progress. We are implementing our detection protocol on a microfluidic platform. 
Expected improvements by miniaturization include (a) more efficient use of expensive reagents 
(i.e. enzymes), (b) increased sensitivity, and (c) faster response time. We have fabricated several 
microfluidic designs on a plastic chip, in which we have been able to accomplish simultaneous 
gas-liquid flow. We have also been able to demonstrate and calibrate electrochemical detection 
of BV 2+ in the micro fluidic platform. 

RELATED ACTIVITIES 

• On Mar. 15, 2005, we hosted a lab tour for Dr. Gary Hunter from NASA Glen Research 
Center while he visited UF campus. 

• During the NASA’s review meeting on May 10-11, 2005, we hosted a lab tour for 7 scientists 
from NASA Glen Research Center and Kennedy Space Center, as well as a separate tour for 
other scientists (Students and Faculties) from Florida Solar Energy Center, Florida State 
University, University of Center Florida, and University of South Florida. 

• A paper on this project has been published in Analytical Chemistry, Vol. 77, 4969-4975, 
2005 . It is entitled “Hydrogen Sensing by Enzyme-Catalyzed Electrochemical Detection”. 

• One provisional patent application has been filed in this reporting period. It is entitled 
“Hydrogen Sensor Using Enzyme-Catalyzed Reaction”, filed on March 16, 2005, with 
application No. 60/662,504. 

• One related proposal to National Science Foundation (NSF) has been recommended for 
funding. 
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INSTRUMENTATION FOR HYDROGEN LEAK DETECTION 
Faculty Lead: Dr. David Hahn, Mechanical & Aerospace Engineering 

PROJECT TITLE: RAMAN SPECTROSCOPY FOR REMOTE 
DETECTION OF HYDROGEN LEAKS 

Task PI: Dr. David W. Hahn, Department of Mechanical and Aerospace Engineering 

PROJECT GOALS 

Leak detection of rocket propellants, notably hydrogen, is an ongoing area of concern for the 
space shuttle program. To aid in safety and efficiency, it is desirable to have remote leak 
detection capabilities that will enable standoff interrogation of various hydrogen handling 
facilities, such as transfer lines and piping, as well as shuttle launch pad detection. This research 
involves the development of diagnostic tools for portable, remote hydrogen leak detection using 
laser-based techniques. Laser-based approaches will include Rayleigh/Mie light scattering and 
Raman spectroscopy. The proposed research will focus on Rayleigh/Mie/Raman integration 
using temporal and spectral data analysis in combination with pulsed laser excitation to 
discriminate from background contaminants and to optimize leak detection at the leak source. 
The overall project goal is the development of a field-portable sensor package for field 
evaluation. 

ACCOMPLISHMENTS 
I. PROJECT SUMMARY 


Leak detection of rocket propellants, notably hydrogen, is an ongoing area of concern for the 
space shuttle program. To aid in safety and efficiency, it is desirable to have remote leak 
detection capabilities that will enable standoff interrogation of various hydrogen handling 
facilities, such as transfer lines and piping, as well as shuttle launch pad detection. This research 
involves the development of diagnostic tools for portable, remote hydrogen leak detection using 
laser-based techniques. A laser-based approach will utilize Raman spectroscopy to detect 
diatomic hydrogen in the gaseous state. The research has focused on Raman spectroscopy using 
temporal and spectral data analysis in combination with pulsed laser excitation to discriminate 
from background contaminants and to optimize leak detection at the leak source. The overall 
project goal is the development, testing, field evaluation, and optimization of a field-portable 
sensor package. 

II. PROJECT OVERVIEW 


Research to date has focused on the development of diagnostic tools for portable and remote 
hydrogen leak detection using laser-based approaches, namely Raman spectroscopy. The 
research has focused on Raman spectroscopy integration using temporal, statistical, and spectral 
data analysis in combination with pulsed laser excitation to discriminate from background 
contaminants and to optimize leak detection at the leak source. 
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The technique of Raman spectroscopy is based on the inelastic scattering of electromagnetic 
radiation due to interactions of the incident radiation with the vibrational energy states of the 
scattering species. The nature of Raman spectroscopy makes selection of the excitation laser 
important. Notably, the process of Raman scattering is a relatively weak one as compared to 
elastic scattering. For example, the Raman scattering cross-section of hydrogen (H 2 ) is 9x1 0~ 30 
cm 2 /sr, which is nearly two orders of magnitude less than the Rayleigh scattering cross-section of 
hydrogen. However, the Raman scattering cross-section scales nominally as the inverse of the 
fourth power of excitation wavelength. It is therefore desirable to use as short a wavelength as 
possible for excitation, although consideration must also be given to fluorescence, which can 
interfere with the Raman-shifted hydrogen signal. Research led to the selection of a frequency- 
tripled (355 nm) Nd:YAG laser as the appropriate laser source. 

III. SUMMARY OF PROGRESS 


1) Optical Configuration : Using the laboratory hydrogen test cell developed last program year, 
tests were performed to determine the optimal excitation wavelength. Using the 6-inch path 
length cell, the hydrogen Raman signal was quantified over a range of hydrogen concentrations. 
All hydrogen concentrations were prepared in nitrogen to avoid the possibility of any explosions 
initiated by the laser pulse. As noted earlier, the response for 355-nm excitation revealed a 16- 
fold increase in response slope, hence all recent work was focused on the use of 355-nm 
excitation. Detailed measurements of the hydrogen Raman signal were recorded as a function of 
hydrogen concentration for two separate laser pulse energies (125 and 250 mJ). In addition, 
Raman data was recorded for two different size collection lens, namely a 2” lens and a 4” lens. 
These data enable additional investigation of the scaling of the Raman signal with laser pulse 
energy and solid collection angle. 

The Raman signal is expected to be directly proportional to the laser pulse energy and the solid 
angle of collection. This latter quantity should scale as the square of the collection lens diameter. 
To quantify this scaling, the Raman signal is presented in Figure 1 as a function of the laser pulse 
energy times the square of the collection lens for a fixed hydrogen concentration of 12.5%. As 
observed in Figure 1, the linear response validates the theoretical scaling, enabling use of the 
data collected to extrapolate to different laser pulse energies and collection optics. This data was 
useful for the design and testing of the prototype field instrument. 
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Figure 1 . Raman response vs. the laser pulse energy times the collection solid angle. 

2) Photon Arrival Rates : The quantification of photon arrival rates for Raman detection of 
hydrogen gas can be predicted based on the optical configuration and the Raman cross-sections 
for hydrogen. The geometry for Raman detection of hydrogen is shown below in Figure 2. 
Theoretically, the photon arrival rate may be expressed as 

S Raman = F 0 a’ N (hv)' 1 dV dfi, 

where N is the hydrogen molecular density, o’ is the Raman scattering cross-section for 
hydrogen at the desired wavelength, hv is the photon energy, dV is the excitation or scattering 
volume, and dQ is the solid collection angle. The Raman scattering cross-section for hydrogen 
is l.OxlO' 29 cm 2 /sr at 355 nm. Using a laser pulse energy of 125 mJ/pulse, a 4” diameter 
collection lens, and a separation distance of 8 ft., and 1% hydrogen in the 6” cell, the calculated 
photon arrival rate is 2x1 0 4 photons per pulse. Assuming typical fiber optic and spectrometer 
losses (20 pm slit size), the estimated photon rate at the photomultiplier tube was about 500 
photons/pulse. 

Raman detection limits were calculated experimentally using the slopes of the Raman signal 
response function, based on laboratory measurements. Specifically, using a 125 mJ/pulse laser, 
which is consistent with the new field-instrument laser, and a 4” collection lens, the Raman 
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signal as a function of hydrogen concentration was evaluated. The data were collected at a 
separation distance (gas cell to collection optics) of 8 feet. Using the slope to extrapolate the 
detection limit at three times the noise limits (3a), the detection limit was calculated at 1% 
hydrogen. 


N (molecules/cm 3 ) 

J o' (cm 2 /sr) 

F 0 (mJ/cm 2 ) 

dQ (sr) 

dV (cm 3 ) 



Figure 2. Collection geometry for Raman photon flux calculations. 

3) Hardware for Field Instrument '. A key component of the recent work was design and 
construction of the remote Raman instrument that is suitable for field evolution. Field 
measurements provide the range of conditions, such as ambient aerosol loadings, varying wind, 
temperature and humidity, which are not readily created in a laboratory setting. The prototype 
instrument was assembled using the following design and components, as illustrated in Figure 3. 
The laser was a BigSky CFR 200 with frequency-tripling package that provides 95 mJ/pulse of 
355-nm laser light. The maximum repetition rate is 10 Hz, and the approximate full-width half- 
maximum pulse width is 10 ns. A rack mountable power supply was purchased. The light is 
collected with a bi-furcated fiber optic and input to two matching custom-designed (Ocean 
Optics) monochrometers, with the output directed to two matched photomultiplier tubes. The 
system was constructed during December 2004. Photographs of the system, including the 
overall optical layout, monochrometers, and system configuration are shown in the following 
figures. 
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Figure 3. Schematic of the field instrument for Raman detection of hydrogen. 



Figure 4. Custom monochrometers for the Raman instrument. 







Figure 5. Photograph of the Raman instrument, showing the optical paths of the output beam 

(blue) and the collected signal (green). 



Figure 6. Photographs of the prototype detector and instrument rack (which includes laser power 

supply). 
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4) Instrument Assessment : Preliminary measurements were conducted in the laboratory using the 
field instrument. No significant hydrogen signal was detected with the configuration under the 
laboratory conditions. This was attributed to the lack of sufficient signal throughput through the 
custom-made monochrometer modules. Because the modules were custom designed, it was fully 
anticipated that adjustments would have to be made to optimize the tradeoff between throughput 
and resolution. The original lO-pm slits were determined to be too narrow. Ocean Optics has 
provided additional entrance and exit slits at no cost to enable optimization of the Raman signal. 
The redirection of the program funding did not enable complete testing of the optimized field 
instrument. It is expected that the current instrument is capable of operating in a field 
environment, although it is difficult to extrapolate laboratory results to the field conditions. For 
example, the laboratory work utilized sealed optical chambers for all work, as the direct venting 
of hydrogen was deemed a safe practice. The test cell windows are known to be a significant 
source of spectral noise; hence moving to an open-air outdoor situation might increase 
instrument performance. Alternatively, the ambient environment offers the potential for other 
sources of spectral noise, such as fluorescence from organic compounds, which may degrade 
overall instrument detection limits. Only a detailed field analysis will answer such questions. 

In summary, the current program demonstrated that Raman spectroscopy utilizing time-resolved, 
pulsed-laser excitation is a viable technology for hydrogen lead detection. Furthermore, current 
off-the-shelf components (lasers, detectors, and signal processors) are available to implement 
this technology, as demonstrated during this project. 
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CRYOGENIC TRANSPORT & STORAGE (CTS) 

Faculty Lead: Dr. James Klausner, Mechanical & Aerospace Engineering 

PROJECT TITLE: TERRESTRIAL CRYOGENIC TWO-PHASE FLOW & 
HEAT TRANSFER 

Task PI: Dr. James F. Klausner, Department of Mechanical & Aerospace Engineering 
Co-I: Dr. Renwei Mei, Department of Mechanical & Aerospace Engineering 

PROJECT GOALS 

Cryogenic fluids flowing through pipelines often operate in a two-phase flow regime. In order to 
design transport systems for two-phase flows it is critical to be able to predict the pressure 
gradients and heat transfer coefficients associated with the fluid and thermal conditions. The 
most reliable two-phase flow models for predicting pressure gradient and heat transfer 
coefficients in pipe flow are specific to a particular flow configuration or flow regime, i.e. 
bubbly, stratified, annular, slug etc. . . Extensive flow regime pattern recognition experiments 
have been carried through with air/water, oil/gas, and refrigerant two-phase flows. From these 
experiments flow regime maps have been constructed. However, the general validity of such 
maps has not yet been demonstrated. Very sparse information on two-phase flow regimes with 
cryogenic fluids is available. Thus, designers of cyrogenic systems rely on generic flow regime 
maps. There exists a significant need to collect flow regime data for the flow of cryogenic fluids 
for a variety of operating conditions and gravitational environments. It is anticipated that the 
data collected from this investigation will be used to explore the applicability of existing flow 
regime maps to cryogenic fluids. The results will provide designers of cryogenic systems the 
means to identify the operating flow regimes. 

In addition, hundreds of two-phase flow pressure gradient and heat transfer models have been 
proposed in the literature. The reason so many different models have been proposed is that each 
is only efficacious over a narrow parameter space. When new data are available that are not well 
predicted from existing models, there is a tendency to develop a new model. There exists very 
little experimental data concerning two-phase pressure gradients and heat transfer for the flow of 
cryogenic fluids, especially in the low turbulent Reynolds number range with low heat fluxes. 
There exists a significant need to compile a data based on cryogenic two-phase pressure gradient 
and heat transfer. There also exists a need to develop mechanistically based pressure gradient 
and heat transfer models that have general validity over a wide parameter space, which will be 
carried out as part of the experimental investigation. It is anticipated that these data will allow 
the development of pressure gradient and heat transfer models that may be used by designers of 
cryogenic thermal transport equipment. 

ACCOMPLISHMENTS 

ABSTRACT 

Cryogenic fluids have found many practical applications in today’s world, from cooling 
superconducting magnets to fueling launch vehicles. In many of these applications the cryogenic 
fluid is initially introduced into piping systems that are in excess of 150 degrees Kelvin higher 
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than the fluid. This leads to voracious evaporation of the fluid and significant pressure 
fluctuations, which is accompanied by thermal contraction of system components. This process 
is known as chilldown, and although it was first investigated more than 4 decades ago, very little 
data are available on the momentum and energy transport during this transient process. 
Consequently, the development of predictive models for the pressure drop and heat transfer 
coefficient has been hampered. In order to address this deficiency, an experimental facility has 
been constructed that enables the flow structure to be observed while temperatures and pressures 
at various locations are measured. This study focuses on the inverse numerical procedure used to 
extract the transient heat transfer coefficient information from the data collected; this information 
is then used to evaluate the performance of various correlations for heat transfer coefficient in 
the flow boiling regime. The method developed utilizes flow structure information and 
temperature measurements, in conjunction with numerical computations for the temperature field 
within the tube wall, to calculate the heat transfer coefficient. This approach allows the 
transition point between the film boiling regime and the nucleate boiling regime to be 
determined, and it also elucidates the variation of the heat transfer coefficient along the 
circumference of a horizontal tube, with the heat transfer on upper portion being significantly 
smaller than that at the bottom. 

INTRODUCTION 

Cryogenic fluids have been used in various applications throughout the years. One application 
that has historically been of great interest is the use of cryogenic fluids for rocket propulsion. 

This interest has been sparked by the fact that “cryogenic propellants are more energetic and 
environmentally friendly than current storable propellant” [1] and the storage systems for these 
cryogenic propellants “have an advantage in reduced weight compared to super-critical tanks” 

[1]. As the propellant tanks are filled, the cryogenic fluid is introduced into a transfer line that is 
in thermal equilibrium with environment, which results in voracious boiling within the line; this 
phenomenon is referred to as line “chilldown” or line “cool-down”. Chilldown is characterized 
by large temperature differences, rapid transients and pressure fluctuations. 

The phenomenon of chilldown is of interest since it directly impacts the design of delivery 
systems for the propellant. For example the magnitude of the pressure oscillation determines the 
thickness of the material used for the transfer lines and the heat transfer rate determines the type 
and thickness of insulation that is used. Hence a proper understanding of the process allows for a 
more economical design of cryogenic propellant delivery systems. 

The vast majority of the studies carried out on chilldown have focused on developing models to 
predict the time taken to complete chilldown, as well as, predict the temperature and pressure 
during the process [2-6]. These studies utilize correlations for the heat transfer coefficient that 
were developed for either steady flows or non-cryogenic flows; neither of which applies to the 
chilldown process and as a result only fair agreement is obtained when compared to experimental 
results. 

The lack of correlations for heat transfer coefficient during chilldown stems from the fact that 
directly measuring the heat transfer coefficient during cooldown is difficult. This study focuses 
on an inverse heat conduction technique that enables the transient heat transfer coefficient during 
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chilldown to be determined. An important feature of this technique is that it accounts for the 
flow structure in calculating the heat transfer coefficient. 

EXPERIMENTAL FACILITY 

The experimental facility constructed to investigate the chilldown process is schematically 
illustrated in Fig. 1. The liquid nitrogen supply is stored in a 1580 kPa (230-psi) dewar. The 
driving potential for the nitrogen flow through the facility is supplied by the dewar pressure. 
Upon exiting the dewar, the nitrogen is subcooled by a shell and tube heat exchanger, which 
works on the Joule-Thomson effect on the shell side, before entering a 12.5 mm I.D. vacuum 
jacketed Pyrex visual test section, where the flow structure is captured using a high resolution 
CCD camera system. The nitrogen then travels through a fixed length of 304 stainless steel pipe 
(with I.D. 12.5 mm, thickness 1.27 mm, conductivity 14.7 W/m-K, and specific heat capacity 
J/kg-K at 30°C) and through the venturi flow meter before being dumped into an expansion tank, 
and then is evacuated from the facility via a ventilation system. This pipe is selected to ensure 
that any circumferential temperature variations are not “smeared out” due to conduction. 
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Figure 1. Cryogenic two-phase flow experimental facility. 


Validyne pressure transducers and type E thermocouples are used to measure the temperature 
and pressure at various locations. A highly sensitive differential pressure transducer measures 
the pressure drop across the visual test section, while inline thermocouples measure the 
temperature of the fluid entering and leaving. The flow then enters into the heat transfer 
measurement section of the facility, located 70 cm downstream of the visual test section. In this 
section two inline thermocouples measure the temperature of the fluid at the inlet and outlet, 
while a series of thermocouples are placed circumferentially on the surface of the pipe and on the 
outer surface of the melamine insulation, so that the temperature profile on the top, bottom and 
either side of the pipe is obtained (see Fig. 2). These thermocouples are kept in place using 
Teflon tape to ensure that adequate thermal contact is maintained. After the heat transfer section, 
the two-phase flow is passed through two heaters, which vaporize any remaining liquid, to 
ensure single-phase vapor enters the venturi. A thermocouple and pressure transducer located 
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prior to the venturi allow for the correction of compressibility effects to be taken into account in 
the mass flow rate computations. A Validyne differential pressure transducer measures the 
pressure drop across the venturi. The analog signals from the instruments are transferred and 
recorded by a 16-bit Measurement Computing data acquisition board and multiplexer system. 




Figure 2. Thermocouple arrangements along the pipe wall and along the insulation of the 

heat transfer section. 


The instruments have been independently calibrated and tested to ensure they function properly. 
The reported accuracy of the Validyne pressure transducers is ±0.25 % of full scale (0-200 psi 
for the static pressure transducers, 0-1.25 psi for the differential pressure transducer across the 
test section, and 0-12.5 psi for the differential pressure transducer across the venturi) and this 
was confirmed in the laboratory. The static pressure transducers are accurate to ± 3.5 kPa while 
the differential pressure transducer across the visual test section is accurate to within ± 0.06 kPa 
and that measuring the pressure drop across the venturi is accurate to within ± 0.35 kPa. Over 
the large temperature range experienced in these experiments the thermocouples used are 
accurate to within ±1.7° C. 

EXTRACTING THE HEAT TRANSFER COEFFICIENT 

The inverse heat conduction method employed in this study varies from previous methods in that 
it does not involve solving a system of least-squares equations. A comprehensive review of 
traditional inverse heat conduction methods is given in [7]. 

The process used here for extracting the transient heat transfer coefficient on the inside of the 
pipe involves a number of iterative steps, which are reported as three major steps. These major 
steps are illustrated in Fig. 3 and are carried out for each instance of time. 

This approach is based on the fact that it is possible to use numerical computations to solve the 
heat conduction equation in the wall of a pipe, and once the initial and boundary conditions are 
known the temperature profile may be calculated. Thus it is possible to guess the heat transfer 
coefficient inside the pipe and knowing the heat transfer coefficient on the outside of the pipe, 
the temperature profile is calculated. The temperature calculated at the outer wall is then 
compared with the measured outer wall temperature; if they match then the guessed heat transfer 
coefficient is taken as the actual heat transfer coefficient otherwise a new guess is made and the 
process is repeated until the computed and measured temperatures match. 
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COMPUTING THE TEMPERATURE FIELD IN THE PIPE WALL 

This procedure encompasses solving the heat conduction equation for the pipe wall. An integral 
part of this process is the implementation of the boundary conditions. The initial condition is 
simple as it is the state in which the pipe is in thermal equilibrium with the environment, a 
constant temperature throughout. 

The unsteady 3-D form of the heat conduction equation in cylindrical coordinates is written as 
follows: 
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where p is the density of the material, c is the heat capacity, t is the time variable, T is the 


temperature, r is the radial coordinate, z is the axial coordinate and ^ is the azimuthal 
coordinate. This equation is non-dimensionalized using the following parameters: 
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where T sat is the saturation temperature of the fluid within the pipe, T w is the temperature on the 
outer wall of the pipe, d is the pipe wall thickness, and the subscript 0 denotes the property is to 
be evaluated at the initial temperature, T 0 . Thus the original equation is thus reduced to (3) 
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It must be noted here that the thermal properties ( k and c ) of the pipe material vary significantly 
over the range of temperatures measured in the experiments. These variations are taken into 
account by curve- fitting the property values as functions of temperature and assigning the values 
based on the temperature at the given point. 

Equation (2) is discretized using a finite volume formulation. A backward Euler scheme is 
applied to the temporal term while a central difference scheme is used for the spatial terms. The 
resulting system of equations is solved using the Alternating Direction Implicit (ADI) method 
described in [8], 

Before the system of equations may be solved, the energy flowing into the system from the 
ambient must be known; hence the heat transfer coefficient for the outside of the pipe insulation 
must be determined. This is accomplished through a steady-state calibration process in which 
cool nitrogen vapor is passed through the heat transfer section at various mass flow rates. A 
relationship between heat flux into the pipe insulation versus temperature difference (ambient 
temperature less the insulation surface temperature) is determined. A constant heat transfer 
coefficient yields a linear relationship as shown in Fig. 4, and the slope is the heat transfer 

ft 

coefficient for the outer surface of the pipe. The measurement of q w is quite difficult due to the 
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small temperature rise in the cryogenic vapor, thus there is some scatter in Fig. 4. However, the 
measured heat transfer coefficient, 4.38 W/m 2 -K, is quite small, compared with convective heat 
transfer on the inner pipe wall. The errors in the calibration have a negligible impact on the 
computed thermal field in the pipe wall. 


Thus the outer pipe boundary condition may be written as, 


, dT 
k — 
dr 


h ou,( T ou< “O’ 


which is then non-dimensionalized to give, 
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Here h is the heat transfer coefficient, and subscripts out and oo represent the outer surface of 
the pipe, and the outer insulation surface, respectively. The insulation temperature is taken to be 
uniform at the measured temperature. 



Figure 3. Flow chart for transient heat transfer coefficient 
extraction process. 
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Figure 4. Calibration curve for determining heat transfer coefficient on the outer surface of the 

pipe - 

The heat transfer coefficient on the inside of the pipe is the quantity that is guessed during the 
calculation. However, because two-phase flow is present for much of the experiment, the flow 
structure is wavy-stratified and there is a significant difference between the temperature at the 
top of the pipe and the temperature at the bottom the pipe. This must be due to circumferential 
variations in the heat transfer coefficient. This problem is dealt with adequately by assuming 
that the heat transfer coefficient is constant in the region where only vapor is in contact with the 
pipe wall. In the region where the liquid film resides, a “geometry influencing factor” developed 
by Liao et al. [9] is used to circumferentially vary the heat transfer coefficient. The “geometry 
influencing factor” is function of the angle 6 shown in Fig. 5 and is determined to be inversely 
proportional to the heat transfer coefficient. 

The value of the geometry influencing factor, F{G) is given by, 
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Figure 5. Film boiling structure in stratified flow (Liao et al. [9]). 


Thus the inner pipe boundary condition may be written as, 
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Here the subscript in and fl represent the inner surface of the pipe, and the fluid flowing 
through the pipe, respectively. The fluid temperature is determined using an internal 
thermocouple that measures the temperature at the inlet of the heat transfer section. 


The heat transfer coefficient along the inside of the pipe varies with the thickness of the vapor 
film as shown in Fig. 6, and is also dependent on the height of the liquid film, which is 
determined using the images captured using the CCD camera. 


Region 1 : h in = h^ constant 



Figure 6. Distribution of heat transfer coefficient on the inside surface of the pipe. 
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Iteration Process for Guessing the Inner Heat Transfer Coefficient 


Any number of methods may be used to obtain a new guess for the inner heat transfer 
coefficient; a systematic method for iterating is recommended. In this work after the initial guess 
is made subsequent guesses are determined using linear interpolation or linear extrapolation. 
Hence the new guess for inner heat transfer coefficient is given by, 


h S' 


pres 


prev 


T 

\ pres 




- T 

known prev 


)+h 


prev 


prev J 


( 8 ) 


s is scaling factor that reduce oscillations and takes the value 0.1 in this work., and the 
subscripts pres , prev , new , and known denote the present quantity, the previous quantity, the 
new quantity, and the known quantity respectively. 

Test for Convergence 

The final step in the process is checking the computed temperature against the measured 
temperature. This is done using the temperature at the top and bottom of the pipe. Once the 
computed temperature is within ±1x1 O' 5 K, the temperatures are considered a match. Care must 
be taken when determining the limit within which to consider the temperatures to be matched; if 
the value selected is too small the computation time increases significantly and the marginal 
improvement in accuracy does not justify the increased computational cost. 

RESULTS AND DISCUSSION 

In order to access the performance and validity of this inverse method, it was first used to 
calculate the heat transfer coefficients for a single-phase flow simulation in which the heat 
transfer coefficient is known and varies in time and secondly it was used to calculate the actual 
heat transfer coefficient for single-phase nitrogen flowing through the experimental facility, with 
the results being compared with the predictions of the Dittus-Boelter correlation for cooling. In 
the first test case, the heat transfer coefficient follows a parabolic path with time. The results of 
these numerical tests are shown in Figures 7, and 8. The second case compares with actual 
experimental data. 

By examining Figure 7 it is seen that a parabolic varying heat transfer coefficient is be captured 
quite accurately with this approach, and from Figure 8 it is observed that the technique agrees 
very well with the predictions of the Dittus-Boelter correlation. 
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Figure 7. Computation of a parabolic varying heat transfer coefficient using the inverse method. 
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Figure 8. Comparison of heat transfer coefficient computed using inverse procedure and the 
Dittus-Boelter model for single-phase nitrogen flow. 
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Figure 9. Computed thermal properties with 
an average mass flux of 229 kg/m 2 -s; (a) 
temperature profile of top and bottom pipe 
wall temperature, (b) heat transfer 
coefficient at the top, (c) heat transfer 
coefficient at the bottom and (d) heat 
transfer coefficient at the bottom with the 
vertical axis re-scaled. 
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Figure 10. Computed thermal properties 
with an average mass flux of 77 kg/m -s; (a) 
temperature profile of top and bottom pipe 
wall temperature, (b) heat transfer 
coefficient at the top, (c) heat transfer 
coefficient at the bottom and (d) heat 
transfer coefficient at the bottom with the 
vertical axis re-scaled. 


Time (sec) 



Time (sec) 

Figure 1 1 . Mass flux variation for experiment with average mass flux of 229 kg/m -s. 


I 



Time (sec) 

Figure 12. Mass flux variation for experiment with average mass flux of 77 kg/m 2 -s. 

Since the technique performed satisfactorily for the single-phase case it was then applied to the 
chilldown experimental data. The results shown here are for two experimental cases, with 
different average mass fluxes; one at 229 kg/m 2 -s and the other at 77 kg/m 2 -s (the actual mass 
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flux values are given in figures 1 1 and 12 respectively). The results are shown in Figs. 9 and 10. 
A comparison of Figs. 9 and 10 leads to the following observations: 

• The chilldown time is much shorter for the higher mass flux case as seen by comparing 
9a and 10a. This is to be expected since the higher mass flux case yields higher heat 
transfer coefficients. 

• The heat transfer coefficient at the bottom of the pipe is at least one order of magnitude 
greater than the heat transfer coefficient at the top, which is seen by comparing Figs. 9b 
and 9d and Figs. 10b and lOd. This is expected, as the lower region of the pipe is where 
the liquid phase resides. 

• The transition between film boiling and nucleate boiling (Leidenfrost point) is clearly 
seen as reflected in the dramatic increase in heat transfer coefficient at the bottom of the 
pipe, as seen in Figs. 9c and 10c. 

An interesting point to note is that the maximum heat transfer coefficient reached in the lower 
mass flux case is actually greater than that reached in the higher mass flux case, see Figs. 9c and 
10c. One possible reason for this is that in the lower mass flux case there may be a greater 
amount of thermal energy remaining in the pipe wall, which drives ebullition from the surface. 

In the higher mass flux regime, ebullition may be suppressed. 

It is also observed that the unsteady nature of the flow is reflected in the heat transfer 
coefficients, which is best seen by comparing the heat transfer (Figs. 9b and 9c) to the 
instantaneous mass flux values (Fig. 1 1) for the experiment with average mass flux of 229 
kg/m 2 -s. 


The initial phases of the two experiments exhibit different behavior, which is seen in Figs. 1 1 
and 12. In the lower mass flux case, a surge of vapor was observed when the nitrogen was 
introduced into the system, which was not observed in the higher mass flux case. These 
behaviors were captured in the heat transfer computations. 

CONCLUSION 

In this study an inverse numerical procedure for predicting the transient heat transfer coefficient 
for two-phase flows has been developed. It is seen to perform well for the single-phase test 
cases. The large circumferential variation of the heat transfer coefficient that occurs during the 
cryogenic chilldown process is clearly highlighted and the Leidenfrost point is clearly identified. 
This technique will facilitate the development of heat transfer coefficient correlations for 
cryogenic chilldown. 


NOMENCLATURE 
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heat flux 
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temperature M 
time [ 5 ] 

temperature difference [A'] 
cylindrical coordinates 
non-dimensional temperature 
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Subscripts 

0 

reference value 

bottom 

bottom of the pipe 

fl 

fluid flowing inside pipe 

in 

inner surface of the pipe 

known 

measured quantity 

new 

next computational step 

out 

outer surface of pipe 

pres 

current computational step 

prev 

previous computational step 

sat 

saturated condition 

top 

top of the pipe 

w 

pipe wall 

00 

outer insulation surface 

Superscript 

c 

dimensionless variable 
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PROJECT GOALS 

The efficient and safe utilization of cryogenic fluids in thermal management, power and 
propulsion, and life support systems during space missions involves the transport, handling, and 
storage of these fluids in reduced gravity. Due to their low boiling points, two-phase flows are 
encountered in most of the cryogenic operations. The complexity of the problem results from the 
intricate interaction of the fluid dynamics and heat transfer, especially when phase-change 
(boiling and condensation) is involved. Because of the large stratification in densities between 
the liquid and vapor phases, the reduced gravity condition in space would strongly change the 
terrestrial flow patterns and accordingly affect the momentum and energy transport 
characteristics. This research focuses on addressing specific fundamental and engineering issues 
related to the microgravity two-phase flow and heat transfer of cryogenic fluids that require 
advanced numerical simulations in concert with validating experimentation. The outcome of the 
research will provide cryogenic practitioners the practical analysis tools required to address 
specific issues impacting the safe and efficient operation of liquid hydrogen cryogenic transport 
systems in reduced gravity. 

ACCOMPLISHMENTS 

ABSTRACT 

For the cryogenic two-phase flow and heat transfer study, we have focused on the liquid 
nitrogen convective film boiling process in both terrestrial gravity and microgravity. Two 
types of experiment were performed. One is aimed at the chilldown phenomena of an 
unheated tube. The other investigates the constant wall heat flux heating condition. An 
experimental apparatus was designed and built. Flow pattern and wall temperature history 
during the chilldown and constant heating condition were obtained. We have measured 
the wall temperature history at three downstream locations. The terrestrial chilldown 
process is divided into three stages based on the temperature of the tube wall. At the 
beginning of the chilldown process, the wall temperature of the test section is near the 
room temperature, which is much high than that of the liquid nitrogen. The initial stage is 
therefore associated with the film boiling that causes the two-phase flow to take the form 
of a vapor core with scattered small liquid fragments. When the wall temperature 
decreases further, the two-phase flow is experiencing a transition from the film boiling to 
the stratified flow. Finally when the wall is chilled down enough, a stable stratified two 
phase flow was observed. Flow visualization under the microgravity condition was 
performed and it was found that suspended liquid chunks were floating in a vapor stream. 

For the constant heating case, two different flow rates were investigated. Both wall 
temperatures and flow visualizations were attempted. The wall temperature and flow 
pattern would achieve relatively steady state condition under a constant heating 
condition. For the lower flow rate case, the two-phase flow pattern for the terrestrial case 
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is influenced by the wall temperature level and the gravity in very similar ways as the 
chilldown case. For the higher flow rate case, the gravity effect is negligible. Therefore, 
the results are similar between the terrestrial and the microgravity cases. 

As a first step towards the film boiling numerical study, a direct numerical simulation of 
bubble dynamics with large liquid-to-vapor density ratio and phase change has been 
developed using a sharp-interface moving boundary numerical method. The numerical 
techniques are based on a fixed-grid, finite volume method capable of treating the 
interface as a sharp discontinuity. The unsteady, axisymmetric Navier-Stokes equations 
and energy equation in both liquid and vapor phases were computed. The mass, 
momentum and energy conditions are explicitly matched at the phase boundary to 
determine the interface shape and movement. The cubic B-spline is used in conjunction 
with a fairing algorithm to yield smooth and accurate information of curvatures. For a 
stationary and growing bubble, the predicted growth rate was found to agree with the 
theoretical limit of t 1/2 . However, for a rising and growing bubble, the predicted growth 
rate was approaching the theoretical limit of t 2/3 for a spherical bubble but did not reach it 
exactly owing to bubble deformation. For a rising and growing bubble, the location of the 
highest heat flux was found and the deviation largely depends on the Reynolds number. 

I. INTRODUCTION 

Role of Cryogenics in Space Exploration 

The extension of human space exploration from a low earth orbit to a high earth orbit, 
then to Moon, Mars, and possibly asteroids and moons of other planets is one of NASA’s 
biggest challenges for this new millennium. Integral to this is the effective, affordable, 
and reliable supply of cryogenic fluids. The efficient and safe utilization of cryogenic 
fluids in thermal management, power and propulsion, and life support systems of a 
spacecraft during space missions involves the transport, handling, and storage of these 
fluids in reduced gravity. The uncertainties about the flow pattern and heat transfer 
characteristics pose a severe design concern. Moreover, the thermo-fluid dynamics of 
two-phase systems in reduced gravity encompasses a wide range of complex phenomena 
that are not understood sufficiently for engineering design to proceed. For the two-phase 
flow and heat transfer study, we have focused on the forced-convection film boiling 
process owning to the low boiling points for the cryogenic fluids. The objective is 
therefore to provide heat transfer coefficients and two-phase flow patterns that are useful 
for the initial phase (film boiling) of the chilldown process of cryogenic fluid flow and 
transport. 

Chilldown process, which is important during initial cryogenic fluid transport, is 
however not fully understood. In this paper, liquid nitrogen chilldown process was 
investigated and divided into three stages based on the flow image and temperature data. 

Due to their low boiling points, boiling and two-phase flows are encountered in most of 
the cryogenic operations. The complexity of the problem results from the intricate 
interaction of the fluid dynamics and heat transfer, especially when phase-change (boiling 
and condensation) is involved. Because of the large stratification in densities between the 
liquid and vapor phases, the reduced gravity condition in space would strongly change 
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the terrestrial flow patterns and accordingly affect the momentum and energy transport 
characteristics. Boiling and two-phase flow behave quite differently when the gravity 
levels are varied. The proposed research will focus on addressing specific fundamental 
and engineering issues related to the microgravity two-phase flow and boiling heat 
transfer of cryogenic fluids that require well-designed and meaningful experimentation. 
The outcome of the research will provide fundamental understanding on the transport 
physics of cryogenic boiling and two-phase flows in reduced gravity. 

Boiling Curve 

A boiling curve shows the relationship between the heat flux that the heater supplies to 
the boiling fluid and the heater surface temperature. According to the typical boiling 
curve as shown in Fig. 1, a chilldown (quenching) process usually starts from point E and 
then goes towards point D in the film boiling regime as the wall temperature decreases. 
Point D is called the Leidenffost point which signifies the m inimum heater temperature 
required for the film boiling. For the film boiling process, the wall is so hot that liquid 
will vaporize before reaching the heater surface which causes the heater to be always in 
contact with vapor. When cooling beyond the Leidenffost point, if a constant heat flux 
heater were used, then the boiling would shift from film to nucleate boiling (some where 
between points A and B) directly with a substantial decrease in the wall temperature 
because the transition boiling is an unstable process. 


Free convection 


Boiling regimes 
Nucleate Transition 


Isolated Jets and 
bubbles columns 



Figure I 


Typical boiling curve 
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The Effects of Gravity 

Because of many difficulties and limitations involved in the drop tower microgravity 
experiments, meaingful terrestrial experiments would therefore serve as a valuable tool to 
provide quidance and check for microgravity experimentation. The following is a scaling 
analysis that examines the gravity effects. For annular-flow film boiling, the effect of 
gravity is assessed based on the ratio of Gr/Re 3 , where Gr is the Grashof number and Re 
is the Reynolds number. The gravity effect is measured by the natural convection 
contribution which is characterized by the Gr while the forced convective film boiling is 
scaled by the Reynolds number. According to Gebhart et al. (Buoyancy-Induced Flow 
and Heat Transport, Hemisphere Publishing Company, 1988), Re 3 is used in the 
denominator when the flow is perpendicular to gravity for a horizontal tube. Re 2 is used 
when the flow is in the same direction of gravity. All the thermal properties are those of 
vapor because of film boiling. Fig. 9 shows a plot for this ratio with the flow velocity 
range of 0-1 m/s. 

Based on Figure 2, we would also perform terrestrial gravity experiment with the vapor 
velocity greater than 10 cm/s which corresponding to Gr/Re 3 less than 0.2. The natural 
convection is negligible for Gr/Re 3 less than 0.225 according to Gebhart ct al. (1988). 
For a vapor velocity of 10 cm/s, the corresponding liquid velocity is approximately 0.02 
cm/s which is the velocity of liquid when it first exits the tank. 



Vapor Velocity [m/s] 

Figure 2. Gr/Re 3 vs. vapor velocity 
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II. LITERATURE REVIEW 

Momentum and energy transport in two-phase flow and boiling under microgravity is a 
relatively new field. This is a complex problem for the scientific community to solve. The 
rewards are rich, however, as new ways to propel future space missions and more 
intensive space activities have been planned. Currently, there exists a limited quantity of 
literature, data, and test facilities to even begin to answer these questions. Microgravity 
cryogenic boiling research was initiated more than fifty years ago, but the progress has 
been limited. The following summarizes the accomplishments. 

II.l Terrestrial Cryogenic Boiling and Two-Phase Flow 

Numerous studies of cryogenic boiling in the one-g environment were conducted in the 
1950’s and 60’s. Reference [1] is a comprehensive review of the experimental studies and 
heat transfer correlations. For the fluids oxygen, nitrogen, hydrogen and helium, it was 
found that for pool boiling, the Kutateladze [3] correlation had the greatest reliability for 
nucleate boiling while the Breen and Westwater [4] correlation was best for film boiling. 
Maximum nucleate flux data were reasonably well predicted by the Kutateladze 
correlation. Although these correlations were selected as the best available, neither has 
particularly good agreement with experimental data. For the case of forced convection 
boiling, Reference [1] reported that no correlation was found to be distinctly better. Some 
simple predictive methods were found to work as well as more complex schemes. In all 
boiling cases, it was questioned as to whether or not the predictive correlations include all 
of the significant variables that influence the boiling process. In particular, it was 
suggested that more detailed and better controlled experiments are needed and that more 
attention to surface and geometry effects is required. 

Reference [2] is another comprehensive review of cryogenic boiling heat transfer 
addressing hydrogen, nitrogen and oxygen. It was reported that nucleate pool boiling 
results cannot be correlated by a single line but cover a range of temperature difference 
for a given heat flux. The spread is attributed to surface condition and geometry, and 
orientation. Maximum heat flux can be reduced by about 50% when going from normal-g 
to near zero-g. Reference [2] reported a fair amount of data for film pool boiling. Film 
boiling heat flux is reduced considerably at near zero-g conditions. Only a very limited 
amount of data is available for subcooled or saturated forced convection boiling and few 
conclusions were drawn. The lack of data for cryogenic forced convection boiling is also 
reported by Brentari and Smith [5]. 

Relatively recent correlations have been published for normal-g saturated flow boiling of 
cryogens [6-9] using the Convection, Boiling and Froude numbers as correlating 
parameters. 
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II. 2 Reduced Gravity Cryogenic Boiling and Two-Phase Flow 

In general, there is little heat transfer data for cryogenic flow boiling in reduced gravity. 
We were able to find just two reports. Adham-Khodaparast et al. [10] have investigated 
the flow film boiling during quenching of R-113 on a hot flat surface. They used 
microsensors to record instantaneous heat flux and heater surface temperature during the 
film boiling on board a KC-135 aircraft. They reported lower heat transfer rates during 
microgravity as compared to normal gravity and contributed that to thickening of the 
vapor layer. The wall superheat and the surface heat flux at the onset of rewetting and 
maximum heat flux were found to increasde with the inlet liquid subcooling, mass flux 
and gravity level. The effect of gravity was determined to be more important for low flow 
rates and less relevant for high flow rates. Antar and Collins [11] reported flow 
visualization and measurements for flow film boiling of liquid nitrogen in tubes on board 
KC-135 aircraft. They were particularly interested in the vapor/liquid flow pattern and 
the thermal charateristics. They identified a new vapor/liquid flow pattern that is unique 
in low gravity. Antar and Collins [11] observed that a sputtering leading core followed by 
a liquid filament annular flow pattern. This new flow pattern is composed of a long and 
connected liquid column that is flowing in the centerof the tube and is surround by a thick 
vapor layer. The vapor annulus that separates the liquid filament from the wall is much 
thicker than that observed in the terrestrial experiment. They attributed the filamentary 
flow to the lack of difference in the speed of vapor and liquid phases. On the heat 
transfer side, they reported that the quench process is delayed in low gravity and the tube 
wall cooling rate was diminished under microgravoty conditions. 

III. EXPERIMENTAL SYSTEM 
Flow Loop 

Figure 3 shows the schematic of the boiling and two-phase flow apparatus and Figure 4 is a 
photographic view of the this apparatus located in two side-by-side cubicles for drop tower 
microgravity experiment. The experimental loop mainly consists of a bellows-controlled 
reservoir, a cryogenic valve, test section inlet portion, test section, test section outlet portion, 
liquid separator, lighting and video system. A vacuum jacket is provided to insulate the glass 
tube test section. 
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Figure 3. Schematic of cryogenic boiling and two-phase flow test apparatus. 


Figure 4. Experimental apparatus located in cubicles 
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Reservoir 

Liquid nitrogen is supplied by a cylindrical stainless steel reservoir that is fitted with a 
bellows driven by a constant speed motor to provide a constant flow rate. Form 
insulation is applied to the reservoir and the whole loop besides the test section. 

Test Section 

Two different test sections were used in the experiments. The first is an unheated Pyrex 
glass tube which was used for the chilldown study while the second test section is a 
heated quartz tube coated with a thin film of Tin Oxide to provide a constant heat flux. 
The non-heated test section is a Pyrex glass tube of 25-cm long. The inner and outer 
diameters of the test section are 11.1 mm and 15.8 mm respectively. At both ends of the 
test section, stainless steel flanges and adaptors are used to connect the test section to the 
test section inlet and outlet portions. Also some traditional methods were adopted to seal 
the glass to the steel. A total of 16 type-T thermocouples are placed on the glass tube, 9 
are embedded very close to the inner surface of the tube wall through drilled holes. These 
9 thermocouples are placed at three downstream locations. At each downstream location, 
there are three thermocouples located circumferentially at equal separation distance. The 
other 7 thermocouples are used to measure the outside wall temperatures of the glass 
tube. Fig. 5 illustrates the locations of the thermocouples. Fig. 6 is a photograph of the 
test tube with attached thermocouples. 



Figure 5. Locations of the 16 thermocouples on test section 
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Figure 6. Photograph of cryogenic two-phase flow test apparatus 


The heated quartz tube has an OD of 9.5mm (3/8”) and ID of 7.0mm. The length of the 
quartz tube is 16.3cm (6.4”). A transparent thin film of Tin Oxide covers 8.6cm of the 
tube. The thin film provides a constant heat flux when heated by a DC power supply. At 
each ends of the thin film an electrode with a length of 1 .3 cm (0.5”) connects the quartz 
tube to the power supply. 

Because its wall is so thin as compared to the previous Pyrex tube, drilling holes on the 
quartz tube for inserting thermocouples is not feasible. To obtain the circumferential 
temperature measurements in this situation, we put three teflon ferrules on three locations 
of the quartz tube. On each of the ferrules, tiny holes are drilled to let the thermocouples 
go through such that they can be fastened. Figure 7 shows the schematic of the heated 
quartz tube and provides the locations of the thermocouples. 



Figure 7. Sketch of the heated quartz tube 
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Flow Section 

Standard 3/8” stainless steel tubes are used for the flow loop besides the test section. 
Swagelok fittings with teflon ferrules were used at the glass to metal connections. 

Cryogenic Valve 

A cryogenic valve with !4” NPT threads on both sides is used to connect the reservoir to 
the test section inlet. The cryogenic valve is electronically controlled. During the 
experiment, after the reservoir is cooled the valve is opened and liquid nitrogen is 
released to chill down the test section. 

Video System 

The video system consists of a high speed CCD camera (CV-730) with a V2” HAD 
Interline Transfer CCD sensor and fluorescent light. The camera’s electronic shutter 
speed was from 1/1,000 to 1/10,000 sec. A Ceptre PD-1212ARD5 (12 VDC 1.2 A) power 
supply was used to provide power to the camera. Fluorescent light can provide 
comparably enough light and generate less heat. 

Data Acquisition Component 

Signals from the thermocouples are sent to the computer through a temperature 
acquisition board plugged into the PCI slot. Video images are also monitored and 
recorded by connected the CCD camera to the computer with a frame grabber board. 
During microgravity tests, long thermocouple extension wires are used to connect the 
thermocouple to the temperature acquisition board, so that temperature data can be 
collected during the drop with comparable accuracy. When the heated quartz tube is used, 
the constant heat flux was maintained by a HP E3 631 A DC power supply. 

For the microgravity test, all the experimental components are put into and fastened to 
two standard rigs. A portable vacuum pump is also attached onboard to maintain the 
vacuum condition for the test section during the microgravity test. Figure 5 shows the 
whole system inside the standard rigs. 

Drop Tower Apparatus 

The microgravity environment was provided by the University of Florida drop tower. The 
drop tower is 5-story high and is equipped with an airbag for deceleration. Figure 8 
shows the overall drop tower schematic. At the current drop height of 15.25 m, UF’s drop 
tower produces 1.7 seconds of free fall. The experiment is actually located in a drag 
shield which prevents the experiment from the aerodynamic drag. As a result, the 
microgravity level is measured between 10‘ 5 to 10' 4 g. 
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Figure 8. A schematic of the drop tower system 

IV. EXPERIMENTAL RESULTS 

Two types of experiments were conducted, the first type is the so called “pipe chilldown” 
experiment, where there is no heater associated with the test tube. During the chilldown, 
pipes experience a fast cooling on the walls. The purpose of the chilldown experiment is 
to provide an overall calibration of the system that includes both flow visualization and 
heat transfer measurements, and also to offer flow patters and heat transfer characteristics 
during pipe chilldown process that takes place in many cryogenic applications. The 
second type experiment is aimed at a fundamental understanding for a steady state 
condition, where the tube wall outer surface is coated with a thin Tin oxide-film to 
maintain the wall at a constant heat flux. This type of experiment is intended for the 
understanding of the two-phase flow and heat transfer under a steady heating boundary 
condition. 

IV. 1. Chilldown Experiment 

During the chilldown experiment, liquid nitrogen enters the Pyrex glass tube that is at the 
room temperature (25 °C). The bulk velocity of the two-phase flow is estimated at 5 cm/s. 
From Fig. 2, it is clear that the two-phase flow investigated here is totally dominated by 
gravity. Because the glass tube is unheated, the tube wall is going through a transient 
cooling process. In the following, the measured wall temperature history is presented first 
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and then flow visualization pictures are given to provide the correlations between the 
wall temperature and flow pattern. 

Tube Wall Temperature History 

Terrestrial Gravity Case 

For the chilldown process, the temperature histories of the tube walls were measured by the 
thermocouples at three downstream locations as shown in Fig. 4. The results are given in Fig. 9. 
Based on the characteristics of the temperature curves, the chilldown process is generally 
divided into three stages. At the beginning of the process, the wall temperature of the test section 
is near room temperature, which is much higher than the temperature of the liquid nitrogen. 
Based on the flow visualization presented later, the initial stage is associated with a hot wall 
which causes the two-phase flow to take the form of a vapor core with scattered small liquid 
chunks. In general, the temperatures measured by thermocouples that are located in the bottom 
of the tube (Ch 12, 15, 9, 6 and 3) are lower than those from thermocouples embedded in the 
upper portion. This is mainly due to the gravity effects that bring direct liquid impingement 
cooling to the bottom walls. 

A distinctive character of the bottom wall temperature history is the sharp increase in the slope 
at the wall temperature of about 190 K. This phenomenon is seen at all three downstream 
locations. The physical explanation is based on the transition of boiling regimes from film 
boiling to nucleate boiling as mentioned earlier and shown in Fig. 1. 

Finally when the wall is chilled down enough such that the wall temperature (< 140 K, point A 
in Fig. 1) is too low to support nucleate boiling, the heat transfer switches to single-phase forced 
convective boiling. 



(a) Temperature profile of section 1. 
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(b) Temperature profile of section 2. 



(c) Temperature profile of section 

Figure 9. Temperature profiles at three downstream locations during the chilldown process. 
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Based on the temperature data, the heat transfer coefficient history for the bottom wall at 
Section 1 is estimated and plotted in Figure 10. The magnitudes of the heat transfer 
coefficient really reflect the boiling regimes, which conforms with the fact that nucleate 
boiling is the most efficient and the film boiling is much less effective. The horizontal 
line indicates the heat transfer coefficient value calculated using the steady state nucleate 
boiling correlation developed by Van Dresar et al. [9], The correlation value sits in the 
middle of our experimental data points which confirms that nucleate boiling is indeed 
responsible for the highest sloped portion of the temperature curve. 



Figure 10. Estimated heat transfer coeficients 


Microgravitv Case 

Because that our drop tower can only provide 1 .7 seconds of microgravity period, it is not 
feasible for the chilldown experiment. However, we can speculate that the wall 
temperature history curve in microgravity would be very similar to those of the two upper 
portion locations (Ch. 11 and 14 or 13 and 16). In microgravity, there would be no 
stratification and perfect inverted annular flow would prevail, which causes the vapor to 
cover all the tube wall and leave liquid in the central core area. 

Flow visualization 

Terrestrial Condition 

The experimental apparatus was first used to provide flow visualization results for the 
chilldown process of a Pyrex glass tube. This initial phase of experimentation would 
provide not only some physical understanding of the two-phase flow characteristics 
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during the chilldown process, but also the performance evaluation for the apparatus 
which includes both image and data acquisition. The chilldown process is divided into 
three stages based on the temperature of the tube wall as discussed above. At the 
beginning of the chill down procession, the wall temperature of the test section is near the 
room temperature, which is much higher than the temperature of the liquid nitrogen. 


The initial stage as shown in Fig. 1 1 is therefore associated with a hot wall which causes 
the two-phase flow to take the form of a vapor core with scattered small liquid chunks. 



Figure 11. Initial stage of chilldown - vapor core Figure 12. Initial stage - inverted 
with scattered liquid chunks. Annular Flow. 



Figure 13. Transition stage Figure 14. Transition stage 


When the wall becomes slightly cooler, the inverted annular flow is seen as given in Fig. 
12 where the liquid is occupying most of the tube with a vapor layer adjacent to the tube 
wall. For both cases as shown in Figures 11 and 12, it is believed that film boiling is the 
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main mechanism (portion between D and E in Fig. 1). When the wall decreases in 
temperature further, the two-phase flows are experiencing a transition from the film 
boiling (inverted annular flow) to the stratified flow. Figures 13, 14 and 15 show the 
transition flows for warm tube walls where liquid chunks are starting to accumulate near 
the bottom side of the tube wall due to the gravity effects. Finally when the wall is chilled 
down enough, a stable stratified two phase flow was observed. As shown in Figures 16, 

17, and 18, a stable and continuous liquid layer is located on the bottom of the tube and a 

pure vapor core is flowing on top of the liquid layer. It is believed that for the stratified 

flow, nucleate boiling (portion between A and B in Figure 1) is taking place between the 

liquid layer and the cold bottom wall while the film boiling is the mode between the » 

vapor core and the tube wall. Therefore under the stratified flow, there is a J; 

circumferential temperature variation on the tube wall. 


Based on our estimation from the video 



Figure 15. Transition stage 


images, the velocity of the vapor is 
around 5 cm/s. Therefore the gravity 
effect overpowers the flow inertia 
according to Figure 2. This is the reason 
why we observed the stratified flow in 
this terrestrial experiment. 



Figure 16. Cold wall - stratified flow. Figure 17. Cold wall - stratified 

flow. 
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Microgravity Condition 

Figure 19 provides a series of four pictures for the chilldown experiment in a 
microgravity condition. We followed a large chunk of liquid as it traveled downstream. It 
is seen that basically it is an inverted annular flow where the vapor is filling the tube with 
liquid chunks in the center. The size of the liquid chunk is decreasing due to vaporization 
as it flows downstream. 




(C) (d) 


Figure 19. Photographs of flow pattern in microgravity 

IV.2. Constant Wall Heat Flux Experiment 

Since the wall is heated at a constant heat flux, the flow and heat transfer inside the tube is 
expected to develop into a steady condition. The applied heat flux is 13002 w/m 2 . Two different 
motor speeds of 5 and 15 rpm were used in the constant heat flux experiment. The bulk velocities 
are estimated at 7.3 and 23.5 cm/s for 5 and 15 rpm, respectively. As a result, the 5 rpm case is 
somewhat affected by the gravity and the 15 rpm case is almost independent of gravity as 
indicated by Fig. 2. For the 5-rpm case, we only have terrestrial results while for the 15-rpm 
case, both terrestrial and microgravity data were obtained in a single run. 

5-rpm TerrestrialEexperiment 

As shown in Figure 7 for the 5 rpm case, wall temperatures were measured at three 
circumferential points for each of the three downstream locations. Because of similarities, only 
the temperatures at first downstream location are plotted and given in Figure 20. Curve 1 shows 
the temperature recorded at the 12 o’clock point and curves 2 and 3 are for the 4 and 8 o’clock 
points, respectively. After 275 seconds, the system has reached relatively steady state. As 
mentioned above, the 5-rpm case is affected by gravity to some extent and the flow visualization 
discussed next supports this prediction. The gravity effects would lead to circumferential 
temperature differences. The temperatures at the two lower points are about 20 K colder, that is 
due to more liquid fragments in the lower half due to gravity. Figure 21 gives the corresponding 
heat transfer coefficients for the same three points. 

Figure 22 provides the flow visualization photographs for the 5-rpm terrestrial experiment. As 
compared to the chilldown experiment where the wall is unheated, the general size of the liquid 
phase is smaller and liquid fragments are more widely dispersed for the heated wall case. The 
effects of gravity are obvious that more liquid pieces in the lower half of the tube. 
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Figure 21. Heat transfer coefficient at three circumferential locations. 
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Figure 22. Flow visualization of heated tube under terrestrial condition. 

15-rpm Terrestrial and Microgravitv experiment 

For the 15-rpm experiment with a heated tube, identical heating was used as the 5-rpm terrestrial 
case. But the flow rate is about three times higher. It is noted that for the combined experiment, 
terrestrial results precede the microgravity counterparts. First, the wall temperature measurement 
results are given in Figure 23. Because of the limitations in a drop tower experiment, 
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temperatures were measured only in two locations, one near the inlet and the other close to the 
outlet. Also we only measured one point (the 12 o’clock position) at each downstream location 
due to the expected circumferential similarity in microgravity. According to the discussion 
above, the gravity effect is supposed to be minimum and indeed the wall temperature 
measurement and the flow visualization both reflect this theoretical prediction. 



Time (s) 

Figure 23. Wall temperature measurement results 

In Fig. 23, the temperature levels registered by the thermocouple near the inlet exhibit only a 
very slight decrease during the microgravity period while near the outlet, a slight increase is 
measured. The temperature difference between the two downstream locations is a general feature 
of constant wall heating. To explain the microgravity results, we need to take a look at the flow 
visualization pictures as shown in Fig. 24. Because of the higher speed in the 15-rpm case, the 
liquid is broken down more in the form of mist. For the terrestrial condition, the mist seen as 
gray-shaded portion blankets the bottom wall as a thin layer. While in microgravity, the flow is a 
perfect inverted annular pattern where the liquid mist is occupying the central core. There are 
two main effects due the absence of gravity. The first is a slight acceleration of the flow as a 
result of the decrease in wall friction due to the replacement of the mist layer by the vapor 
adjacent to the wall. The second is the loss of the circumferential conduction in the wall. In 
terrestrial gravity, there is a circumferential temperature gradient in the wall that results in the 
conduction flow of heat from the top half to the bottom half. The slight decrease in the wall 
temperature near the inlet is thought to be due to the first effect as a faster flow gives rise to a 
higher heat transfer coefficient, while the slight increase near the outlet is due to the second 
effect that stops the flow of heat from the top of the wall. 
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Before Drop: 1-g During Drop: micro-g 

Figure 24. Flow visualization of constant wall heat flux tube 

V. NUMERICAL SIMULATION 

As a first step towards the numerical simulation of film boiling, a detailed treatment of the 
liquid-vapor interface with phase change should be performed. We have chosen a single 
bubble during a dynamic translational motion as the model case. In the numerical 
simulation, a single bubble is rising through a superheated liquid, so the bubble is growing 
due to boiling and experiencing a shape change due to hydrodynamic force. A fixed-grid 
direct numerical simulation method has been developed and verified for single bubble 
boiling and condensation phenomena. The original contribution of this research is the first 
introduction of a sharp interface between the vapor and liquid phases, which allows the 
simulation of density ratios up to the order of 1000. The interface location is obtained as 
part of the solution which accounts for all coupled dynamic and thermal effects. 

For a stationary and growing bubble, the predicted growth rate was found to agree with the 
theoretical limit of t 1/2 where t is the time. However, for a rising and growing bubble, the 
predicted growth rate was approaching the theoretical limit of t 2/3 for a spherical bubble but 
did not reach it exactly owing to bubble deformation. 

V.l Numerical Method 

In the present approach, we employ a combined Eulerian-Lagrangian strategy. The fixed 
Cartesian grid is used as the Eulerian framework of the algorithm to facilitate the field 
equation computation. Within this framework, separate marker points, connected by piece- 
wise polynomials, are adopted to represent the interface to form the Lagrangian portion of 
the method. With moving boundaries, the motion of the interface is tracked through the 
translation of the marker points over the stationary, Cartesian grid. In each phase, a finite 
volume, fractional step method is employed to solve the coupled governing equations of 
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momentum, energy and mass transfer. A cut-cell approach is developed to handle arbitrary 
intersections between an interface and the grid line. In the interface region, the grid will be 
recombined to form non-Cartesian cells. Consistent interpolation formulas are chosen for 
estimation of the fluxes along any of the cell surfaces. Both inviscid and viscous terms can 
be handled to maintain a globally second-order accurate algorithm. Here, this approach is 
further extended to treat liquid-vapor interaction. The schematic of the computational 
domain is given in Fig. 25. 

V.2 Results 

Fig. 26 shows the bubble shape changes for varying density ratios without phase change. 
The system condition of the bubble are Re (Reynolds number) = 100, We (Weber number) 
= 4, Fr (Froud number) = 1.0. Fig. 27 represents the effects of Reynolds number and the 
Weber number. The density ratio is 1605 and the dynamic viscosity ratio is 22. Fig. 28 is a 
summary of the bubble shape change for various Reynolds and Weber numbers. Fig. 29 
shows the bubble shapes for a rising bubble with phase change with density effects. The 
system condition is Re = 100, We = 4, Fr = 1.0. Fig. 30 shows the streamline structures for 
the same cases in Fig. 29. 



Figure 25. Schematic of computational system 
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Figure 27. Buoyancy Driven Bubble without Phase Change: 
Reynolds number and Weber number Effect 
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Figure 28. Summary of bubble shapes for different Re and We. 
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Figure 30. Streamline structures of the bubbles. 


V.3 Conclusion of Numerical Stud> 


■ A direct numerical simulation capability has been developed for bubble dynamics 
with phase change. The flow, thermal fields and interface motion are all obtained as 
part of the solution process. 

■ The methodology is unique in terms of the combined capabilities of resolving sharp 
interfaces, honoring mass conservation, handling large property ratios and phase 
change. 

■ Major dimensionless parameters, with and without phase change, are varied to 
probe the physical and numerical issues involved. 

■ For a stationary and growing bubble, the predicted growth rate was found to agree 
with the theoretical limit of t 1/2 where t is the time. However, for a rising and 
growing bubble, the predicted growth rate was approaching the theoretical limit of 
t 2/3 for a spherical bubble but did not reach it exactly owing to bubble deformation. 
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PROJECT TITLE: CHILL DOWN PROCESS OF HYDROGEN 
TRANSPORT PIPELINES 

Task PI: Dr. Renwei Mei, Department of Mechanical & Aerospace Engineering 
Co-I: Dr. James Klausner, Department of Mechanical & Aerospace Engineering 

PROJECT GOALS 

The objective of this research task is to experimentally and computationally study the unsteady 
dynamics of a liquid hydrogen wave front translating down a pipeline during the chill down 
mode and develop a comprehensive computational model to predict the associated flow fields, 
thermal fields, and residence time. In order to accomplish the objective, the following tasks are 
planned: a) experimentally study the dynamics of a propagating liquid hydrogen wave front, b) 
experimentally measure the heat transfer rate associated with film boiling beneath the wave 
front, c) experimentally investigate transient nucleate boiling heat transfer associated with liquid 
film flow, d) develop a comprehensive computational model to predict the flow and temperature 
fields associated with propagating liquid hydrogen waves, and e) develop engineering models for 
cryogenic practitioners to predict the required chill down residence time under a variety of 
different operating conditions. 

ABSTRACT 

A pseudo-steady model has been developed to predict the chill-down history of pipe wall temperature in 
the horizontal transport pipeline for cryogenic fluids. A new film boiling heat transfer model is 
developed by incorporating the stratified flow structure for cryogenic chill-down. A modified nucleate 
boiling heat transfer correlation for cryogenic chill-down process inside a horizontal pipe is proposed. 
The efficacy of the correlations is assessed by comparing the model predictions with measured values of 
wall temperature in several azimuthal positions in a well controlled experiment by Chung et al. (2004). 
The computed pipe wall temperature histories match well with the measured results. The present model 
captures important features of thermal interaction between the pipe wall and the cryogenic fluid, 
provides a simple and robust platform for predicting pipe wall chill-down history in long horizontal pipe 
at relatively low computational cost, and builds a foundation to incorporate the two-phase hydrodynamic 
interaction in the chill-down process. 

INTRODUCTION 

The cryogenic chill-down is encountered in many applications but is of particular importance in 
cryogenic transportation pipelines. For example, in rockets or space shuttle launch facility, cryogenic 
liquids as fuel are filled from the storage tank to the internal fuel tanks of a space vehicle through a 
complex pipeline system. To avoid evaporated fuel entering the space vehicle, a cryogenic chill-down 
prior to the filling is required to reduce the pipe wall temperature to the saturation temperature of the 
cryogenic liquid. 

Cryogenic chill-down involves complicated hydrodynamic and thermal interactions among liquid, 
vapor, and solid pipe wall. There exist few basic experimental studies and modeling efforts for chill- 
down of cryogenic fluids. Studies on cryogenic chill-down started in 1960’s accompanying the 
development of rocket launching system. Early experimental studies were conducted by Burke et al. [1], 
Graham [2], Bronson et al. [3], Chi and Vetere [4], Steward [5] among other researches. Bronson et al. 
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[3] studied the flow regime in a horizontal pipe during the chill-down by liquid hydrogen. The results 
revealed that the stratified flow is prevalent in the cryogenic chill-down. 

Flow regimes and heat transfer regimes in the horizontal pipe chill-down were also studied by Chi and 
Vetere [4]. Information of flow regimes was deduced by studying the fluid temperature and volume 
fraction during the chill-down. Several flow regimes were identified: single phase vapor, mist flow, slug 
flow, annular flow, bubbly flow, and single-phase liquid flow. Heat transfer regimes were identified as: * 

single-phase vapor convection, film boiling, nucleate boiling, and single-phase liquid convection. 

Recently, Velat et al. [6] systematically studied cryogenic chill-down in a horizontal pipe. Their study 
included: visually recording the chill-down in a transparent Pyrex pipe, which is used to identify the 
flow regime and heat transfer regime; collecting temperature histories at different positions of wall in 
the chill-down; recording pressure drop along the pipe. Chung et al. [7] conducted a similar study on the 
nitrogen chill-down at relatively low mass flux and provided the data needed to assess various heat * 
transfer coefficients in the present study. 

Burke et al. [1] developed a crude chill-down model based on the one-dimensional heat transfer through 
the pipe wall and the assumption of infinite heat transfer rate from the cryogenic fluid to the pipe wall. 

The effects of flow regimes on the heat transfer rate were neglected. Graham et al. [2] correlated heat 
transfer coefficient and pressure drop with the Martinelli number [8] based on their experimental data. 

Chi [9] developed a one-dimensional model for energy equations of the liquid and the wall, based on 
film boiling heat transfer between the wall and fluid. An empirical equation for predicting chill-down 
time and temperature was proposed. 

Steward [5] developed a homogeneous flow model for the chill-down. The model treated the cryogenic 
fluid as a homogeneous mixture. The continuity, momentum and energy equations of mixture were 
solved to obtain density, pressure and temperature of mixture. Various heat transfer regimes were 
considered: film boiling, nucleate boiling and single-phase convection heat transfer. Separate treatment 
of different heat transfer regimes resulted in a significant improvement in the prediction of the chill- 
down process. Homogeneous mixture model was also employed by Cross et al. [10] who obtained a 
correlation for the wall temperature in the chill-down with an oversimplified heat transfer model of the 
heat transfer between the wall and the fluid. 

Stratified flow regime, which is the prevalent flow regime in horizontal chill-down, was first studied by 
Chen and Baneijee [11-13]. They developed a separated flow model for the simulating cool-down by a 
stratified flow in a hot horizontal pipe. Both phases were modeled using one-dimensional mass and 
momentum conservation equations. The wall temperature was computed using a 2-dimensional transient 
heat conduction equation. Their prediction for wall temperature agreed well with their experimental 
results. Although a significant progress was made on handling the momentum equations, the heat 
transfer correlations employed were not as advanced. 
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Figure 1. Schematic of chill-down heat transfer regimes in a horizontal pipe. 

Typical chill-down involves several heat transfer regimes as shown in Fig. 1. Near the liquid front is the 
film boiling. The knowledge of the heat transfer in the film boiling regime is relatively limited, because 
i) film boiling has not been the central interest in industrial applications; and ii) high temperature 
difference causes the difficulties in experimental investigations. For the film boiling on vertical surfaces, 
early work was reported by Bromley [14], Dougall and Rohsenow [15] and Laverty and Rohsenow [16]. 
Film boiling in a horizontal cylinder was first studied by Bromley [17]; and the Bromley correlation was 
widely used. Breen and Westwater [18] modified Bromley’s equation to account for very small tubes 
and large tubes. If the tube is larger than wavelength associated with the Taylor instability, the heat 
transfer correlation is reduced to Berenson’s correlation [19] for a horizontal surface. 

Empirical correlations for film boiling were proposed by Hendrick et al. [20] [21], Ellerbrock et al. [22], 
von Glahn [23], Giarratano and Smith [24], These correlations relate a simple or modified Nusselt 
number ratio to the Martinelli parameter. Giarratano and Smith [24] gave detailed assessment on these 
correlations. All these correlations are for steady state cryogenic film boiling. They may not be suited 
for transient chill-down application. 

When the pipe chills down further, film boiling ceases and transient boiling occurs, followed by nucleate 
boiling. The heat transfer in the transient boiling is even more complicated and it is usually assumed that 
the boiling switches from film boiling to nucleate boiling right away. The position of film boiling 
transitioning to the nucleate boiling is often called rewetting front, because from that position the cold 
liquid starts touching the pipe wall. Usually the Leidenfrost temperature indicates the transition from the 
film boiling to the nucleate boiling. 

The study on the convection nucleate boiling is extensive. A general correlation for saturated boiling 
was introduced by Chen [25]. Gungor and Winterton [26] modified Chen’s correlation and extend it to 
subcooled boiling. Enhancement and suppression factors for macro-convective heat transfer were 
introduced. Gunger and Winterton’ s correlation can fit experimental data better than the modified 
Chen’s correlation [27] and Stephan and Auracher correlation [28]. Kutateladze [29] and Steiner [30] 
also provided correlations for cryogenic fluids in pool boiling and forced convection boiling. Although 
they are not widely used, they should be more applicable for cryogenic fluids as the correlation was 
directly based on cryogenic conditions. 
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As the wall temperature drops further, the nucleate boiling is replaced by pure convection. The 
convection heat transfer can be modeled using Dittus-Boelter equation [31] for the fully developed 
turbulent pipe flow or corresponding laminar heat transfer equation for laminar pipe flow [31]. The 
condition in which nucleate boiling switches to single-phase convection heat transfer is that all nucleate 
sites are suppressed. 

Although two-fluid model can describe the fluid dynamics aspect of the chill-down process, it suffers 
from computational instability for moderately values of slip velocity between two phases which limits 
its application. To gain fundamental insight into the thermal interaction between the wall and the 
cryogenic fluid and to be able to rapidly predict chill-down in a long pipe, an alternative pseudo-steady 
model is developed. In this model, liquid wave front speed is assumed to be constant and is the same as 
the bulk liquid speed [32]. It is also assumed that steady state thermal fields for both the liquid and the 
solid exist in a reference frame that is moving with the wave front. The heat transfer between the fluid 
and the wall is modeled using different heat transfer correlations depending on the operating heat 
transfer regime at a given location. Various improvements on the correlations are introduced, including 
the development of a new film boiling heat transfer coefficient. The governing equation for the solid 
thermal field becomes a parabolic equation that can be efficiently solved. It must be emphasized that a 
great advantage of the pseudo-steady model is that one can assess the efficacy of the film boiling model 
independently from that of the nucleate boiling model since the down stream information in the nucleate 
boiling regime cannot affect the temperature in the film boiling regime. In another word, even if the 
nucleate boiling heat transfer coefficient is inadequate, the film boiling heat transfer coefficient can still 
be assessed in the film boiling regime by comparing with the measure temperature for the right period of 
time. After the satisfactory performance is achieved for the film boiling regime, the nucleate boiling heat 
transfer model can be subsequently assessed. In the Results section, those detailed assessments of the 
heat transfer coefficients are provided by comparing the computed temperature variations with the 
experimental measurements of Chung et al. [7]. Satisfactory results are obtained. 

NOMENCLATURE 
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Ri and R 2 inner and outer radius of pipe 
Ra Rayleigh number 


Boiling number 
solid heat capacity 
heat capacity 
pipe diameter 
thickness of pipe wall 
gravity 

heat transfer coefficient 

pool boiling heat transfer coefficient 

convection boiling heat transfer coefficient 

latent heat 

Jacob number 

thermal conductivity 

effective thermal conductivity 

Nusselt number 

pressure 

radius of pipe 
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Subscripts 

0 

i and o 

1 

v 

w 

sat 

Superscripts 


Reynolds number 
Peclet number 
Prandtl number 
suppression factor 
temperature 

Leidenffost temperature 

transition temperature between nucleate boiling to convection heat transfer 

room temperature 

time 

velocity 

vapor film velocity 
vapor film coordinates 
transformed coordinate 
cylindrical coordinates 
liquid volume fraction 
Martinelli number 
vapor film thickness 
emissivity 

azimuthal coordinate 
viscosity 

dimensionless temperature; azimuthal coordinate 
liquid surface tension; Stefan Boltzmann constant 

characteristic value 

inner and outer pipe 

liquid 

vapor 

wall 

saturated 

dimensionless variable 


FORMULATION 


In pseudo-steady chill-down model, it is assumed that both the liquid and its wave front moves at a 
constant speed U which is taken from estimated experimental condition in the present model. Thus the 
main emphasis of the present study is on the modeling of the heat transfer coefficients in different heat 
transfer regimes and the computation of the thermal field within the solid pipe. 


Solid Heat Transfer 

The thermal field inside the solid wall is governed by the three-dimensional unsteady energy equation: 
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Since the wave front speed U is assumed to be a constant, it can be expected that when the front is 
reasonably far from the entrance region of the pipe, the thermal field in the solid is in a steady state 
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when it is viewed in the reference frame that moves with the wave front. Thus, the following coordinate 
transformation is introduced, 

Z = z + Ut. (2) 



Figure 2. Coordinate systems: laboratory frame is denoted by z; moving frame is denoted by Z. 


Eq. (1) is then transformed to: 
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For further simplification, the following dimensionless parameters are introduced, 
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where k 0 is the characteristic thermal conductivity, and c 0 is the characteristic heat capacity. Eq. (3) can 
be normalized as 
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where p c = is the Peclet number. 

K 

Under typical operating condition for the cryogenic chill-down process, Pc ~ 0(10 ) - 0(10 ). The first 
term on the RHS of Eq. (5) is small compared with the rest of terms and thus can be neglected. Eq. (5) 
becomes 
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( 6 ) 


which is parabolic. Hence in the Z' -direction, only one condition is needed. In the ^direction, periodic 
boundary conditions are used. On inner and outer surfaces of the pipe wall, where r=R/ and R 2 , proper 
boundary conditions for the temperature are required. 


For convenience, Z' =0 is set at the liquid wave front. In the region of Z' <0, the inner wall is exposed to 
the pure vapor. Although there may be some liquid droplets in the vapor that can cause evaporative 
cooling when the droplets deposit on the wall and the cold flowing vapor can absorb some heat from the 
wall, the heat flux due to these two mechanisms is small compared with the heat transfer between liquid 
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and solid wall in the region of Z’>0. Hence, heat transfer for Z’<0 is neglected and it is assumed that 
9 = 1 at z' = 0 . The computation starts from the Z'=0 to Z' — » oo , until a steady state solution in the Z' - 
direction is reached. An implicit scheme in the Z' - direction is employed to solve Eq. (6). 


Liquid and Vapor Flow 

The two-phase flow is assumed to be stratified as was observed in [7], Both liquid and vapor phases are 
assumed to be at the saturated state. Liquid volume fraction determining which part of wall is in contact 
with the liquid or the vapor, is specified at every cross-section along the Z' -direction based on 
experimental information. The visual studies [6] [7] show that the liquid volume fraction increases 
gradually, rather than abruptly, near the liquid wave front and becomes almost constant during most of 
the chill-down. Hence, the following liquid volume fraction a as a function of time is used in the 
computation of the solid- fluid heat transfer coefficient, 


V^o 


t <t« 


( 7 ) 


a = a 0 t>t 0 

where t 0 is characteristic chill-down time, and a 0 is characteristic liquid volume fraction. Here the time 

when the nucleate boiling is almost suppressed and slope of the wall temperature profile becomes flat is 
set as characteristic chill-down time. 


The vapor phase velocity is assumed to be a constant. However, it was not directly measured. It is 
computationally determined by trial-and-error by fitting the computed and measured wall temperature 
variations for numerous positions. 

Heat Transfer between Cryogenic Fluid and Solid Pipe Wall 

During the chill-down, the fluid in contact with the pipe wall is either liquid or vapor. The mechanisms 
of heat transfer between the liquid and the wall and between the vapor and the wall are different. Based 
on experimental measurements and theoretical analysis, liquid-solid heat transfer accounts for a majority 
of the total heat transfer. However, the prediction for heat transfer is much more complicated than the 
heat transfer between the vapor and the wall. The heat transfer between the liquid and the wall is 
discussed first. 

Heat Transfer between Liquid and Solid wall 

Heat transfer between liquid and solid wall includes film boiling, nucleate boiling, and single-phase 
convection heat transfer. The transition from one type of boiling to another depends on many 
parameters, such as wall temperature, wall heat flux, and various properties of fluid. For simplicity, a 
fixed temperature approach is adopted to determine the transition points. That is, if the wall temperature 
is higher than the Leidenfrost temperature, film boiling dominates. If the wall temperature is between 
Leidenfrost temperature and a transition temperature, T 2 , nucleate boiling takes place. Here, the transient * 
boiling regime between film boiling and nucleate boiling is neglected. The reason is mainly the 
difficulty associated with the determination of the conditions for the transient boiling to occur and the 
lack of a reliable transient heat transfer correlation. If the wall temperature is below the transition 
temperature T 2 , pure convection heat transfer dominates. The values of Leidenfrost temperature and 
transition temperature are determined by matching model prediction with the experimental results. 
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Film boiling heat transfer 

Because of the nature of high wall superheat in chill-down, the film boiling plays a major role in chill- 
down in terms of the time span and in terms of the total amount of heat removed from the wall. 
Currently there exists no specific film boiling correlation for chill-down applications with such high 
superheat. If one uses conventional film boiling correlations, necessary modifications for cryogenic 
application must be made for chill-down. 


One of cryogenic film boiling heat transfer correlations was provided by Giarratano and Smith [24] 

I *Bo~ 0A =f( Zll ) (8) 


f Nu ^ 


\ Nu calcJ 


where Nu, is the Nusselt number for the forced convection heat transfer. In this correlation, the heat 


transfer coefficient is the averaged value for the whole cross section. Similar correlations for cryogenic 
film boiling also exist in literature. The correlations were obtained from measurements conducted under 
steady states. The problem with the use of this steady state film boiling correlation is that it does not take 
into account the change of flow regime as encountered in the chill-down. For example, for the same 
quality, the heat transfer rate in annular flow is much different from that in stratified flow, while those 
empirical correlations cannot take such difference into account. 


Furthermore, in this study, local heat transfer coefficient is needed in order to incorporate the thermal 
interaction with the pipe wall. Since the two-phase flow regime information is available through 
visualization in the present study, the modeling effort needs to take into account the knowledge of the 
flow regimes. 


There are several correlations for the film boiling based on the analysis of vapor film boundary layer and 
stability of the thin vapor film, such as Bromley’s correlation [17] and Breen and Westerwater’s 
correlation [18] for film boiling on the outer surface of a hot tube, Frederking and Clark’s [33] and 
Carey’s [34] correlations for film boiling on the surface of a sphere. However, none of these was 
obtained for cryogenic fluids or for film boiling on the inner surface of a pipe. 

A new correlation for the film boiling in cryogenic chill-down inside a tube is presented here. The 
schematic of film boiling inside a pipe is shown in Fig. 3 with a cross-sectional view. Bulk liquid is near 
the bottom of the pipe. Beneath the liquid is a thin vapor film. Due to buoyancy force, the vapor in the 
film flows upward along the azimuthal direction. Heat is transferred through the thin vapor film from the 
solid to the liquid. Reliable heat transfer correlation for film boiling in pipes or tubes requires the 
knowledge of the thin vapor film thickness which can be obtained by solving the film layer continuity, 
momentum, and energy equations. 
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Figure 3. Schematic diagram of film boiling in stratified flow. 


To simplify the analysis for the vapor film heat transfer, it is assumed the liquid velocity in the 
azimuthal direction is zero and the vapor flow in the direction perpendicular to the cross-section is 
negligible. It is further assumed that the vapor film thickness is small compared with the pipe radius and 
vapor flow is in steady state, incompressible and laminar. The laminar flow assumption can be 
confirmed post priori as the Reynolds number. Re, based on the film velocity and film thickness is 
typically of 0(1 o° ~ 10 2 )- In terms of the x- & y-coordinates and (u, v) velocity components shown in Fig. 


3, the governing equations for the vapor flow are similar to boundary-layer equations: 
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where subscribe v refers to properties of vapor. 


( 11 ) 


Because the length scale in the azimuthal (x) direction is much larger than the length scale at the normal 
(y) direction, the v-component may be neglected. Furthermore, the convection term is assumed small 
and is neglected. The resulting momentum equation is simplified to 
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The vapor pressure can be evaluated by considering the hydrostatic pressure from liquid core as: 
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where (j> is angular position where the film merges with the vapor core. The momentum equation 
becomes 


(P; ~ Pv 
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The vapor velocity boundary condition is u = 0 at y = 0 and u-u, = 0 at y = 8 . The vapor velocity 
profile is: 
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The mean velocity u is 
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The energy and mass balances on the vapor film requires that 
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Neglecting the convection, the vapor energy equation is: 

d 2 T 


dy 2 


= 0 . 


The following linear temperature profile is thus obtained, 
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Substituting the temperature and velocity profiles into Eq. (17) yields 
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Eq. (20) has analytical solution: 

where Ja is Jacob number and Ra is Raleigh number: 
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describes the geometric dependence of vapor film thickness. 
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The mean velocity u as a function of 6 is thus 


u = 


{T.-Tsa, tP,~P,)g R 


F^e) Sin(0). (25) 

J 

Curves for f{o) andF 2 (6>)sin<? based on numerical integration are shown in Fig. 4. The vapor film 
thickness has a minimum at 9 = 0 and is nearly constant for 9 < n / 2 . It rapidly grows after 9 > nil. The 
singularity at the top of tube when G -> n is of no practical significance since the film will merge with 
the vapor core at the vapor- liquid interface. The vapor velocity is controlled by f(g) 2 sin 9 which is zero 
at the bottom of the pipe and increases almost linearly in the lower part of the tube where the vapor film 
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thickness does not change substantially. In the upper part of the tube, due to the increase in the vapor 
film thickness, the vapor velocity gradually drops back to zero at the top of the tube. Thus a maximum 
velocity may exist in the upper part of the tube. 

The local film boiling heat transfer coefficient is easily obtained from the linear temperature profile. 
It is 


h 



0.6389 


K f Ra' 
DF(e){ja y 


4 


(25) 



Figure 4. Numerical solutions of the vapor thickness and velocity influence functions. 

Nucleate boiling and convection heat transfer 

For the nucleate flow boiling heat transfer, Gungor and Winterton’s correlation [26] is widely used due 
to its relatively better accuracy in predicting heat transfer coefficients. However, a closer examination on 
this correlation shows that it is based mainly on the following parameters: Pr, Re, and quality x. Similar 
to the development of conventional film boiling correlations, these parameters all reflect overall 
properties of the flow in the pipe and are not directly related to flow regimes. Thus it cannot be used to 
predict the local heat transfer coefficients in chill-down. 

Chen’s [25] flow boiling correlation is based on separating the heat transfer to micro- and macro- 
convection heat transfer. Micro-convection heat transfer represents the contribution from boiling heat 
transfer, and macro-convection represents the contribution from the forced convection heat transfer. 
However his correlation fits best for annular flow because of the derivation of suppression factor. In 
stratified flow regime, which is common in cryogenic chill-down, Chen’s correlation may not be 
applicable. 

Several correlations have been tried in this study, including Gungor and Winterton’s correlations [26], 
Chen’s correlation [25], and Kutateladze’s correlations [29]. None gives a satisfactory heat transfer rate 
that is need to match the experimentally measured temperature histories in [7] in the nucleate boiling 
regime. Among them, Kutateladze’s correlation gives more reasonable result. In this correlation, the 
total heat transfer coefficient h is 

h = Kony+ h pool > (26) 

where h conv is given by Dittus-Boelter equation for fully developed pipe flow, 
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in which AT is wall superheat. 


(27) 

(28) 


Kutateladze’s correlation [29] was proposed without considering the effect of nucleate site suppression. 
This obviously leads to an overestimation of the nucleate boiling heat transfer rate. Hence a modified 
version of Kutateladze’s correlations (28) is used, 

h = h conv +S*h pool , (29) 

with S being the suppression factor and h pool is given by (28). 


When AT drops to a certain range all the nucleate sites are suppressed. The heat transfer is dominated by 
single phase forced convection. The heat transfer coefficient can then be predicated using Dittus-Boelter 
equation, Eq. (27), if flow is turbulent, or Eq. (30), if flow is laminar. 

Kom ~ 4.36* A:, ID. (30) 

Heat Transfer between Vapor and Solid Wall 

The heat transfer between the vapor and the wall can be estimated by treating the flow as a fully 
developed convection flow, neglecting the liquid droplets that are entrapped in the vapor. The heat 
transfer coefficient of vapor forced convection flow is 

h v = 0.023 *Re“- 8 Pr v 0 - 4 *£„/£> (turbulent flow) (31) 

or 

h v = 4.36 * kJD (laminar flow) (32) 

Heat Transfer between Solid Wall and Environment 

For a cryogenic flow facility, although serious insulation is applied, the heat leakage to environment is 
still considerable due to the large temperature difference between the cryogenic fluid and the 
environment. It is necessary to evaluate the heat leakage from the inner pipe to environment in order to 
make realistic assessment of the model prediction with the experimental results [7]. 

Vacuum insulation chamber between the inner and outer pipes is used in cryogenic transport pipe [7], as 
shown in Fig. 5. Radiation heat transfer exists between the inner and outer pipe. Furthermore, the space 
between the inner and outer pipe is not an absolute vacuum. There is residual air that causes free 
convection between the inner and outer pipes driven by the temperature difference of the inner and outer 
pipe. 

The radiation between the inner pipe and outer pipe becomes significant when the inner pipe is chilled 
down. The heat transfer coefficient is proportional to the difference of the fourth power of wall 
temperatures. Exact evaluation of the radiation heat transfer between the inner and outer pipe is a 
difficult task. Hence a simplified model based on the overall radiation heat transfer between long 
concentric cylinders with constant temperature at inner pipe and outer pipe [31] is used to evaluate the 
heat transfer rate at every axial location of the pipe. It is not quantitatively correct, but can provide 
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reasonable estimate for the magnitude of the radiation heat transfer between pipes through the vacuum. 
The local radiation heat transfer rate per unit area on the surface of inner pipe q' rad is 



Figure 5. Schematic of vacuum insulation chamber. 
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where the cr is Stefan Boltzmann constant, (r„ %) and (r 0 , s 0 ) are the radius and emissivity of inner pipe 
and outer pipe, respectively, and T wa ii is the local inner wall temperature, T 0 is the room temperature that 
is assumed constant in the entire outer pipe. Here the emissivity is also assumed to be constant during 
the entire chill-down process. 


For the free convection heat transfer due to the residual air in 
pipe and outer pipe, Raithby and Hollands’ correlation [35] 
average heat transfer rate per unit length of the cylinder is 


2 nk 


q frc 


eff 
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the vacuum chamber between the inner 
is used for the heat transfer rate. The 

(34) 


where T is assumed constant on the inner and outer wall along the azimuthal directions, and k eff is the 
effective thermal conductivity. 


RESULTS AND DISCUSSIONS 

In experiment by Chung et al. [7], liquid nitrogen was used. The flow regime is revealed to be stratified 
flow by visual observation, as shown in Fig. 7, and wall temperature history in several azimuthal 
positions is measured by thermal couples and recorded on a computer. 

Experiment of Chung et al. 

In experiment of Chung et al. [7], a concentric pipe test section (Fig.6) was used. The chamber between 
the inner and outer pipe is vacuumed, but about 20% air remained. The inner diameter (I.D.) and outer 
diameter (O.D.) of inner pipe are 11.1 and 15.9 mm, and I.D. and O.D. of outer pipe are 95.3 and 
101.6mm, respectively. Numerous thermal couples were placed at different locations of the inner pipe. 
Some were embedded very closely to the inner surface of the inner pipe while others measure the 
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outside wall temperature of inner pipe. Experiments were carried out at room temperature and 
atmosphere pressure. Liquid nitrogen flows from a reservoir to the test section driven by gravity. As 
liquid nitrogen flows through the pipe, it evaporates and chills the pipe. Some of the typical visual 
results are shown in Fig. 7. The measured average liquid nitrogen velocity is U~5 cm/s. Vapor velocity 
is not measured in this experiment. In this study, it is determined through trial-and-error by fitting the 
computed and measured temperature histories. The characteristic liquid volume fraction is 0.3 from the 
recorded video images. The characteristic time used in this computation is ? 0 = ioo s. The Leidenfrost 
temperature for the nitrogen is around 180 K; hence the temperature when the film boiling ends and 
nucleates boiling starts is set at 180 K. The transition temperature at which the nucleate flow boiling 
switches to purely convection heat transfer is 140 K based on experimental results. The material of the 
inner pipe and outer pipe used in the experiment of [7] are Pyrex glass with emissivity of 0.82 at room 
temperature. 


S'tondord fo 



Figure 6. Schematic of Chung et al.’s cryogenic two-phase flow test apparatus [7]. 



Figure 7. Experimental visual observation of Chung et al.’s cryogenic two-phase flow test [7]. 


Comparison of Pipe Wall Temperature 

In the computation, there are 40 grids along the radial direction and 40 grids along azimuthal direction 
for the inner pipe (Fig. 8). The results of temperature profile at 40X40 grids and higher grids resolution 
show that 40X40 grids are sufficient. Figures 9-10 compare the measured and computed wall 
temperature as a function of time at positions 11, 12, 14 and 15 shown in Fig. 8. For the modified 
Kutateladze’s correlation a suppression factor 0.01 is used. Vapor velocity is 0.5m/s based on the best 
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fit. The overall temperature histories agree well in the film boiling stage. Thus, the film boiling heat 
transfer coefficient based on the first principle and incorporated the flow structure gives very reliable 
prediction. It must also be noted that the value of the Leidenfrost temperature does not affect the 
computed temperature history prior to that transition point since the governing equation is parabolic. 
The good agreement before the Leidenfrost temperature is reached is entirely due to the superior 
performance of the new film boiling heat transfer coefficient. 



Figure 8. The computational grids arrangement and positions of thermal couples. 

During the stage of the rapidly decreasing wall temperature after the Leidenfrost temperature, the 
computed wall temperature drops slightly faster than the measured value. The rapid decrease in the wall 
temperature is due to initiation of nucleate boiling which is more efficient for heat transfer than the film 
boiling. Reasonable agreement between the computed and measured histories in this nucleate boiling 
regime is due to: i) the good agreement already achieved in the film boiling stage; ii) valid choice for the 
Leidenfrost temperature that switches the heat transfer regime correctly; and iii) appropriate 
modification of Kutateladze’s correlations. 

In the final stage of chill-down, the wall temperature decreases slowly, and the computed wall 
temperature shows the same trend as the measured one but tends to be a little lower. Fig. 1 1 shows the 
temperature distribution of a given cross-section at different times of chill-down. Because the upper part 
of pipe wall is exposed to nitrogen vapor, cooling effect is much reduced. 


104 


UNIVERSITY OF 

^FLORIDA 



Figure 9. Comparison of wall temperature of transducer 12 and 15, which is at the bottom of pipe, 
between the experiment and computation. For comparison purpose, the solution in the Z-direction is 
converted to time t. 



Figure 10. Comparison of wall temperature of transducer 1 1 and 14, which are at the upper part of the 
pipe, between the experiment and computation. For comparison purpose, the solution in the Z-direction 
is converted to time t. 
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Figure 11. Wall temperature distribution across the section at t=0, 50, 100 and 300 seconds. 

Discussions and Remarks 

In Fig. 9, the wall temperature based on the film boiling correlation of Giarratano and Smith [24] is also 
shown. Apparently, the correlation of Giarratano and Smith [24] gives a very low heat transfer rate so 
that the wall temperature remains high. This comparison confirms our earlier argument that correlations 
based on the overall flow parameter, such as quality and averaged Reynolds number, are not applicable 
for the simulation of unsteady chill-down. 

The nucleate flow boiling correlations of Gungor and Winterton [26], Chen [25], Kutateladze’s 
correlations [29] are also compared in this study. Gungor and Winterton’ s correlation fails to give a 
converged heat transfer rate during the transition from the film boiling to the nucleate boiling. Chen’s 
correlation overestimates the heat transfer rate, and cause unrealistically large temperature drop on the 
wall, which results in the halt of the computation. Only Kutateladze’s correlation gives a reasonable heat 
transfer rate. However, the temperature drop near the bottom of the pipe is still faster than the measured 
one as shown in Fig. 9. This may be due to the fact that most of nucleate boiling correlations were 
obtained from experiments of low wall superheat. However, in this cryogenic chill-down, wall superheat 
is much higher than that in the normal nucleate boiling experiments. Another reason is that the original 
Kutateladze’s correlation does not include suppression factor. This leads to overestimating heat transfer 
coefficient. The modified Kutateladze’s correlation with suppression factor gives reasonable chill-down 
result in Fig.9. 

Further examination of Fig. 9 and Fig. 10 indicates that although we have considered the heat leak from 
the outer wall to the inner wall through radiation and free convection, the computed temperature is still 
lower than the measured temperature during the final stage of chill-down, where the heat transfer rate 
between fluid and wall is low due to lower wall superheat and possibly heat leak. The temperature 
difference between the computed and measured values at position 12 and position 15 suggests that there 
may be additional heat leak, which affects the measurements but is not taken in account in the present 
modeling. 

In this study, the pseudo-steady chill-down model is developed to predict the chill-down process in a 
horizontal pipe in stratified flow regime. This model can also be extended to describe annular flow chill- 
down in the horizontal or the vertical pipe with minor changes on the boundary condition for solid 
temperature. It can also be extended to study the chill-down in the slug flow as long as we specify the 
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contact period between the solid and the liquid or the vapor. The disadvantage of the current pseudo- 
steady chill-down model is that the fluid interaction inside the pipe is largely neglected and both vapor 
and liquid velocities are assumed to be constant. Compared with a more complete model that 
incorporates the two-fluid model, the present pseudo-steady chill-down model requires more 
experimental measurements as inputs. However, the pseudo-steady chill-down model is computationally 
more robust and efficient for predicting chill-down process. It provides overall reasonable results for the 
solid wall temperature. While a more complete model for chill-down process that incorporates the mass, 
momentum, and energy equations of vapor and liquid is being developed to reduce the dependence of 
the experimental inputs for the liquid velocity and trial-and-error for the vapor velocity, the present 
study has revealed useful insights into the key elements of two-phase heat transfer encountered in the 
chill-down process which have been largely ignored. It also laid the necessary modeling foundation for 
the incorporation of the two-fluid model . 


CONCLUSIONS 

A pseudo-steady chill-down computational model has been developed to understand the heat transfer 
mechanisms of cryogenic chill-down and predict the chill-down wall temperature history in a horizontal 
pipeline. The model assumes the constant speed of the moving liquid wave front, and steady state of the 
thermal field in the solid in a reference frame that moves with the liquid wave front, and saturated 
temperature of both liquid and vapor so that the 3-D unsteady problem can be transformed to a 2-D, 
parabolic problem. A new film boiling heat transfer coefficient in the cryogenic chill-down condition is 
developed using the first principle and incorporating the stratified flow structure. The existing nucleate 
boiling heat transfer correlations may not work well under cryogenic condition. A modified 
Kutateladze’s correlations with suppression factor adequately describes heat transfer coefficient. With 
the new and modified heat transfer correlations, the pipe wall temperature history based on the pseudo- 
steady chill-down model matches well with the experimental results by Chung et al. [7] for almost the 
entire chill-down process. The pseudo-steady chill-down model has captured the important features of 
thermal interaction between the pipe wall and the cryogenic fluid. 
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PROJECT TITLE: COMPOSITE PRESSURE VESSELS FOR CTS 

Task PI: Dr. Bhavani Sankar, Department of Mechanical & Aerospace Engineering 
Co-I: Dr. Peter Ifju, Department of Mechanical & Aerospace Engineering 

PROJECT GOALS 

Honeycomb sandwich panels find a wide range of structural applications due to their high 
strength and stiffness-to-weight ratio compared to other material systems. Structural 
honeycombs are composed of two thin, stiff, strong sheets serving as the primary load carrying 
elements and a thick layer of low-density cellular material, providing shear resistance and 
stiffness. Current use of honeycomb material ranges from secondary structure in commercial 
aircraft to primary structure in military aircraft, helicopters, and even space shuttles. New 
challenges to damage tolerance emerge as more frequent use of this material tales place in 
advanced aerospace applications. 

As the cost of space operations increase, new designs for reusable launch vehicles have come to 
be light as a necessary solution for future space missions. In this research, interfacial mixed- 
mode fracture of the face-sheet and the core is evaluated both experimentally and analytically in 
environments similar to those encountered in cryo-tanks used in reusable launch vehicles. 
Interfacial fracture toughness tests are performed at near liquid nitrogen temperature as well as 
room temperature to determine temperature effects on the critical strain energy release rate of the 
sandwich system. The tests will be modeled analytically using ABAQUS finite element software 
to check agreement between analytical and experimental results, and develop a predictive 
capability for future RLV designs. In addition to studying the interfacial fracture behavior of 
sandwich composites, it is necessary to understand the micro cracking of face sheets that may 
cause hydrogen leakage. 

The micro cracking behavior of the liquid hydrogen composite tank at cryogenic temperature has 
been investigated using the micromechanics methods. The micro cracks on composite laminates 
cause the leakage of liquid hydrogen because thermal stress in the composite laminates increases 
due to the difference of the coefficients of thermal expansion (CTE) between fiber and matrix at 
cryogenic temperatures. The micro cracks normal to the interface between composite laminates 
are propagated repeatedly without rupturing the fiber. The interfacial failure between orthotropic 
materials due to micro cracking is investigated. The stress singularity near the micro crack 
perpendicular to an interface between two dissimilar orthotropic materials can be obtained using 
the Stroh’s method. The finite element analysis was performed to verify the analytical result 
under identical boundary condition. In the experimental approach, the four-point bending test is 
conducted with carbon composite laminate. 

The singularity result from the analytical approach can be integrated into the experimental 
approach to predict the interfacial fracture behavior near crack-tip. Furthermore, the relationship 
between crack density and strain energy release rate can be obtained using the finite element 
method. The results will be used to suggest efficient and safe designs for composite hydrogen 
storage systems. Analysis of a hydrogen storage system under more realistic loading conditions, 
modeling progressive damage, and permeability modeling will be conducted in the future. 
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ACCOMPLISHMENTS 

ABSTRACT 

A micromechanics method is developed to investigate microcrack propagation in a liquid 
hydrogen composite tank at cryogenic temperature. The unit cell is modeled using square and 
hexagonal shapes depends on fiber and matrix layout from microscopic images of composite 
laminates. Periodic boundary conditions are applied to the unit cell. The temperature dependent 
properties are taken into account in the analysis. The laminate properties estimated by the 
micromechanics method are compared with empirical solutions using constituent properties. The 
micro stresses in the fiber and matrix phases based on boundary conditions in laminate level are 
calculated to predict the formation of microcracks in the matrix. The method is applied to an 
actual liquid hydrogen storage system. The analysis predicts micro stresses in the matrix phase 
are large enough to cause microcracks in the composite. 

Stress singularity of a transverse crack normal to a ply-interface is investigated to predict the 
fracture behavior at cryogenic conditions using analytical and finite element analysis. When a 
transverse crack touches a ply-interface of a composite layer with same fiber orientation, the 
stress singularity is equal to Vz. When the transverse crack propagates to a stiffer layer normal to 
the ply-direction, the singularity becomes less than Vz and vice versa. Finite element analysis is 
performed to predict the fracture toughness of a laminated beam subjected to fracture loads 
measured by four-point bending tests at room and cryogenic temperatures. As results, the 
fracture load at cryogenic temperature is significantly lower than that at room temperature. 
However, when thermal stresses are taken into consideration, for both cases of room and 
cryogenic temperatures, the difference of the fracture toughness becomes insignificant. The 
result indicates fracture toughness is a characteristic property, which is independent to 
temperature changes. 

The experimental analysis is performed to investigate the effect of cryogenic cycling on 
permeability for various composite material systems. Textile composites have lower 
permeability than laminated composites even with increasing number of cryogenic cycle. Nano- 
particles dispersed in laminated composites do not show improvement on permeability. The 
optical inspection is performed to investigate the microcrack propagation and void content in 
laminated composites and compared the microscopic results before and after cryogenic cycling 

THE FULL REPORT FOR THIS SECTION IS ATTACHED 
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TITLE: ZERO BOIL-OFF (ZBO) PRESSURE CONTROL 

Task PI: Dr. Jacob N. Chung, Department of Mechanical & Aerospace Engineering 

PROJECT GOALS 

A unique facet of storing liquid cryogens is that a sizable amount of heat will continually flow 
from the surrounding ambient to the low temperature fluid. Despite the fact that excellent 
multilayer insulations have been developed for storing cryogens, there will always be a finite 
amount of heat flow into the storage vessels, especially in the vicinity of supporting structures. 

As the heat accumulates, it will eventually vaporize a finite amount of liquid. Due to the large 
volumetric expansion from vaporization, the pressure in the vessel will significantly increase 
unless the vapor is exhausted from the vessel. Thus, there exists continual purging of excess 
vapor from the cryogenic storage tank. This is not efficient, and depending on the duration of 
storage, a significant amount of cryogens can be wasted. In order to preserve the cryogens in 
storage for a long duration space flight, an active heat removal system is required to balance the 
heat that permeates into the storage vessel. Such a system would eliminate the need for purging 
the cryogen vapor from the storage tank. We expect to develop a ZBO system model that would 
determine whether active mixing of the cryogen is required for ZBO and whether the benefits 
outweigh the complications. Next the model would provide a comparison study between a 
sensible heat transport system and a latent heat transport system. We also will design an optimal 
latent heat transport system and determine the appropriate degree of sub cooling. 

ACCOMPLISHMENTS 

ABSTRACT 

This report documents the research work performed under the NASA hydrogen funding for the 
past two years. This project was focused on the evaluation of the feasibility of a cryogenic 
storage system based on phase-change heat transfer processes in microgravity. The basic 
mechanisms were designed following the fundamental concept of zero-boil-off pressure control. 
To maintain zero-boil-off, the heat transfer from the tank depends on bubble removal from the 
tank wall, entrainment of the bubbles by the swirling jet flow and eventual total condensation of 
the entrained bubbles by the cryocooler. The conceptual design was evaluated by scaled-down 
laboratory models, mathematical models and numerical simulations. 

I. INTRODUCTION 

A unique facet of storing liquid cryogens is that a sizable amount of heat will continually flow 
from the surrounding ambient to the low temperature fluid. Despite the fact that excellent 
multilayer insulations have been developed for storing cryogens, there will always exist a finite 
amount of heat flow into the storage vessels, especially in the vicinity of supporting structures. 
As the heat accumulates, it will eventually vaporize a finite amount of liquid. Due to the large 
volumetric expansion from vaporization, the pressure in the vessel will significantly increase 
unless the vapor is exhausted from the vessel. Thus there exists continual purging of excess 
vapor from the cryogenic storage tank. This is not efficient, and depending on the duration of 
storage, a significant amount of cryogens can be wasted. In order to preserve the cryogens in 
storage for a long duration space flight, an active heat removal system is required to balance the 
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heat that permeates into the storage vessel. Such a system would eliminate the need for purging 
the cryogen vapor from the storage tank. We expect to develop a zero-boil-off (ZBO) system 
model with phase-change heat transfer that would determine whether active mixing of the 
cryogen is required for ZBO and whether the benefits outweigh the complications. Next the 
model would provide a comparison study between a sensible heat transport system and a latent 
heat transport system. We also will design an optimal latent heat transport system and determine 
the appropriate degree of subcooling. 

Since the main mechanism for controlling the boil off relies on the total condensation of the 
vapor generated as a result of the heat leakage through the tank insulation, it is very important to 
make sure that the mixing jet is properly designed such that it can efficiently entrain the vapor 
and mix it with the subcooled liquid for a complete condensation. In a terrestrial condition, the 
vapor bubbles would detach from the tank wall after they grow to a certain size and rise up to the 
top of the tank due to gravity where all the bubbles would gather and coalesce to form an 
overhead vapor void. So mixing jets are designed to break up the large vapor chunks and then 
mix the vapor fragments with the subcooled liquid for condensation. In microgravity, however, 
the vapor bubbles would tend to grow and stick around the tank surface where they are nucleated 
because of the surface tension dominance in the absence of buoyancy in space. Therefore, in 
microgravity the mixing jets must be properly designed such that they can first shear off the 
vapor bubbles from the wall, entrain them, and then mix them with the subcooled liquid for 
condensation. Please see Figure 1 for a schematic that explains this proposed process in 
microgravity. Based on the microgravity scenario, two tasks of the ZBO project are given as 
follows : 

Task 1 : To investigate the removal (shear-off) characteristics of a vapor blob from the heater 
surface by a parallel flow which simulates a wall jet during the jet mixer operation. The expected 
outcome is a relationship between the threshold jet flow rate and the wall heat flux for bubble 
removal from the tank wall. 

Task 2 : To study the bubble entrainment capability of a swirling jet flow, vapor bubble 
transport in the jet flow, and condensation rate of vapor bubble. 

For the microgravity experiment, we will use the currently available 1.8-second drop tower at the 
University of Florida. The 1.8-second drop tower is believed to be able to provide sufficient 
preliminary data for the vapor blob shear-off experiment because that the formation of bubble 
blobs of a size around 1-2 cm only takes about less than 0.7 second under certain conditions and 
it only needs less than 1 second for the planned slowest flow to travel through the test section. 
According to a report submitted to NASA by Merte, Lee and Keller (1995), their microgravity 
experiment of pool boiling with R-l 13 on board Space Shuttle STS 60 showed that bubbles grew 
to a size of 5 mm in just 0.05 second under the heat flux of 2 W/cm 2 . Our own experiment also 
revealed a similar phenomenon. An experimental apparatus for performing Task 1 has been 
competed and calibration and testing of the loop has started. Figure 2a is a schematic of the 
experimental apparatus and Figure 3a is a photograph of the apparatus. The apparatus is 
composed of a flow loop with the velocity maintained by a variable-speed pump. The square test 
section is made with transparent plexiglass for visualization and image acquisition by CCD and 
high-speed camera. A downward-facing gold-film heater is located on the top surface of the test 
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section. The downward-facing heater is designed to simulate zero-gravity in a terrestrial 
condition. Justification is given in the Appendix. For Task 2, the experimental apparatus is under 
design and construction. A schematic of that is shown in Figure 2b. Figure 3b shows the design 
of Task 2 apparatus. The Task 2 apparatus is intended to study the entrainment and condensation 
of bubbles by the jet flow. 


II. EXPERIMENTAL APPARATUS 


A flow controlled loop was built for Task 1 as shown in Figure 3a. The experimental fluid is 
refrigerant-113. As shown in the Appendix, scaling analysis based on three dimensionless 
parameters (Ri, Bo, We) proves that the flow loop built can be used to model cryogenic 
convective heat transfer in microgravity conditions for some special cases. A gold-film heater 
has been used to model heat leaks in cryogenic propellant tank in space. Based on the scaling 
analyses, the heater was placed in an inverted horizontal position. Upon bubble formation on the 
heater surface, the pump, used to drive the fluid, would be turned on to the desired flow rate. 



Stream Lines 


Figure 1 . Jet mixer interaction with the vapor blobs. 
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Figure 2b. Schematic of Task 2 experimental apparatus 
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Figure 3a. Photograph of the Task 1 experiment 



Figure 3b. Design of Task 2 experimental apparatus 
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The test section, where the heater was placed, was made of clear polycarbonate. A high speed 
CCD camera was used to view the formation and departure of the bubbles. The flow loop is 
approximately 16"X17"X20" and designed to expand into NASA's standard drop tower cubicles' 
dimensions. 

FLOW LOOP COMPONENTS 
Test Section 

The flow loop's test section, shown in Figure 4 was made of clear 0.5" thick polycarbonate. The 
test section is 8"X3"X2". A 2"X0.25" deep square hole was cut on the test section to place the 
base of the heater. Below the square hole, a 1.75" diameter whole by 0.25" deep was cut to place 
the heater, which is attached to the square base, and an O-ring, to prevent leaks. The test section 
was attached at the seams with leak proof adhesive that could withstand the fluid's boiling 
temperature of 48 °C and high impacts, for experiments in the drop tower. 



Figure 4. Loop’s Test Section Drawing. 


Pump 

The pump used is the following: In-Line Circulator Pump - Centrifugal Type: Wet Rotor Model 
Number: Astro 30B Flange Size: .75,1,1.25,1.5 In. Horsepower: 1/25 Voltage: 115 Phase: 1 
Housing Material: Bronze; Flow rate 0-20 GPM. Figure 5 shows the pump’s picture. 
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Figure 5. Pump. 


Heater 

The heater used in the loop is a 1"X1"X0.125" piece of polycarbonate with a thin gold film layer 
on top. The heater was attached to a 1.75" diameter by 0.25" deep piece of polycarbonate that 
served as the base for the heater. The round base of the heater was connected to a 2" square 
(0.25" thick) piece of polycarbonate in order to connect and secure the base of the heater to the 
test section. Four 0.125" diameter holes were drilled through the square piece for power and 
voltage leads. 

The gold film was used to transfer heat from a power supply to the fluid via the lead brass 
screws. The gold film was sputtered on the polycarbonate piece to a thickness of 400 Angstroms 
(40 nm) in the center and 5000 Angstroms (500 nm) on the leads of the screws. The Hummer 
VII sputtering system with etch/plate control was used to sputter the gold film to the 
polycarbonate section. Cathode sputtering was the method used with the Hummer VII. This 
method uses high voltage (>10,000 volts) across two electrodes. An ion bombardment vaporizes 
the cathode metal and deposits it as a crystalline, fine-grained coating on material placed near the 
anode. To connect the gold film surface to the lead screws, silver filled epoxy was applied from 
the top of the lead screws to the surface of the heater. The heater was then baked under 1 70 °F 
for two hours to allow for the epoxy to set. 

Pipes 

Standard 1" NPT-PVC pipes were used to connect the loop together. Teflon tape was used at the 
connections to prevent leakage under moderate pressure. 

Fluid 

The fluid used is refrigerant-113 (l,l,2-trichloro-l,2,2-trifluoroethane 99.9%) used to simulate 
cryogenic fluids, specifically liquid hydrogen, at near saturation conditions. The fluid was 
purchased from ChemNet in Sarasota, Florida. The trade name for R-113 is 
Trichlorotrifluoroethane. Below are some specifications for R-l 13: 


Molecular Weight 187.38 
Boiling Point 1 1 8°F (48°C) 

Vapor Pressure 6.46 Psia at 77°F (25°C) 
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Liquid Density 

1,51 1 kg/m 3 

Vapor Density 

3.48 kg/m 3 

Surface Tension 

15.9E-03 N/m 

Viscosity 

497E-06 N-s/m 2 


Heater Calibration 

Based on the research of Oker and Merte (1981), it was found that there is a linear relationship 
between the heater resistance and surface temperature. It was also found that small shifts in 
calibration would happen with time, but the slope of the line would remain the same. To verify 
Oker and Merte's results, a calibration was done on our heaters. The gold film heater was 
calibrated over a temperature range of 18°C to 80°C. The calibration was performed with a 
double-heating system. The inner container size is 6”X6”X6” and is made of 0.25” thick 
aluminum plate with 6.125” thick and 1” fins protruding inward for maximum heat transfer. 
There were two 0.5” fins that protrude out of the bottom to prevent direct contact with the outer 
container. A Coming PC-35 heater was used to heat up the outer container and the water 
between the two containers. The inner container was filled with Dow Coming 200 fluid. The 
heater calibration assembly contains a base for the heater and the thermocouple was immersed in 
the fluid. The heat output from the PC-35 heater transfers through the outer container and water 
uniformly into the inner container. 

Data was taken when the heater was turned off and the temperature and resistance were steady, 
the process was repeated until all the required data was taken. The temperature data was taken 
with a type T thermocouple and a digital thermometer. The resistance data was taken using a 
Keithley 197 Autoranging Microvolt DMM with a double-ohm method to get rid of the 
connecting resistance of the leads. The current leads of the gold-film heater were connected to 
the input side of the DMM, and the voltage leads were connected to the Q sense sides of the 
DMM, as shown in Figure 6. 


Figure 6. Heater Calibration Schematic. 
The calibration curve is shown in Figure 7 with a linear relation: 



T=3 1 7.3295*R-894. 1 226 


L 
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Figure 7. Heater Calibration Curve. 


Because of the small shifts in calibration, a single point calibration for finding the y-intercept is 
necessary before each experiment. 

Video System 

To observe the bubbles coming from the heater surface, a video system was setup. The video 
system consisted of a high speed CCD camera (CV-255E) with a Vz n HAD Interline Transfer 
CCD sensor. The camera’s electronic shutter speed was from 1/1,000 to 1/10,000 sec. A Ceptre 
PD-1212ARD5 (12 VDC 1.2 A) power supply was used to provide power to the camera. A 
piece of Teflon was used as insulation to prevent the heat of the lights to affect the fluid. The 
Teflon also provided a better contrast background to see the bubbles. Figure 8 shows the 
schematic. 
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Power Supplies 

A Hewlett Packard power supply (model # 6263B) provided electrical power to the heater with 
variable voltage input between 0-20V. The power supply to the camera and pump was a 
standard 110 VAC outlets. 


III. EXPERIMENTAL PROCESS 

Several procedures needed to be performed before beginning each experiment. First, the single 
point heater calibration must be completed. Second, the desired heat flux must be set via the 
power heater’s power supply. Finally, the desired flow rate is chosen. 

Single Point Calibration 

A single point calibration was performed before each experiment. A digital thermometer was 
connected to the lower thermocouple in the test section to take the temperature data. A double- 
ohm method was used to eliminate the connection resistance. 


Temperature and resistance reading were taken once the temperature and resistance were 
stabilized. The y-intercept was calculated based on the linear relation, 


T, = ^R, +b 
r dR r 


where T r is the temperature reading, R r is the resistance reading. The slope, dt/dR, is known 
from the calibration of the heater. This calibration procedure has to be done before each 
experiment because the y-intercept varies between experiments. 


Heat Flux 

The heat flux was obtained by varying the voltage on the power supply to the heater. The heat 
fluxes chosen for the experiment varied between 1 W/cm 2 to 3 W/cm 2 . The heat flux from the 
heater was calculated from the measured values of voltage and current by Equation 2, 


q 


tl 


V-I 

A 


where V is voltage, I is current, and A is the heater’s surface area (1” X 1”). The current reading 
was calculated from the measured values of voltage and current by Equation 3, 

K 

where V s is the voltage across the shunt resistance and R s is the resistance of the shunt resister. 


Flow System 

According to the numerical analyses’ results, the flow rates for the experiments conducted were 
between 0-2 USGPM. These flow rates yielded maximum heat transfer to condensate the 
bubbles, minimum kinetic energy from the fluid to the bubbles, and allowed for microgravity 
simulation in earth’s gravity. The pump was turned on and allowed to reach steady state; 
approximately in 2 minutes. Several flow rates were used for each heat flux chosen. 
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Experimental Procedure 

The experimental procedure is summarized as follows: 

1. Setup required software (LabView) and hardware (power supply, heater, pump, and CCD 
camera). 

2. Select the desired heat flux. 

3. Perform the single point calibration 

4. Select the desired flow rate. 

5. Allow for steady state to be reached. 

6. Activate LabView to record data. 

7. Turn off all equipment once data has been collected. 

Data Reduction 

The raw data recorded by LabView was copied into spreadsheet. The heater temperature was 
converted from resistance to temperature by the calibration curve. The bulk liquid temperature 
was measure with type-T thermocouples and a cold junction compensator. The voltage reading 
was converted to temperature with the aid of a reference table provided by Omega. After all data 
was converted to readable data, curves were plotted and results were analyzed. 

IV. EXPERIMENTAL RESULTS 
Flow Visualization 

Flow visualization results have been obtained by our research group for three heating levels and 
one flow rate. The flow rate corresponds to a bulk velocity of 7.7 cm/s. Both the 1-g and 
microgravity results are shown together. The no-flow and with flow photographs are displayed 
side by side with each other for a direct comparison. Figures 9, 10 and 11 show the results for a 
heat flux level of 2.9 W/cm 2 . For a higher heating rate, Figures 12 and 13 provides visualization 
photographs for a heat flux level of 4.6 W/cm 2 . Figure 14 and 15 represent the highest heating 
rate of 7.6 W/cm 2 . As seen from these photographic results, the velocity of 7.7 cm/s is able to 
sweep off about 30% of the bubbles for the 2.9 W/cm 2 case. No bubble was removed for both 4.6 
W/cm 2 and 7.6 W/cm 2 cases. It is also noted that the sizes of the bubbles increase with the 
heating rate. It is confirmed that in microgravity, the generated bubble size can be as large as 1 
cm after only 0.067 second for the lowest heating level of 2.9 W/cm 2 . 
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Figure 9. Heat flux level of 2.9 W/cm 2 
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Velocity = 0 cm/s 


Velocity = 7.7 cm/s 


Figure 10. Heat flux level of 2.9 W/cm 
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Velocity = 0 cm/s 


Velocity = 7.7 cm/s 


Figure 1 1 . Heat flux level of 2.9 W/cm' 
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Figure 13. Heat flux level of 4.6 W/cm 2 
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Figure 14. Heat flux level of 7.6 W/cm 
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Figure 15. Heat flux level of 7.6 W/cm' 
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The Heater Surface Roughness Effect 

The effects of the heater surface roughness on the bubble removal are presented next. Two 
heaters with different roughness levels are shown in Figures 16a and 16b. Heater #2 is rougher 
than heater #1. Figure 17 gives the threshold values for flow rates and heat fluxes for the two 
heaters. It is noted that the rougher the heater surface, the more power is available for the surface 
to hold on to the bubble. Therefore, the threshold flow rate is higher for the rougher surface at 
the same heat flux. For the threshold heat flux, a higher flow rate is needed for the rougher 
surface. While for the threshold heat flux, at the same flow rate, the heat flux is higher for the 
less rough surface. 



Figure 16a. Rough heater #1. 



Figure 16b. Rough heater #2. 
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Figure 17. Threshold flow rates and heat fluxes. 

V. THERMODYNAMIC AND MATHEMATICAL MODELING 
System Model 

In order to evaluate the performance and optimize the compact storage module, a model that 
correctly describes the actual physical process is needed. Figure 18 shows the schematic of the 
model that corresponds directly to the compact storage system shown in Figure 2b. For the 
model, it is seen that the heat input to the system is represented by a heater (yellow plate) and the 
heat sink (cryocooler in blue) is the condenser. 


Inflow: 77K Outflow 
! 0.25m/s 


Heater VA 


Model 


Cryocooler 70K 


Adiabatic 
Walls 





Figure 18. Schematic of the Model 
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The thermodynamic and heat transfer diagram for the model is given in Figure 19 where the 
system is composed of an evaporator, a condenser and a jet pump. Further more, the heat transfer 
processes in the evaporator and the condenser are given in Figure 19 using electrical equivalent 
circuits, respectively. 



Figure 19. Thermodynamic and Mathematical Model 
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The evaporator represents the boiling process on the wall that uses the heat entering the storage 
tank from outside through the insulation. There are three heat transfer resistances involved. They 
are convection resistance on the outer surface, conduction resistance of the wall insulation and 
boiling resistance on the inner surface. For the condenser, there are two resistances, one with the 
condensation heat transfer and the other with the cryocooler. 

The governing equations are given below : 



evap 


(T„-T m ) 

1 L 1 

1 1 

K- h evap 


Q x is evaporation rate, 


A evap is evaporator surface area. 



cond 


(K-T c ) 

1 1 

1 _ — 

h cond 


Q\ ~ Qi 


Q 2 is condensation rate, 


A cond is condenser surface area. 


Comparison Study 

For the sample calculation, nitrogen is used as the storage liquid. The insulation layer is 10 cm 
thick Kapton with thermal conductivity of 0.2 W/mK. The ambient temperature is set at 300K. 
The jet flow rate is assumed to be O.lkg/s. The cryocooler heat transfer coefficient is estimated to 
be at 10 W/m2K. 

The overall heat transfer coefficients for the evaporator and the condenser are presented in 
Figure 20 with single-phase and two-phase flow regimes indicated in the figures. Whereas the 
effectiveness values for the evaporator and the condenser are plotted in Figure 21 . 


{ 
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Figure 20. Overall heat transfer coefficients for the evaporator and the condenser. 
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Figure 21. Heat exchanger effectivenesses for the evaporator and the condenser. 
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The performance evaluation for the entire system is illustrated in Figures 22 and 23 where the 
overall system effectiveness as a heat exchanger and the cryocooler size are provided. We see 
that the average effectiveness is 0.2 for the single-phase system while it is around 0.24 for the 
two-phase system. The improvement is approximately 20%. On the cryocooler size study, the 
average is 1.4 m for the single-phase system and is 1.15 m for the two-phase system, 
respectively. The reduction in the cryocooler size is also 20% for the two-phase system. 

Numerical Simulation 

Based Figure 18, numerical simulations using the commercial code FLUENT were also 
performed. The results are plotted in Figure 23 for the flow streamlines and in Figure 24 for the 
bubble trajectories, respectively. The streamline pattern basically verified the design objective 
that the flow would form a strong recirculation around the cryocooler for bubble entrainment and 
condensation. In Figure 24, in addition to the bubble trajectories, the plot also shows the bubble 
temperatures. In these preliminary calculations, we can see that the bubbles are cooled down as 
they are entrained in the recirculation. 



Figure 22-a. Overall effectiveness of the system. 
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Figure 22-b. Cryocooler size 
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Figure 23. Flow streamlines. 
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Figure 24. Bubble trajectories and their temperatures. 
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APPENDIX: SCALING ANALYSIS 

In addition to the drop tower work, we also plan to run terrestrial experiments that can simulate 
the microgravity condition. The following is a study that provides criteria on how to perform 
terrestrial experiments that would simulate microgravity conditions. 

Basis for simulating microgravity with Terrestrial Experiment 

Scaling Analysis for Convective Boiling Heat Transfer 

Flow boiling is a complicated subject because that it is involved with both convection and 
boiling. First we would investigate the forces that are present in the flow boiling process. The 
flow boiling behavior is governed by the inertia force of the convection flow and the dynamics 
behavior of the vapor generated through nucleation and vaporization, which itself is a 


138 



UNIVERSITY OF 



consequence of the forces acting at the interface between the liquid and vapor phases. These 
would be more precisely expressed in terms of local force balance. 


Flow Boiling over a Flat Plate (applicable to vapor shear-off experiment for ZBO) 

For the current evaluation, we would only consider the global nature of the forces as suggested 
by Merte et al. (Criteria for approximating certain microgravity flow boiling characteristics in 
earth gravity. Proceedings of the Microgravity Transport Processes in Fluid, Thermal, Biological 
and Materials Sciences, Banff, Alberta, Canada, 2001) 


The potential forces are listed below : 

I. Buoyancy - arising from gravity due to density difference between the vapor and 
liquid. 

II. Momentum associated with the liquid flow that is imposed on the vapor bubbles in the 
form of drags. 


III. Interfacial (surface tension) - considered here as liquid/vapor and liquid/solid/vapor. 

The former Primarily controls the liquid/vapor surface geometry, while the latter 
controls the adhesion or departure of the vapor bubble from the solid heater surface. 

VI. Viscous shear- arising principally from the liquid viscosity. 

V. Lift- arising from the externally imposed liquid flow around the vapor bubble and 
it acts normal to the flow. 

VI. Internal vapor pressure- generated by the mean liquid/vapor interface 
temperature, provides the driving force for vapor bubble growth /collapse dynamics. 

For the purpose of assessing the degree to which flow boiling over a flat surface in earth gravity 
may approximate that in microgravity, advantage will be taken of ratios of certain of the above 
forces, in particular those associated with the removal of the vapor bubble from its nucleation 
site on the heater in the inverted horizontal or nearly horizontal position relative to earth gravity. 
Since the bulk liquid flow is parallel to the heater surface, only those components of the forces 
that are parallel to the heater surface will be considered, including buoyancy. Because of the 
absence of terrestrial gravity in space, the normal component of buoyancy as a vapor bubble 
removal mechanism is eliminated by the inversion of the heater such that the removal of bubble 
would be effected by other forces. However, the influence of lift, which is also normal to the 
heater surface, must be given consideration, and could be of special significance in microgravity. 
For spherical particles moving in an inviscid fluid parallel to a plane wall, it has been 
demonstrated that the lift moves the sphere toward the wall. With a viscous fluid, on the other 
hand, the resulting velocity gradient in duct flow causes the sphere to migrate away from the 
stationary wall. This away lift force is important for bubble removal in microgravity. In earth 
gravity, buoyancy acting to push the vapor bubble towards the wall for the case considered here, 
order of magnitude analysis shows that such away lift can be neglected relative to the buoyancy 
for the conditions involved in the current experiment. In addition, the internal vapor pressure 
giving rise to the bubble dynamics will be neglected for present purposes because of its minor 
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effects as compared to other forces. This leaves the forces related to buoyancy, momentum (or 
inertial drag), surface tension, and viscous shear. 

Several dimensionless numbers, defined in terms of ratios of these forces acting on the vapor 
bubbles, will now be examined to estimate the conditions under which certain aspects of the flow 
boiling process in microgravity might be approximated in earth gravity. These are: 


Richardson number, Ri, = Buoyancy/Inertia Force 
Weber number, We, = Inertia Forced/Surface Tension 
Bond number. Bo, = Buoyancy/Surface Tension 
Capillary number, Ca, = Viscous Force/Surface Tension 
Reynolds number, Re, = Inertia Force/Viscous Force 
Ca = We/Re 

The Capillary number can be expressed as the ratio of the Weber number to the Reynolds 
number: If the Viscous Drag in the Capillary number were to be replaced by the Inertial Drag, 
which would correspond to a Reynolds number of unity, the Capillary number would become 
equivalent to the Weber number. Since lift is being neglected, the Capillary number will be given 
no further consideration here. Because of the involvement of both vapor and liquid phases, a 
two-phase Richardson number, Ri 2lj) , has been defined in which the buoyant and drag forces 

acting on a vapor bubble are defined in terms of the bubble radius, the relative velocity between 
the bubble and the bulk liquid, U re i, and coefficients, Cb and Cd, given as: 

F b = c b (All)xR 3 g e kp (Al) 

F d = c d xR 2 p ( {U 2 re i / 2) (A2) 

The buoyancy coefficient, Cb, accounts for the difference between a spherical bubble of radius, 
R, and that of a non spherical bubble of equal mean radius, and Cd is the drag coefficient, 
determined empirically. The drag force can be considered to be a manifestation of the net bulk 
flow momentum effects acting on the vapor bubble. The characteristic bubble radius, R, in Eqs. 
(Al) and (A2) is taken to be the departure bubble radius Rd- In the above, g e is the effective 
gravity depending on the inclination angle between the heater and the gravity. A p is the 
difference between the liquid and vapor densities. p t is the liquid density and U re i is the relative 
velocity between the liquid and vapor bubble. If the bubble is stationary with the heater, then 
U rel is equal to the liquid velocity, U t . Once the vapor bubble has departed its local nucleation 

site, it can be considered to no longer directly influence the local heat transfer processes. An 
expression for the departure radius was obtained from a force balance on a vapor bubble growing 
from its nucleation site, utilizing Eqs. (Al) and (A2), along with a surface tension force, F s , at 
the solid/liquid/vapor contact line, as: 

F s =C s <jR (A3) 
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The surface tension coefficient, C s , was obtained empirically by Kirk and Merte (A study of the 
relative effects of buoyancy and liquid momentum in forced convective nucleate boiling. Final 
Report, NASA Grant NAG3-1310, 1992) from the balance of buoyancy and surface tension 
forces for the case of pool boiling by equating Eq. (A3) to Eq. (Al) with Cb=l. Using the 
measured departure radius of 0.50mm for R113, the surface tension coefficient was found to be 
approximately unity. 

An estimate of the drag Coefficient, Cd, was determined aslo by Kirk and Merte (1992) 
empirically from a balance of the buoyant and drag forces, Eqs. (Al) and (A2), parallel to a 
heating surface at which an immediately adjacent spherical bubble was moving. Measurements 
of the bubble rise velocities, sizes and liquid velocities were made for various orientations of the 
flat heating surface relative to the horizontal face-up position, such that g e , in Eq. (Al) is 
replaced by: 

g e = gSinO (A4) 


♦ 


Where g is earth gravity and 6 is the inclination angle (0° heater facing upward and 1 80° heater 
facing downward). Over a range of bubble Reynolds numbers from 24 to 920, Kirk and Merte 
(1992) found that the drag Coefficient, Cd, is roughly constant and equal to 0.56. This compares 
well with the turbulent drag coefficient of 0.50 for a solid sphere. 

Taking the ratio of Eqs. (Al) and (A2) as the two-phase Richardson number, together with the 
coefficients given above: 

R ht= F B /F D ( A5 ) 


R h</> ~ 


8gC b (Ap/p e )R 


3 c U 
^d^rel 


(A6) 


Assuming that both buoyancy and drag act to detach the bubble and the surface tension tends to 
hold the bubble, then at the bubble departure, we can write : 


F b + F d -F s = o 


(A7) 


Next we solve Eq. (A7) to obtain the departure radius, R<j, as follows : 


R _ ( 3C d PfUh ¥ L . 64C,C b o&pgj 


l6C d Apg 


2>nC d p f U 


rel 


(A8) 


Substituting Rd for R in Eq. (A6), we obtain the following 


^=(i/2)(Jl + 64C ^ < f^ -l) 
* 1 3<plC 


(A9) 
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The best estimates for Q,, C s and Cd are 1.0, 1.0 and 0.56, respectively based on most boiling 
conditions as discussed above. With these constants inserted, Eq. (A9) becomes : 


Ri 2 s = (1 / 2 )(. |1 + 


2\.loApg e 

Pforel 


- 1 ) 


(A10) 


By the same approach, we can show that the two-phase Weber number (inertia drag/surface 
tension), We 2<jl and the Bond number (buoyancy/surface tension), Bo as follows : 


We u =(0.1/ M)(Vl + 21.7M - 1) 

5o = 0.42(21.7 -4^.) 

M 


(All) 
(A 12) 


Where M = 8e f P ° . 

P(Urel 

The liquid velocity U t is taken for U rel and g e is given by Eq. (A4). 

Eqs. (A10), (All) and (A12) are plotted in Figures Al, A2 and A3 for R-113, nitrogen and 
hydrogen at atmospheric pressure, for three orientation angles of 90°, 135° and 175° and the 
range of velocities used here are 0 to 25 cm/s. 

Outcome of the Scaling Analysis 

For the inclination angle of G = 180°, the heater is upside down and facing downward, we have 
the following : 

Rht 0 
We 2 ^\ 

Bo ->■ 0 

Therefore, when both the Ri 2(j> and the Bo are zero the gravity effects are excluded. For 
convective flow boiling, if we intend to simulate a microgravity condition in earth gravity we 
only have to focus on the Ri 2(/) and the Bo because that the buoyancy force is only associated with 

these two numbers. As for the We 2tj> it only gives an indication on the extent of bubble shape 

deviation from a sphere. It is generally agreed that a bubble can be approximated as a sphere if 
the We 2rj> is less than four. 

For the inclination angle of 180°, theoretically the buoyancy is totally irrelevant. But in actual 
experiments, the gravity would still play some role due to the following reasons : 

1. The low-level vibrations and shaking of the laboratory floor that are inherent in most 
buildings would produce buoyancy force purterbations to the boiling experiments. 


142 


49& UNIVERSITY OF 

^FLORIDA 

2. The upward buoyancy force that pushes the bubble against the heater surface could cause 
the bubble to flatten and have larger base area on the heater surface. As a result, more 
inertia force might be required to detach the bubble. 

Based on these two possible buoyancy effects, we would use the inclination angle of 175° to 
select the conditions that would allow us to simulate microgravity in earth gravity. In other 
words, the 5° deviation from the upside down (180°) would hopefully represent the two 
buoyancy purterbations mentioned above. 

For our experiments with R-113, we may conclude that if the liquid velocity is greater than 15 
cm/s that would result in both the Ri 2/j) and the Bo less than one which allows us to do earth 
gravity experiment to simulate microgravity conditions. 

The other important information we leam from the scaling analysis is that, for all three numbers, 
R-113 and Nitrogen give almost identical values for the velocity ranges selected. It basically 
verify that we will be able to use the R-113 results to predict the nitrogen phenomena for cases 
that are dominated by the Ri 2 ^ , the Bo and the We 2<f) . Based on the values of the Richardson 
number, the buoyancy is negligible for the liquid velocity higher than 20 cm/s. 



Relative Velocity [m/s] 


Figure A1 . Richardson number vs. flow velocity 
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Figure A2. Weber number vs. flow velocity 
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Figure A3. Bond number vs. flow velocity 
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PROJECT TITLE: FLUID DISTRIBUTION FOR IN-SPACE CRYOGENIC 
PROPULSION 

Task PI: Dr. William Lear, Department of Mechanical and Aerospace Engineering 

PROJECT GOALS 

The ultimate goal of this task is to enable the use of a single supply of cryogenic propellants for 
three distinct spacecraft propulsion missions: main propulsion, orbital maneuvering, and attitude 
control. A fluid distribution system is sought which allows large propellant flows during the first 
two missions while still allowing control of small propellant flows during attitude control. 
Existing research has identified the probable benefits of a combined thermal 
management/power/fluid distribution system based on the Solar Integrated Thermal Management 
and Power (SITMAP) cycle. Both a numerical model and an experimental model are 
constructed in order to predict the performance of such an integrated thermal 
management/propulsion system. This research task provides a numerical model and an 
experimental apparatus which will simulate an integrated thermal/power/fluid management 
system based on the SITMAP cycle, and assess its feasibility for various space missions. Various 
modifications are done to the cycle, such as the addition of a regeneration process that allows 
heat to be transferred into the working fluid prior to the solar collector, thereby reducing the 
collector size and weight. Fabri choking analysis was also accounted for. Finally the cycle is to 
be optimized for various space missions based on a mass based figure of merit, namely the 
System Mass Ratio (SMR). 

The theoretical and experimental results from these models are be used to develop a design code 
(JETSIT code) which is able to provide design parameters for such a system, over a range of 
cooling loads, power generation, and attitude control thrust levels. The performance gains and 
mass savings will be compared to those of existing spacecraft systems. 

OVERVIEW 

After an extensive review of fluid distribution technology, it was decided that the most 
promising avenue of research for this project would involve a combined thermal 
management/power generation/ propellant distribution system that is based on the SITMAP 
cycle. The SITMAP cycle integrates a vapor compression (refrigeration) thermal management 
subsystem, which uses a jet pump rather than a compressor, and a power subsystem that is 
essentially a Rankine cycle in which the jet pump replaces a turbine. In the current 
implementation, the working fluid would be one or more cryogenic propellants. The SITMAP 
cycle’s physical setup within a spacecraft could provide a constant flow of propellant to the 
spacecraft’s extremities, and this flow could easily be tapped into at an ideal location in order to 
bleed off propellant for attitude control purposes. 

This leads to another benefit of the SITMAP cycle for fluid distribution purposes. Because of 
the SITMAP cycle operational characteristics; there is a region in the cycle where pressurized 
superheated vapor exists - after the waste heat recovery unit and/or solar collector and before the 
turbine inlet. This part of the cycle could easily be tapped and the superheated vapor used 
directly as a thrust source for attitude control. 
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Using the fluid in its highest energy state within the SITMAP cycle would optimize its 
effectiveness as a propellant and could eliminate the need for any type of burner, resulting in a 
simple and more reliable means of spacecraft attitude control. Figure 1 illustrates such a system. 

\ Qrad 


Pump 



Figure 1: SITMAP-based integrated thermal management/power generation/ propellant 
delivery system. 


A parallel University of Florida research project has focused on the numerical and experimental 
analysis of a SITMAP cycle that does not include a mass flow bleed. The results of this parallel 
research have been and will continue to be utilized in the development of the modified SITMAP 
cycle (with a mass flow bleed, or bleeds). Both computer codes and experimental apparatus 
have been constructed which will be progressively modified for use in the analysis of a modified 
SITMAP cycle. 


Research work done on the SITMAP cycle is comprised of two parts, an analytical part, and an 
experimental part. Following are summaries of both parts of the program. 


ANALYTICAL PROGRAM 

The analytical part went through different stages starting with the development of a computer 
code (JetSit, which is short for Jet-pump and SITMAP) for the thermodynamic simulation of the 
cycle. An expression was then needed to measure the mass based performance of the cycle. For 
that purpose a System Mass Ratio (SMR) expression was derived and incorporated in the 
simulation code JetSit. Further development of the code involved incorporating thermodynamic 
properties software to calculate thermodynamic properties instead of the use of a data file. The 
software used is called REFPROP and is made by NIST. Recuperation was then included in the 
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SITMAP cycle as an effort to decrease the weight of the system (see Figure 2). Fabri choking is 
an important phenomenon that might take place in the mixing chamber of the Jet-pump. Fabri 
choking refers to conditions when the primary flow expands in the mixing chamber constricting 
the available flow area for the secondary stream, causing it to accelerate. It is possible for the 
secondary stream to reach sonic velocity, therefore causing the secondary mass flow rate to 
become independent of downstream conditions. The Fabri choking analysis was included in the 
JetSit cycle simulation code to make sure that the input entrainment ratios are physically 
possible. 



Figure 2. Schematic of the Solar Integrated Thermal Management and Power (SITMAP) 
cycle with regeneration 

The current and most recent stage of the analytical part in the SITMAP program is the 
optimization of the cycle to minimize the mass for different space missions. To achieve this, an 
optimization program was incorporated in the cycle simulation code. The optimization routine is 
written by Dr. Leon Lasdon of the University of Texas in Austin and it utilizes a Generalized 
Reduced Gradient algorithm, and is hence called GRG2. 

The most results of the analytical part of the program are presented below. Appendix A includes 
the more recent work that hasn’t been published yet, including the Fabri choking, and 
optimization analyses. 


Results of the preceding analysis are presented in Tables 1 through 5. Each of the five sets has 

(P . / P\ 

the same primary inlet pressure, while varying the secondary pressure. The ratio v ' V! ,bo is 
the break-off pressure ratio. This is the ratio below which the flow operates either in the mixed 
regime or in the saturated supersonic regime (inlet choking regime), and above which the flow 
operates either in the mixed regime or the supersonic regime (Fabri choking regime). 

There are some general trends that can bee seen in all data sets. First, once the pressure ratio 
Ppi/Psi is larger than its break-off value, cjiFabri becomes the upper limit (supersonic flow 
regime). When Ppi/Psi is less than its break-off value the limiting value for the entrainment ratio 
is cfiinlet choke (saturated supersonic regime). For any line of data in all data sets, if the 
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entrainment ratio drops below the limiting value, the jet-pump will be operating in the “mixed 
regime”, where cj) is dependent on the downstream conditions (or the back pressure). It should 
also be noted that as long as the ratio Ppi/Psi is above the break-off ratio, the pressure ratio 
Pse/Pne drops below unity, which is important for the primary flow to expand into the secondary 
causing Fabri choking to take place. 


The second general trend is the increase of the limiting value of the entrainment ratio with higher 
secondary pressure. This should be expected because since the primary inlet pressure is fixed, a 
higher secondary stagnation pressure corresponds to a lower backpressure. The lower 
backpressure allows for more secondary flow entrainment before choking occurs. 


Data sets 1 to 3 show the effect of the ratio of primary and secondary stagnation temperatures on 
the maximum entrainment ratio. Figure 3 is a graphical representation of Tables 1 to 3. It can be 
seen that higher entrainment ratios can be achieved for higher stagnation temperature ratios. The 
effect of Tpi/Tsi on maximum <\> is less significant at higher Ppi/Psi ratios, but the same trend still 
holds. 


Tables 3 and 4 show the effect of varying the primary nozzle geometry. Figure 4 shows the 
same effect in graphical form. It can be seen that lower Ant/Ane ratios (i.e. higher Mne) allow 
more secondary flow entrainment. This makes sense, since the entrainment mechanism is by 
viscous interaction between the secondary and primary streams. Therefore, faster primary flow 
should be able to entrain more secondary flow. 

Tables 4 and 5 show the effect of the area ratio Ane/Ase. Figure 5 shows that lower primary-to- 
secondary area ratios allow for more entrainment. This trend is expected since a lower area ratio 
means more area for the secondary flow and thus more secondary mass flow rate. Figure 6 
shows the variation of the compression ratio with Ppi/Psi, for different jet-pump geometries. The 
compression ratio was calculated for maximum allowable entrainment ratio. It can be seen that 
as the ratio Ppi/Psi increases the compression ratio increases as well, which is expected. 


Table 1: Data set # 1. Tpi = 200, Tsi = 100, Ant/Ane = 0.25, Ane/Ase = 0.1 


Ppi 

Psi 

Ppi / Psi 

PHIinlet 

PHIfabri 

Pse 

Pse/Pne 

(Ppi / Psi)bo 

Pne 

3200000 

1600000 

2 

144.996 


752031.170 

7.30495 

16.51 

102948 

3200000 

1066666.667 

3 

85.465 


766654.157 

7.447 


Mne 

3200000 

640000 

5 

11.1029 


366170.684 

3.55685 


3.07605 

3200000 

320000 

10 

5.58389 


169866.433 

1.65002 


3200000 

213333.3333 

15 

3.67901 


113017.917 

1.09781 


3200000 

160000 

20 


2.683032 

100556.41 

0.97677 


3200000 

128000 

25 


2.083661 

86249.331 

0.83779 


3200000 

106666.6667 

30 


1.676613 

76063.5258 

0.73885 


3200000 

91428.57143 

35 


1 .383796 

68217.1583 

0.66264 


3200000 

80000 

40 


1.164268 

61920.2561 

0.60147 


3200000 

71111.11111 

45 


0.9936368 

56761.0127 

0.55136 


3200000 

64000 

50 


0.8574538 

52445.6223 

0.50944 



Table 2: Data set # 2. Tpi = 100, Tsi = 100, Ant/Ane = 0.25, Ane/Ase = 0.1 
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Ppi 

Psi 

Ppi / Psi 

PHIinlet 

PHIfabri 

Pse 

Pse/Pne 

(Ppi / Psi)bo 

Pne 

3200000 

1600000 

2 

23.2413 


752031.17 

1.8963 

4.62411 

396579 

3200000 

1066666.667 

3 

13.6961 


766654.081 

1.93317 


Mne 

3200000 

533333.3333 

6 


1.416579 

376543.452 

0.94948 


2.57789 

3200000 

320000 

10 


0.6935064 

264842.461 

0.66782 



3200000 

213333.3333 

15 


0.3329346 

195893.323 

0.49396 




Table 3: Data set # 3. Tpi = 400, Tsi = 100, Ant/Ane = 0.25, Ane/Ase = 0.1 


P p, 

Psi 

P./P. 

^ inletchoke 

^ Fabri 

<t> 

P /P 

se / ne 

Pde/Psi 

(Pde/Psi ) bo 

P„e 

3200000 

1600000.000 

2 

213.904 


213.904 

7.968 

0 

18.009699 

94375.5814 

3200000 

640000.000 

5 

16.3609 


16.36093 

3.880 

1.0485 


M ne 

3200000 

320000.000 

10 

8.22779 


8.227792 

1.800 

1.19303 


3200000 

213333.333 

15 

5.42427 


5.424266 

1.198 

1 .34552 


2.957 

3200000 

160000.000 

20 


3.9997 

3.999688 

0.878 

1 .49797 



3200000 

128000.000 

25 


3.1214 

3.121448 

0.773 

1.65079 



3200000 

106666.667 

30 


2.5321 

2.532101 

0.691 

1.80439 



3200000 

91428.571 

35 


2.1111 

2.111098 

0.627 

1.95876 



3200000 

80000.000 

40 


1.7958 

1.795751 

0.574 

2.11383 



3200000 

71111.111 

45 


1 .5509 

1.55091 

0.543 

2.26951 



3200000 

64000.000 

50 


1.359 

1.359 

0.529 

2.427 



Table 4: Data set # z 

. Tpi 

= 400, Tsi = 100, Ant/Ane = 0.6, Ane/Ase = 

= 0.1 


P P, 

P S1 

P./P. 

^ inletchoke 

^ Fabri 

<!> 

P / P 

se / ne 

Pde/Psi 

(Pde/Psi ) bo 

Pne 

3200000 

1600000.000 

2 

89.060 



1.816 

0.637 

4.433 

414056.656 

3200000 

640000.000 

5 


6.738 

6.704 

0.995 

1.187 


M 

ne 

3200000 

320000.000 

10 


3.052 

3.052 

0.575 

1.508 


3200000 

213333.333 

15 


1.810 

1.810 

0.420 

1.842 


1.996 

3200000 

160000.000 

20 


1.194 

1.194 

0.334 

2.183 



3200000 

128000.000 

25 


0.832 

0.832 

0.278 

2.527 



3200000 

106666.667 

30 


0.594 

0.594 

0.239 

2.876 



3200000 

91428.571 

35 


0.428 

0.428 

0.210 

3.228 



3200000 

80000.000 

40 


0.305 

0.305 

0.187 

3.583 



3200000 

71111.111 

45 


0.2109 

0.210895 

0.165 

3.941 



3200000 

64000.000 

50 


0.137 

0.137 

0.153 

4.301 



Table 5: 

Data set # 5. Tpi 

= 400, Tsi = 100, Ant/Ane = 0.25, Ane/Ase 

= 0.3 


Pp, 

P S1 

P./P. 

Pi/ si 

^ inletchoke 

*t ) Fabri 

<t> 

P / P 

se / ne 

Pde/Psi 

(Pde/Psi ) bo 

Pne 

3200000 

1600000.000 

2 

71.301 


71.301 

7.968 

0.715 

18.009699 

94375.581 

3200000 

640000.000 

5 

5.459 


5.459 

3.876 

1.131 



3200000 

320000.000 

10 

2.745 


2.745 

1.800 

1.658 


M ne 

3200000 

213333.333 

15 

1.810 


1.810 

1.198 

2.185 


2.957 

3200000 

160000.000 

20 


1.008 

1.008 

1.112 

2.712 



3200000 

128000.000 

25 


0.789 

0.789 

0.942 

3.108 



3200000 

106666.667 

30 


0.631 

0.631 

0.850 

3.507 



3200000 

91428.571 

35 


0.513 

0.513 

0.772 

3.908 



3200000 

80000.000 

40 


0.421 

0.421 

0.706 

4.311 



3200000 

71111.111 

45 


0.401 

0.401 

0.650 

4.561 



3200000 

64000.000 

50 


0.348 

0.348 

0.602 

4.716 
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Figure 3: Maximum entrainment ratio variation versus the stagnation pressure ratio for different 
stagnation temperatures ratio. Based in (J)inlctchokc and c()Fabri only 



Figure 4: Maximum entrainment ratio variation versus the stagnation pressure ratio for 
different primary nozzle geometry. Based in cjiinlctchokc, and (j)Fabri only. 
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Figure 5: Maximum entrainment ratio variation versus the stagnation pressure ratio for 
different primary to secondary nozzles area ratios. Based in (J)inlctchoke, and 
(jiFabri only. 



15 25 35 45 55 

Ppi / Psi 


Figure 6: Compression ratio variation versus the stagnation pressure ratio 


EXPERIMENTAL PROGRAM 

The initial consideration in the experimental design process was to conceptually develop the best 
method of designing an ejector test rig that will be used to study the phenomena of two-phase 
flow mixing and compression. The primary design criterion established stated that the ejector 
primary and secondary nozzle inlet states were to be varied over a wide range of pressures, 
temperatures, and qualities in order to obtain a general database. The wide range of conditions is 
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vital in thoroughly investigating the effects to two-phase flow. A second design requirement 
dealt with the ability to change the area ratios of the primary and secondary inlet and outlet 
nozzles. The simulation program JetSit aided in the design of the ejector nozzles. 


The first and most important design feature to accommodate was how to vary nozzle inlet states. 
Conceptually, some kind of loop must be designed to heat and cool the primary and secondary 
fluid flow. This heating and cooling process would necessarily incorporate the use of one or 
more heat exchangers. In order to gain control of the pressure, a valve (e.g. adjustable throttling 
valve) would be needed. The working fluid must also be circulated through the loop; 
consequently, a pump would be needed. All of the components listed above are basic parts of a 
system operating on the SITMAP cycle. From this reasoning it was resolved that in order to gain 
relatively good control of the experimental apparatus, all components of that system would need 
to be incorporated. Since it was deemed advantageous to also obtain system-level data for the 
SITMAP cycle, the choice of ejector rig design was driven towards implementing that cycle 
directly. 


Another important design consideration was the selection of a working fluid for the two-phase 
flow ejector test rig. The main criterion used in the fluid selection involved the temperature and 
pressure ranges where the fluid would be in the two-phase region. For safety reasons it was 
decided that the loop would be operated near ambient temperature and pressure. Information 
about individual fluids was found by examining their Material Safety Data Sheets (MSDS). 
Those sheets provide information that cover hazardous ingredients, health effects, handling 
details, fire and explosion hazard data, first-aid procedures, and critical fluid properties. An 
important critical fluid property found in the MSDS is the boiling point of a fluid. Since the 
fluid will be operated in the two-phase region, only fluids with a boiling point temperature 
slightly above ambient temperature were considered. Based on this information, R-141b was 
selected. Once the working fluid was selected, the components were sized. Shown below in 
Figure 7 is a schematic of the final layout of the experimental setup. 



Figure 7. Schematic of the experimental testing apparatus 
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The most recently published paper on the experimental work is presented in appendix B. 


APPENDIX A 


NOMENCLATURE 
Latin Symbols 

A cross-sectional area, m 2 

a speed of sound, m/s 

COP coefficient of performance, dimensionless 

h specific enthalpy, kJ/kg 

M Mach number, dimensionless 

m mass, kg 

m mass flow rate, kg/s 

P pressure, MPa 

P r compression ratio 

Q' heat transfer rate, kW 

r heat exchanger pressure ratio 

T temperature, °C 

V velocity, m/s 

w' work rate, kW 

Greek Symbols 

c() entrainment ratio, m s / rii p 

p density, kg/m 

Subscripts 
de diffuser exit 

ei evaporator inlet 

evap evaporator 

me mixing chamber exit 

n2 primary flow passage where fabri choking occurs 

ne primary nozzle exit 

nt primary nozzle throat 

p primary flow 

pi primary nozzle inlet 

pe pump exit 

pump mechanical pump 

rad radiator 

re radiator exit 

s secondary flow 

s2 secondary flow passage where fabri choking occurs 

sc solar collector 

se secondary flow exit 
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si secondary flow inlet 

ss state immediately downstream of shock wave 

t turbine 

ti turbine inlet 

ts isentropic turbine exit state 


FABRI CHOKING ANALYSIS 


c 


si 


s2 


se 



pi nt ne 


n2 


] 


de 


si 




D 


Figure 1. Schematic for the Jet-pump with constant area mixing, showing the Fabri choked state s2. 


The momentum equation for the control volume over the mixing chamber can be written as 


P se Ae + P ne A ne ~ P s2 A s2 ~ P n2 A n2 = ™ P K 2 + ™s V s 2 ~ ™p V ne ~ ™s V se 


(la) 


dividing by m yields 


~{ P se A se + P ne A ne ~ P s2 A s2 ~ P n2 A n2 )~{K2~Ke) + h Fabri ( V s2 ~ Ke ) 

m p 

. , _ ( P se A se + P ne A ne ~ P s2 A s2 ~ P n2 A i,2 ) IK?! ~ P ne ) 

"* PM = *,(V A -V M ) (mJ 




Fabri 


( A A A A ^ 

p p ^ne _ p 211 2 _ p »2 

se ne . s2 . n2 . . 

V A se A se A ne A se 


o V V -V ) 

r ne r ne ^ \ r s2 r se ) 


(K2~Ke) 

(K 2 -Ke) 


(lb) 

(lc) 


(Id) 


The iteration scheme starts by guessing a value for P se , knowing that s se = s sj , that defines the 
state (se). From the energy equation 
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F,,=[2(*„-A,,)] K 

then <f> Fahn can be calculated as 

l _ PsAseAe 

YFabri ~ ' 


( 2 ) 


(3) 


P nAne A ne 

It should be noted that the area ratio A ne / A se is an input to the SITMAP code. 

Then a guess is made for P s2 , and s s2 = s se , that defines the state (s2). The velocity V s2 can be 
obtained from the energy equation between se and s2 


V s2 = 


( 


V 


2 \ 


h - h 2 + 

se s2 /*> 

v 1 J 


V 


(4) 


calculate M s2 = — , and check if it is equal to 1 . If not another value for P s2 is guessed till 


a, 


s 2 


^ 2 = 1 - 


The area ratio A sJ / A se can be calculated from the continuity equation between se and s2, 


A s2 _ PseKe 
A se Ps?Vs2 


(5) 


knowing that for constant-area mixing A ne + A se = A s2 + A n2 , then 


A n2 _ ^ _j_ A 

A., 


f A A A 

A s 2 A se 

\ Ae Aw J 


( 6 ) 


From equation Id another value for (j) Fuhri can be obtained. Iterate on P se till the values for (f> Fahn 
from equations (Id) and (3) match. 

<j> Fahn is the maximum possible entrainment ratio for a given geometry and inlet states. Therefore 
the simulation code JetSit sets (f> = </> Fabrl , if the input value is greater than <j> Fahri . 


OPTIMIZATION ANALYSIS 

As mentioned before the optimization analysis is performed by incorporating a Generalized 
Reduced Gradient algorithm into the simulation code. The optimization algorithm is written by 
Dr, Leon Lasdon of the University of Texas in Austin, and is called GRG2. 
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All the input variables to the simulation program JetSit are the variables in the optimization 
process. A list of these variables is 


• Jet-pump primary inlet pressure, P p i 

• Jet-pump primary inlet entropy, S pi . 

• Jet-pump secondary inlet pressure, P s i 

• Jet-pump secondary entropy, S sl 

• Jet-pump entrainment ratio, (j) 

• Area ratio Ant/Ane in the jet-pump primary nozzle 

• Area ratio Ane/Ase in the jet-pump 

• Turbine inlet pressure, P u 

The constraints on the optimization process are 

• Pti - Pte >0 

• Ppe - Pre >0 

• Ppi - Psi >0 

• Pjpe - Pse >0 

• Qevap ^ 0 

• Qboiler ^ 0 

• Qcond ^ 0 

• SMR > 0 (objective function, to be minimi z ed) 

The objective function is the SMR, which is to be minimized for space applications. All the 
optimization variables are changed within a specified range to obtain the minimum value for the 
objective function, SMR, provided all the aforementioned constraints are satisfied. 
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ABSTRACT 

The authors are presently involved in developing 
a design code to optimize an active space thermal 
management system that includes as a key 
component an ejector, which operates with fluids in 
the two-phase regime. In order to validate this code, 
and for other applications of two-phase ejectors, a 
comprehensive experimental data set is needed for 
this device. This paper deals with the conceptual 
design and implementation of a constant-area ejector 
experimental rig intended to provide the required 
data set. The system has been designed to 
implement the same thermodynamic cycle as the 
proposed thermal management system, allowing a 
preliminary performance database to be developed 
upon testing, in addition to the ejector data. The 
ejector itself will be an interchangeable part in this 
system, allowing geometrical variables to be 
manipulated. 


INTRODUCTION 

Thermal management of spacecraft relies on 
rejection of heat via radiation, a process that can 
result in large radiator mass, depending upon the 
heat rejection temperature. For some missions, it is 
advantageous to incorporate an active thermal 
management system, allowing the heat rejection 
temperature to be greater than the load temperature. 
This allows a reduction of radiator mass at the 
expense of additional system complexity. A 
particular type of active thermal management system 
is based on a thermodynamic cycle, developed by the 
authors, called the Solar Integrated Thermal 
Management and Power (SITMAP) cycle. This 
system has been a focus of the authors’ research 
program in the recent past (see Fig. 1). The system 
requires no moving parts, which decreases the 
vibration level and enhances reliability. 
Compression of the refrigerant working fluid is 
accomplished in this scheme via an ejector. 
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However, a key uncertainty, addressed in the current 
paper, is the design of the ejector to accommodate 
the two-phase flow, which occurs in many operating 
regimes of interest within the device. 

S Qrad 



Figure 1. Schematic of the Solar Integrated Thermal 
Management and Power (SITMAP) cycle 

Ejectors have been used in a wide variety of 
applications for over a century. Examples of 
applications for ejectors include vacuum pumps in 
the food industry, power stations, and in the 
chemical industry; ejector systems used in the 
aircraft industry for thrust augmentation; and steam- 
jet ejectors used in refrigeration. The ejector, shown 
schematically in Fig. 2, has the advantages of 
simplicity, no moving parts, low cost, and reliability. 
Other advantages of the ejector are decreases in both 
vibration and weight. 

Literature dealing with ejectors is abundant, 
especially when dealing with single-phase flow. 
Single-phase ejectors exhibit a wide range of fluid 
phenomena and have historically incorporated the 
use of empirical methods in their design. However, 
for the case of two-phase flow in ejectors, there is 
little information readily available. More 
information about two-phase ejectors can be found in 
the works of Sherif et al. [1], More research is 
needed in this area to be able to understand and 


design the ejector component, which would enable 
the potential mass savings of the SITMAP system. 



Figure 2. Schematic of ejector (non-constant area mixing) 

The work to be presented in this paper is based 
on research conducted by Nord et al. [2], 
Freudenberg et al. [3], and Kandil et al. [4]. Nord et 
al. [2] developed the combined power and thermal 
management cycle for onboard spacecraft 
applications referred to above (i.e. SITMAP cycle). 
This is essentially an integrated vapor compression 
cycle and a Rankine cycle with the compression 
device being an ejector instead of the regular 
compressor. The power subsystem is a Rankine 
cycle, which drives the system. The ejector acts as 
the joining device between the thermal and power 
subsystems, by mixing the high-pressure flow from 
the power subsystem with the low-pressure flow 
from the refrigeration subsystem providing a 
pressure increase in the latter. Nord et al. [2] used 
Refrigerant 134-a as the working fluid in their 
analysis. The mechanical power produced by the 
turbine can be used to drive the mechanical pump as 
well as other onboard applications. This allows a 
system operating on the SITMAP cycle to be solely 
driven by solar thermal input. 

Examples of work relevant to the SITMAP 
investigation include Bredikhin et al. [5], 
Cunningham and Dopkin [6], Cunningham [7], Eiger 
et al. [8], Fabri and Paulon [9], Fabri and Siestrunk 
[10], Fairuzov and Bredikhin [11], Holladay and 
Hunt [12], Holmes et al. [13], Jiao et al. [14], Lear et 
al. [15], Marini et al. [16], Neve [17]and Sherif et al. 
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[1]. Literature dealing with thermally actuated 
cooling systems includes those by Kakabaev and 
Davletov [18], Chen [19], Lansing and Chai [20], 
and Chai and Lansing [21], Many systems dealing 
with power and thermal management have been 
proposed for which the System Mass Ratio (SMR) 
analysis developed by Freudenberg et al [3] can be 
used, including absorption cooling systems and 
solar-powered vapor jet refrigeration systems. 
Examples of those systems are found in the works of 
Abrahamsson et al. [22], Alefeld and Radermacher 
[23], Anderson [24], Chai and Lansing [21], Chen 
[19], and Lansing and Chai [20]. 


This paper deals with an experimental 
investigation of a novel active thermal management 
system that we have proposed for use in space due to 
its potential for low weight and high reliability. 
Accordingly, this paper will present the design of an 
experimental facility for measuring the performance 
of constant-area ejectors with two-phase flow. 
Significant modeling has been performed in order to 
design this device, and the modeling approach will 
be described, along with the design choices in the 
construction of the facility. The purpose of the 
experimental investigation is two fold: to develop an 
experimental data set for two-phase ejectors (the key 
uncertainty in designing the proposed thermal 
management system) and to develop a proof-of- 
concept experiment for the thermal management 
system itself. The experimental results will be used 
to both calibrate the aforementioned model and also 
provide a high-quality design database of the global 
performance parameters of two-phase ejectors. 

NOMENCLATURE 
Latin Symbols 

2 

A cross-sectional area, m 

D diameter 

h specific enthalpy, kJ/kg 
m mass, kg 

m mass flow rate, kg/s 

P pressure, MPa 

P r compression ratio 

T temperature, °C 

V velocity, m/s 


Greek Symbols 

(j) entrainment ratio, m s / m p 

Subscripts 


de diffuser exit 

me mixing chamber exit 
ne primary nozzle exit 

nt primary nozzle throat 

p primary flow 

pi primary nozzle inlet 

s secondary flow 

se secondary flow exit 

si secondary flow inlet 

EXPERIMENTAL DESIGN 

The initial consideration in the experimental 
design process was to conceptually develop the best 
method of designing an ejector test rig that will be 
used to study the phenomena of two-phase flow 
mixing and compression. The primary design 
criterion established stated that the ejector primary 
and secondary nozzle inlet states were to be varied 
over a wide range of pressures, temperatures, and 
qualities in order to obtain a general database. The 
wide range of conditions is vital in thoroughly 
investigating the effects to two-phase flow. A 
second design requirement dealt with the ability to 
change the area ratios of the primary and secondary 
inlet and outlet nozzles. The computer program 
JetSit, written by Nord et al. [2], aided in the design 
of the ejector nozzles. JetSit, described below, is a 
computer program designed to model the SITMAP 
cycle and calculate the key ejector geometry 
parameters. 

The first and most important design feature to 
accommodate was how to vary nozzle inlet states. 
Conceptually, some kind of loop must be designed to 
heat and cool the primary and secondary fluid flow. 
This heating and cooling process would necessarily 
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incorporate the use of one or more heat exchangers. 
In order to gain control of the pressure, a valve (e.g. 
adjustable throttling valve) would be needed. The 
working fluid must also be circulated through the 
loop; consequently, a pump would be needed. All of 
the components listed above are basic parts of a 
system operating on the SITMAP cycle. From this 
reasoning it was resolved that in order to gain 
relatively good control of the experimental 
apparatus, all components of that system would need 
to be incorporated. Since it was deemed 
advantageous to also obtain system-level data for the 
SITMAP cycle, the choice of ejector rig design was 
driven towards implementing that cycle directly. 

Another important design consideration was the 
selection of a working fluid for the two-phase flow 
ejector test rig. The main criterion used in the fluid 
selection involved the temperature and pressure 
ranges where the fluid would be in the two-phase 
region. For safety reasons it was decided that the 
loop would be operated near ambient temperature 
and pressure. Information about individual fluids 
was found by examining their Material Safety Data 
Sheets (MSDS). Those sheets provide information that 

cover hazardous ingredients, health effects, handling details, fire and 

explosion hazard data, first-aid procedures, and critical fluid 
properties. An important critical fluid property found in the MSDS is the 
boiling point of a fluid. Since the fluid will be operated in the two-phase 
region, only fluids with a boiling point temperature slightly above ambient 

temperature were considered. Based on this information, R- 

141b was selected. 


Experimental Setup 

Once the working fluid was selected, the 
components were sized. Shown below in Fig. 3 is a 
schematic of the final layout of the experimental 
setup. It should be noted that both the boiler and 
evaporator have a controlled, measurable, and 
variable heat input. The heat input to the boiler can 
be varied from 0 to 5 kW, while that to the 
evaporator can be varied from 0 to 3.3 kW. This 
allows the ejector’s primary and secondary inlets to 
be varied over a wide range of temperatures and/or 
qualities. It should also be noted that the pump 
design incorporates a bypass loop. This loop allows 
the pump to be operated at different speeds without 
changing the total mass flow rate through the rest of 
the piping system. Varying the pump speed allows 
the primary nozzle inlet pressure (P p i) to be varied 
directly. The addition of the bypass loop also gives 
us the capability to vary the secondary nozzle inlet 
pressure (P s i) over a larger range. With this 
configuration, P pi can be varied over a range of 0.5 
MPa to 0.83 MPa, while P s i can be varied over the 
range of 0.15 MPa to 0.5 MPa. The condenser 
shown in Fig. 3 is a water-cooled condenser. The 
use of a condenser assures that the refrigerant at the 
pump inlet is always liquid. The charge or amount 
of refrigerant in the system is also adjustable. 
Different charge levels produce different mass flow 
rates and the pressures at all state points except 
during the phase change processes. 
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Pressure Jet Pump Pressure 



The arrangement of the aforementioned 
components used in conjunction with proper 
instrumentation allows the inlet conditions of the 
ejector to be varied over a wide range of conditions. 

EJECTOR 

Once the operating range of the testing apparatus 
was established, the ejector was designed. The 
ejector analysis is not included and the reader is 
referred to Nord et al. [2] for a detailed description 
of that part. Nord et al. [2] presented an iterative 
scheme that solves for the ejector exit state given the 
primary and secondary inlet states and the 
entrainment ratio. They then presented another 
iterative scheme to solve for all the other states in the 
SITMAP cycle given the output of the ejector 
analysis. The JetSit program referred to above was 
developed to efficiently perform the above ejector 
analysis and create performance graphs. JetSit’s 


main functions are to calculate the ejector geometry 
and diffuser exit state and to use those results to 
calculate a solution to the SITMAP cycle. The 
program reads refrigerant properties from a data file, 
performs the necessary interpolations, uses the 
properties in the pump analysis schemes, and writes 
the results to a file. It should be noted that the JetSit 
code assumes constant-pressure mixing, whereas we 
use it for a guide for a constant-area mixing design. 
This is mentioned in more detail in the future work 
section. It should also be mentioned that JetSit is 
capable of analyzing ejector solutions for all flow 
regimes, including saturated flow. 

The operating range of the testing apparatus was 
mapped into the inputs for the JetSit analysis. P s ; 
was varied over the range of 0.15 MPa to 0.5MPa, 
while X s ; was varied from 0.5 to 1.0. The 
entrainment ratio (j) was varied from 0.2 to 2.0. The 
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primary inlet state was always in the superheated 
vapor region. Table 1 below shows some of the 
SITMAP inputs and their corresponding outputs. 
The ejector component outputs are the ejector 
compression ratio, the throat-to-exit area ratio of the 
primary nozzle, and the primary exit-to-secondary 
exit area ratio. 


Table 1. Sample Inputs and Outputs of the 
Ejector Analyses 



Secondary Nozzle Mixin S Chamber 


Figure 4. Ejector design 


SITMAP Inputs 


SITMAP Outputs 


Cas 

T P i 

Pp. 

x si 

Psi 

<t> 

Pr 

A nt /A 

A n e/ A 

1 

13 

.6 

.7 

.15 

1. 

1.65 

.694 

.222 

2 

13 

.6 

1 

.1 

1. 

1.57 

.694 

.194 

3 

13 

.6 

.7 

.1 

.2 

2.7 

.694 

1.223 

4 

13 

.6 

1 

.1 

.2 

2.65 

.694 

1.067 


Once the analysis has been performed over the 
range of specified conditions, the appropriate area 
ratios were determined. 

Shown below in Fig. 4 is a schematic of the 
conceptual design of the ejector. The design is 
centered on the idea of a removable primary nozzle. 
In order to get multiple primary-to-secondary area 
ratios, multiple removable primary nozzles will need 
to be designed. The key to this design is to make it 
where the primary nozzle is easy to get to and easy 
to change. Making the nozzle easy to change will 
speed up the test phase. 


Once the basic ejector geometry has been fixed, the 
removable nozzle dimensions could be designed. 
This part of the design started by first fixing the 
mixing chamber diameter. Next the mixing chamber 
diameter was set equal to the outside diameter of the 
secondary nozzle, while the inside diameter of the 
secondary nozzle was set equal to the outside 
diameter of the primary nozzle. A schematic of this 
is shown in Figure 5. By using the above design 
geometry it is possible to achieve different area 
ratios A nt /A ne and A ne /A se by simply changing the 
primary nozzle. Using this technique, ten nozzles 
were designed to incrementally cover the desired 
range of area ratios. 

Outside Dse=Dme 
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EXPERIMENTAL PROCEDURE 


FUTURE WORK 


An experimental procedure has been developed 
in order to systematically test the ejector over the 
desired range of conditions. The procedure involves 
a test plan using the variable heat inputs along with a 
range of valve settings to achieve specified pressures 
and temperatures. This allows the pump to be tested 
over the widest range of inlet conditions possible. 
The procedure is to be repeated for every nozzle, so 
that the geometrical effects of the primary nozzle 
expansion ratio and the primary-to-secondary area 
ratio may be studied. The data obtained will be 
reduced using a generalized equation of state, so that 
the results will have maximum generality. Thus it is 
expected that the data obtained in this program will 
be applicable with reasonable accuracy to two-phase 
ejectors operating with other working fluids. 

SUMMARY 

An experimental rig has been designed in order 
to test constant-area ejector performance over a 
range of primary and secondary inlet conditions, 
concentrating on operation in the two-phase regime. 
The design allows P p ; to vary from 0.5 MPa to 0.83 
MPa while P s i can be varied over the rage of 0.15 
MPa to 0.5 MPa. Since the heat input to both the 
primary and secondary flow is controlled, the fluid 
for both streams can be liquid, two-phase, or 
superheated vapor. In addition, several ejector 
geometry design variables may be adjusted, allowing 
for multiple ejector area ratios to be tested. All 
support structure and flow components have been 
acquired, and the test facility is in the final stages of 
assembly. The system has been designed to 
implement the SITMAP cycle, allowing a 

preliminary performance database to be developed 
upon testing. The flexibility of this experimental test 
rig will also allow a high-quality design database of 
the global performance parameters of two-phase 
ejectors to be obtained. 


The JetSit code will be modified from its current 
form in order to accommodate certain analytical 
needs. The first change involves matching code 
inputs with experimental inputs. Currently JetSit is a 
design-oriented code and will need to have an 
analysis mode added to it. The current design mode 
only takes the ejector inlet conditions as inputs and 
allows the ejector geometry to automatically adjust 
to its optimum design point. A new analysis mode 
will need to have the ejector geometry as well as the 
inlet state points as inputs. By having the geometry 
as an input a direct comparison between the 
theoretical and actual results for a given geometry 
can be made. 


The second change that will need to be made 
involves the governing equations. Currently the 
JetSit code is based on a constant-pressure mixing 
ejector; whereas, in this experimental setup, the 
ejector is constant-area. The current form of the 
code was used only as a guide to get a range of 
viable area ratios. The main difference in the ejector 
analysis is found in the mixing chamber. In the case 
of constant-pressure mixing P ne = P se ; this must be 
changed in order to account for the case where P„ e ^ 
P se . Using a control volume approach and 
comparing the momentum equation for each case 
help realize the differences in the ejector analysis. 
The constant-pressure control volume as shown in 
Fig. 1 yields Equation (1) and the constant-area 
control volume shown in Figure 5 yields Equation 
(2). In analysis mode, the modified JetSit code will 
implement Equation (2) in the mixing section. 


(F/ ?e f r ) m; ™ p Uc 0 Ue + (l + 4*)^ pV me 

( 1 ) 


{P„e - P me )A me + ( P se - P ne )A se - 


-m p V ne -<f> m p V se + (l + <j>)m p V me 

( 2 ) 
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A series of tests must also be conducted in order 
to characterize the experimental apparatus and to aid 
in the interpretation of the ejector performance data. 
The tests will be performed under certain operating 
conditions in order to get data in the specified 
regions. The results obtained from the tests will be 
used to validate the results from the ejector analysis. 
Shown below in Figs. 6 and 7 are examples of some 
graphical methods that will be used to compare the 
test results to the analytical model. A graph similar 
to Fig. 6 will be used to show the effects that varying 
the secondary inlet will have on the ejector pressure 
ratio for fixed ejector geometry, over a range of 
entrainment ratios. 


3 1 

Psi=.15MPa 

Psi=.3MPa 
Psi=.5MPa 


1 - 



0.5 4 1 , , , 

0 0.5 1 1.5 2 2.5 

<l> 


Figure 6. Performance curves of constant 
A nt /A ne , A ne /A se , T pi , P pi , and X si for various 
design points 


A graph similar to Fig. 7 will be used to show 
how the predicted ejector pressures ratio relates to 
the actual ejector pressure ratio. This will be used 
mainly for the validation of the model. 



Figure 7. Predictions of performance relative 
to experimental findings 
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PROJECT GOALS 

The proposed research will investigate the stability and cryogenic properties of solid propellants 
that are critical to NASA’s goal of realizing practical propellant designs for future spacecraft. 

We will determine the stability and thermal properties of a solid hydrogen-liquid helium 
stabilizer in a laboratory environment in order to design a practical propellant. In particular, we 
will explore methods of embedding atomic species and metallic nano-particulates in hydrogen 
matrices suspended in liquid helium. We will also measure the characteristic lifetimes and 
diffusion of atomic species in these candidate cryofuels. 

The most promising large-scale advance in rocket propulsion is the use of atomic propellants; 
most notably atomic hydrogen stabilized in cryogenic environments, and metallized-gelled liquid 
hydrogen (MGH) or densified gelled hydrogen (DGH). The new propellants offer very 
significant improvements over classic liquid oxygen/hydrogen fuels because of two factors: (1) 
the high energy-release, and (ii) the density increase per unit energy release. These two changes 
can lead to significant reduced mission costs and increased payload to orbit weight ratios. An 
achievable 5-10% improvement in specific impulse for the atomic propellants or MGH fuels can 
result in a doubling or tripling of system payloads. 

The high-energy atomic propellants must be stored in a stabilizing medium such as solid 
hydrogen to inhibit or delay their recombination into molecules. The goal of the proposed 
research is to determine the stability and thermal properties of the solid hydrogen- liquid helium 
stabilizer. Magnetic resonance techniques will be used to measure the thermal lifetimes and the 
diffusive motions of atomic species stored in solid hydrogen grains. The properties of metallic 
nano-particulates embedded in hydrogen matrices will also be studied and analyzed. Dynamic 
polarization techniques will be developed to enhance signal/noise ratios in order to be able to 
detect low concentrations of the introduced species. The required lifetimes for atomic hydrogen 
and other species can only be realized at low temperatures to avoid recombination of atoms 
before use as a fuel. 

Goals for 2005 

I. Manufacture and characterization of particulates. The goal is to test different techniques for 
producing the solid hydrogen matrices; slow adiabatic cooling of gas mixtures versus rapid 
injection of hydrogen gas into a pressurized cold helium gas cell. We propose to determine 
the particle size, diffusion constants and lifetimes of species introduced into matrices 
prepared by each method to determine optimum conditions for producing a stable cryogenic 
environment. 
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II. Introduction of metallic species into the hydrogen particles . We will test methods of 
introducing metals into the hydrogen matrices by using two different techniques: (i) radio- 
frequency (RF) discharge dissociation of molecular hydrogen in the cryogenic injection path, 
and (ii) evaporation of metallic atoms from hot filaments in the hydrogen/helium gas stream 
in the condensation stream. Magnetic resonance (MR) techniques are the most effective for 
studying the atomic species in-situ because the MR frequency of each species is unique and 
can be studied independently. By this means we can determine the stability and obtain 
information about the surroundings of different atomic species. 


ACCOMPLISHMENTS 

Samples of finely divided hydrogen and methane suspended in liquid helium have been made. 
Two different methods were tested: condensation from gas mixtures onto cold surface, and 
injection of Hi/He gas mixtures into cold helium gas 


The thermal characteristics, notably the relaxation times and thermal diffusion were determined 
form NMR measurements of the relaxation times T 2 and Ti that are directly related to the critical 
path of thermal stability and the internal diffusion of the host matrix. 


Methane samples were explored and found to be significantly more favorable than hydrogen 
slush., The two cryogenic matrices were compared by measuring Ti and T 2 on CH4/He slush and 
comparing with H2/He slush. 


IMPORTANT CONCLUSIONS 

Based on our experimental findings we conclude that solid H 2 is not suitable as a host matrix due 
to the strong surface interactions of ortho-hydrogen molecules with helium at the surface of the 
hydrogen granules. Since solid H 2 has a strong surface scattering due to its large electric 
quadrupole moment and its rotational motion that persists to low temperatures, we explored 
solid methane which has no quadrupole moment as a new matrix. 

The results determined from the analysis for the microscopic correlation times x leads to the 
following conclusions: 

For hydrogen slush the NMR relaxation rate is consistent with scattering at grain boundaries 
and/or surface (due to large quadrupole moment Q of hydrogen). 

For methane slush the NMR relaxation rate is consistent with the limiting factor in 
thermalization being internal diffusion as opposed to surface scattering. 

In terms of a host matrix for atomic propellants methane is therefore a better host than hydrogen 
because the grains of methane are better isolated from the helium bath. 


NEW RESULTS ON CH 4 / He SLUSH 

From spin-spin relaxation times for freshly prepared methane-helium slush at 4.2K we found two 
different exponents which show the existence of two time scales in the system. As methane has a 
very complex ground state at low temperature we carried out our measurement on CH 4 /Hc slush 
as the sample was aging with time. The new results show more than one time scales as the 
sample was aged. 
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Figurel. Spin-spin relaxation times for methane-helium slush at 4.2K for 5days old sample. Two 
different exponents show existence of two time scales in the system. 
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Figure2. Spin-spin relaxation times for methane-helium slush at 4.2K for 21 days old sample. 
Three exponents show three different time scales in the system. 
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PROJECT GOALS 

The research is devoted to developing efficient techniques to measure the ratio of the ortho and 
para hydrogen concentrations in gaseous samples. Special cells have been designed to measure 
the thermal conductivity of the gas mixtures at temperatures in the vicinity of 100 K. The 
thermal conductivity is very different for the two molecular species and this enables one to 
determine the ratio of ortho to para hydrogen very reliably. The effort has been focused on the 
need for a fast, high accuracy measurement using a simple compact cell. Tests of the design have 
been successful, and the last trimester of research was dedicated to optimizing an electronic 
readout package and demonstrating the sensitivity, stability, and speed of response of the 
ratiometer. We have also developed the basic design for a cryogen free ratiometer. We have also 
explored the development of high efficiency ortho-para converters. 

Methods of measuring the ortho-para hydrogen ratio is very important for the manufacture, 
transport and storage of liquid hydrogen, and therefore critical for the space missions of NASA 
where liquid hydrogen fuels are employed. Because of the large heat of ortho-para conversion 
(670 J/g), incomplete conversion during liquefaction results in severe losses due to boil-off 
during storage and transport. Most manufacturers use a two-stage catalytic converter, at 77 K, 
and at 20.3 K in the liquefying process. At 300 K, the equilibrium concentration is 75% ortho, 
while at the boiling point (20.3 K), the equilibrium ortho concentration is 0.21%. The 
conversion at 77 K can produce 50% para-TB, and the remaining conversion at 20.3 K consumes 
a significant fraction (~ 50%) of liquid hydrogen. From these considerations, one of the most 
important considerations in the production and storage of liquid hydrogen is the need for 
efficient cost-effective means of conversion from ortho to para hydrogen. 1 We are investigating a 
new conversion process using a highly efficient ortho-para hydrogen converter in the form of 
Cr 203 powder as the active element. The oxide is contained in a continuous heat exchanger that 
allows the heat of conversion to be removed by the hydrogen gas flow. 

ACCOMPLISHMENTS 

1 . We have modified the electronic readout circuit for the ortho-para ratiometer to optimize 
the speed of response of the thin-film circuit used to measure the thermal conductivity of 
the gas in the ortho-para hydrogen cell. The circuit diagram is given in Figure 1. The 
results of tests for the circuit to control current to maintain the resistance (temperature 
indicator) at a constant value is given in Table 1. The cell is shown in Figure 2. 

2. We have changed one of the critical parameters of the ortho-para hydrogen ratiometer, 
namely the gap between the thin film sensor and the cold plate from 0.75 mm to 1.575 
mm. The results for the new configuration are shown in Figure 3(c). The data is 
consistent with the change in geometry. The thermal conductivity of nitrogen gas has also 
been measured to provide a calibration test. The results show a similar pressure 
independent behavior for a wide range of pressure as seen in studies of hydrogen for 
similar current settings. The thermal conductivity data has shown a wide range of 
pressure independence up to 1 .4 bar. This region of independence is consistent with 
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theory and indicates the design is successful in eliminating the effects of convection. The 
new results for the test are shown in Figure 3(a) and an expanded view is shown in Figure 
3(b). We will also carry out measurements for different gases and different mixtures of 
ortho and para hydrogen. 

3. Research has been initiated to test the efficiency of CT 2 O 3 as a catalytic agent for rapid 
ortho-para conversion. An ultra-fine chromic oxide powder (~1 micrometer diameter) is 
being considered for temperatures ranging from 4 to 25 K. (Figure 4) 


Table 1: Bridge Controller Performance (Maintain Resistance as a Constant) 


Operating Point: 
Current Gain: 

6.9 

4.0 

Cell Current: 
Cell Resistance: 
Cell Power: 

70.01 mA 
10.05 Q 
49.26 mW 

Cell Current: 
Cell Resistance: 
Cell Power: 

69.89 mA 
10.09 Q 
49.29 mW 

Cell Current: 
Cell Resistance: 
Cell Power: 

70.13 mA 
10.01 Q 
49.23 mW 


Operating Point: 
Current Gain: 

6.9 

8.0 

Cell Current: 
Cell Resistance: 
Cell Power: 

69.71 mA 
10.05 Q 
48.84 mW 

Cell Current: 
Cell Resistance: 
Cell Power: 

69.39 mA 
10.09 Q 
48.58 mW 

Cell Current: 
Cell Resistance: 
Cell Power: 

70.03 mA 
10.01 Q 
49.09 mW 


172 


UNIVERSITY OF 

FLORIDA 



Current Supply 

Note PrograrmrgCcmronGm 


Voltage Supply 


J4 



n 


ORG 


, *15V 

YEL 


1 

GRN 


3 -15V 


tit 
HEADER 3 



Layout & Assembly Schematic 


1 Hie 

Thermal Conductivity Cell - Resistance Control Circuit 


'M 

Document Number 

Rev 

A 

SULLIVAN 

i»p 1 1. 


Figure 1: Schematic Circuit Diagram of Bridge Control for a Constant Resistance 
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Figure 2. Ratiometer Design and Construction 
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Measurement by the Constructed Ratiometer 


Input Power vs Pressure (Add Hydrogen & Pump out) 



♦ Pump Out Hydrogen 
■ Add Hydrogen 


Pressure (bar) 


Figure 3(a) 

Expanded View 



Figure 3(b) 
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Hydrogen &Nitrogen Gas Measurement 

(Gap between film and cold plate changed to 1.575 mm) 



Figure 3(c) 
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Figure 4. C^CE as a catalytic agent for rapid ortho-para conversion 
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PROJECT GOALS 

The major impetus in the recent surge in the study of solid hydrogen and in particular the study 
of chemical species isolated in solid hydrogen has been the possible use of the material as rocket 
fuel. The incentive to this was the suggestion that metal atoms could be used as performance 
improving additives to the current state of the art liquid oxygen/ liquid hydrogen propellant 
system. The major point of interest is the specific impulse I sp of the propellant system, which is 
a measure of the momentum transferred to the vehicle per mass of the working fluid expelled 
from the rocket engine, and is proportional to the mean velocity of the exhaust. The specific 
impulse is proportional to the specific enthalpy of the chemical reaction occurring in the engine, 
and so can be changed by introducing dopants. Table 1 shows the theoretical specific impulses 
and their percentage change for 5% mole concentration of atoms in solid hydrogen propellants 
compared with the current liquid hydrogen system, reprinted from Fajardo et al (1). 

Table 1: Theoretical specific impulses of atom doped solid hydrogen propellants 


Additive 

Isp 

% improvement 

None 

389 

- 

H 

407 

+5 

Li 

401 

+3 

Be 

451 

+16 

B 

470 

+21 

C 

469 

+21 

N 

414 

+6 

0 

411 

+6 

F 

398 

+2 

Na 

373 

-4 

Mg 

398 

+2 

Al 

425 

+9 

Si 

432 

+11 

Ti 

404 

+4 


Clearly the systems of interest in this study are beryllium, boron, carbon, silicon and possibly 
aluminum. The major point of note is that this is a theoretical study and assumes that there is no 
interaction between the atomic species and the solid hydrogen. Several questions are raised by 
this study. First, can the atomic species of interest be actually produced and entrapped in solid 
hydrogen at 2 or 4 K? Second, to what degree does the atom interact with the solid hydrogen, 
and are molecular species produced? 
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These questions may be answered by reviewing the open academic literature considering either 
studies of species in solid hydrogen or matrix isolation studies in noble gas matrices where the 
matrix is doped with a small amount of hydrogen. 


In our earlier reports we reviewed the open academic literature considering the entrapment of 
either boron or carbon in hydrogen matrices or in noble gas matrices where the matrix is doped 
with a small amount of hydrogen. In review, although it appears that both boron and carbon 
vaporizing species interact with hydrogen, boron vaporizes as an atomic or dimeric species 
whereas carbon vaporizes in a plethora of carbon chain species the nature of which is being 
studied extensively by many research groups using matrix isolation. See, for example, studies by 
Larson (2) and Meirs (3), including detailed studies using hydrogen as the isolating matrix. 


ACCOMPLISHMENTS 


Boron in Solid Hydrogen 

Boron and carbon are considered of interest due to the high percentage improvements in the I sp . 
Our first candidate for study was be boron. It is well known that boron vaporizes as boron atoms 
and the dimeric species B 2 (4), and it is believed that boron atoms may be easily trapped in solid 
hydrogen. A theoretical study (5) indicated that the reaction path 


B + H 2 — » BH 2 

was unlikely to occur. To date, two major studies have vaporized boron into solid hydrogen. 
Fajardo et. al. (6) used ultraviolet \ visible spectroscopy to study the process and noted that 
although the resulting spectra were similar to those observed in rare gas matrices, there were 
significant frequency shifts. In the field of matrix isolation spectroscopy, major frequency shifts 
between matrices can sometimes be an indicator of the isolated species interacting with the 
surrounding matrix. The second study vaporizing boron into solid hydrogen was done by 
Weltner et. al . (7) using electron spin resonance (esr) as the identifying technique. The sole 
species identified was boron tetrahyride (BH 4 ); no boron atoms were observed. The esr spectrum 
of boron atoms is complex but very distinctive and has been the subject of much study (8,9). 

There is clearly a conflict in the open literature, with one study indicating the possibility of boron 
atoms, and the other study concluding that all of the boron reacts with hydrogen to form BH 4 . 
Supplementary clues are obtained by the studies of Andrews et al (10,1 1). In a detailed study 
Andrews and Tague (10) vaporized boron in argon matrices doped with amounts of hydrogen 
and found that a plethora of boron hydride species were formed. Increasing the hydrogen 
content and trapping in neon resulted in the production of the species BH 4 and BH 4 ‘ (1 1); all of 
these species were characterized using infrared spectroscopy. The literature evidence suggests 
that boron readily reacts with hydrogen to possibly form a variety of species but predominately 
hydrogen boron tetrahydride. In addition, boron readily reacts with trace impurities in the matrix 
gas (see for example discussion in ref (10)). This is known and has been studied extensively, but 
may increase the difficulty in interpreting results if ignored. 

Carbon in Solid Hydrogen 

Studies of the vaporization of graphite and the subsequent entrapment of the vapor in hydrogen 
have followed the pattern of the noble gas matrices. Van Zee et. al. (12) vaporized graphite and 
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using esr spectroscopy identified the species C 4 in solid hydrogen, whereas by using deuterium 
obtained spectra for C4, C6, Cs and C10 . Using infrared spectroscopy Shida et. al. (13) identified 
the species C3 C5 C9 and other clusters. It is interesting to note that no hydrocarbon species were 
observed. Tam et. al. (14) vaporizing graphite into hydrogen matrices also observed the carbon 
species C3 ,C 6 ,and C9, as well as the hydrocarbon species CH 4 and CH 3 . It was noted that the 
CH4: C3 ratio increased dramatically as the laser ablation energy increased. This is consistent 
with the mass spectrometric study by Kato et. al. (15) of the vaporization of graphite in the 
presence of hydrogen and deuterium gasses. The conclusion of this study was that hydrocarbon 
species originate mainly from the reaction of carbon species and the thermally decomposed 
hydrogen atoms produced by the intense vaporization pulses. 

The major conclusion from the literature on carbon species entrapped in solid hydrogen is that 
there is a general consensus on which carbon species are isolated and that the infrared spectra are 
of a very high quality. This was the main reason that we chose to make our initial studies on the 
boron / hydrogen system. 


Apparatus 

In our previous reports we described the construction and initial testing of our cryostat. In 
summary, cool down of the polished aluminum support from room temperature to circa 4 K was 
achieved routinely in less than ten minutes. Figure (1) shows a cross sectional view of the 
lower section of the cryostat showing the positioning of the boron sample relative to the polished 
aluminum sample support and the access window for the laser. Figure (2) shows a typical 
experiment in progress. The cryostat (a) is mounted on a custom made optical bench, which 
enables consistent positioning of the cryostat relative to the external sampling port of the Bruker 
infrared spectrometer (b). The movable pumping station (c), consisting of a diffusion pump and 
a backing rotary pump, is connected to the cryostat via a flexible vacuum line. The computer 
console (d) runs the spectrometer and analyzing spectra. The panel (e) is the hydrogen handling 
and purification system. Liquid helium is transferred from the liquid helium dewar (f) via the 
flexible transfer line (g) to the cryostat; exhaust gas is collected by the flexible recovery line (h), 
which is connected to the permanent laboratory recovery system (i). The Nd:YAG laser (j) with 
steering optics have been temporarily moved aside. The movable pumping station with flexible 
vacuum line (c) and the flexible recovery line (h) are features that enable the rotation of the top 
section of the cryostat and the inner cold finger under cryogenic and high vacuum conditions 
The spectrometer used is a Bruker IFS-1 13 with p.c. control, the laser used is a Nd:YAG pulsed 
q switched and frequency doubled, and is used under conditions known to vaporize the vast 
majority of materials. 

Deposition of Solid Hydrogen at 4 K 

To produce a sample for study, a solid hydrogen film was formed on the polished aluminum 
support by cooling the cryostat down to 4 K and then slowly bleeding hydrogen into the 
hydrogen inlet line using a needle valve. Figure (3) shows a typical infrared spectrum obtained 
after about twenty minutes deposition, without utilizing the purification system. The infrared 
spectrum exhibits the fundamental vibrations for carbon dioxide and monomeric water (16), in 
addition to weak features labeled (*) due to hydrocarbons. All three species are known to be 
trace impurities in many gasses including hydrogen, and carbon dioxide and water have been 
characterized in solid hydrogen (17,18). The bands assigned to carbon dioxide and water are 
particularly striking as those molecules are very strong infrared absorbers and are often used as 
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internal calibrants. Utilizing the cold trap in the purification apparatus removes all of these 
features apart from a very weak absorbance at about 2345 cm' 1 , which is the strong fundamental 
for carbon dioxide. 


Vaporization and Trapping of Boron in Solid Hydrogen at 4 K 

As shown in figure (1) the sample of solid boron is mounted in the vacuum chamber in such a 
way that the plume of laser ablated material is mixed with the large excess of hydrogen before 
being deposited on the cold polished support. Typically a thin bed of solid hydrogen is laid 
down on the polished surface prior to commencement of laser ablation. Infrared spectra are 
acquired every ten minutes of laser ablation in order to observe band growth and to assist 
spectral assignment. Occasionally laser ablation lasts up to one hour but is usually twenty to 
thirty minutes in duration. A typical experiment lasts several hours with no noticeable 
deterioration in spectral quality or relative band intensity changes. 

Figure (4) shows a typical early infrared spectrum produced after the deposition of boron vapor 
and hydrogen without purifying the hydrogen, as usual residual impurities water, carbon dioxide 
and hydrocarbons are observed. However, two new bands are observed in the region 2250 to 
2000 cm ' 1 that are within the B-H stretching region. The feature at about 2017 cm ' 1 is identified 
as a hydrolysis species such as HOBO. Later experiments exhibit spectra typical of figure (5) in 
which the feature at 2017 cm ' 1 is absent but a feature at about 474 cm ' 1 consistent with the 
species HBO is observed. In all of the experiments involving boron, a sharp band at 2138 cm ' 1 is 
observed. When this band is strong a feature at 1305 cm ' 1 appears to be growing with it. 
Purifying the hydrogen prior to deposition eliminates the infrared absorptions due to water and 
virtually eliminates the feature at 2345 cm ' 1 assignable to the strong fundamental of carbon 
dioxide. However, the features at 2138 and 1305 cm ' 1 are still produced. 

We tentatively assign the infrared absorption at 2138 cm ' 1 to the chemical species BH 4 in 
addition to the feature at 1305 cm 1 . Andrews and Wang (11) assign four infrared absorptions to 
BH4 in neon. However, it is believed they have a distorted tetrahedral system consisting of two 
short and two long B-H bonds with the bond angle between the two long B-H bonds a mere 46.7 
degrees. This structure is consistent with the structure proposed by Weltner et. al. (7) who 
characterized BH4 in hydrogen at 2K. We are proposing that BH4 in hydrogen at 4K is more like 
the classic tetrahedral with all four B-H bonds essentially equivalent. It is known that solid 
hydrogen at 4K is a lot softer than the same solid at 2 K, which would enable the species to relax 
to the tetrahedral we are postulating. The difference in softness between 2 and 4K is one of the 
reasons that the majority of researchers operate at 2 K. Clearly further work is required to clarify 
the situation. 

Vaporization and Trapping of Carbon in Solid Hydrogen at 4 K 

A small chunk of high purity graphite was mounted in our sample chamber in a similar manner 
to that of boron, was ablated using the Nd:YAG laser, and the resulting vapor isolated in a 
hydrogen matrix. Figure ( 6 ) shows a typical infrared spectrum after about 30 minutes of 
vaporization / deposition. By comparison with the earlier literature studies, we can readily 
identify the carbon species C 3 ,Ce ,C 5 ,and C9 , as well as a very weak absorbtion assignable to 
the species CFL . It is important to note that we did not observe the species CH4 which has been 
extensively studied in hydrogen matrices (19) and was observed by Tam et al (14). Of interest 
are the two weak features marked * and may be the species C(H 2 ). This species has been the 
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subject of extensive interest and questions have been recently raised ( 20 ) concerning how atomic 
carbon interacts with molecular hydrogen, whether it forms an adduct species CXFL ), or if it 
goes all the way and inserts into the hydrogen molecule to form the molecular species CFL . The 
positions and relative intensities of the two infrared bands will help to tackle this question as a 
dipole induced by the carbon atom in the hydrogen molecule may induce infrared activity. 


Continuous Flow Liquid Helium Cryostat 

We have constructed a ‘bare bones’ continuous flow liquid helium cryostat designed for infrared 
studies of solid hydrogen and isolated molecular species at either 2 or 4 K. A detailed 
description of the cell construction and testing has been published (21). The use of a heater we 
installed will enable studies up-to about 6 K, which opens up the possibility of diffusion and 
stability studies. To date, all of our studies have been at 4K which is a testament to the skill and 
expertise of the Physics Department Machine Shop and Cryogenic Services in the construction 
and design of the cell. 

We have been able to routinely trap high temperature chemical species in solid hydrogen at 4 K 
using this apparatus and produce stable matrices for study at leisure. The majority of researchers 
usually operate at 2 K when trapping species in hydrogen due to the fickle nature of solid 
hydrogen. A temperature of 2 K is obtained by pumping on the liquid helium coolant with a 
roughing vacuum pump resulting in evaporative cooling, but at the expense of higher liquid 
helium consumption. 

In this project we have demonstrated the vaporization of both boron and carbon and trapped the 
resulting vapors in solid hydrogen as initially proposed by us. In addition, preliminary analysis 
of the infrared spectra, indicate that molecular species such as boron hydrides and carbon chains 
or clusters are formed rather than solely atomic species. 

This work was briefly mentioned in an article in Cold Facts, the magazine of the Cryogenic 
Society of America Inc., authored by Van Sciver (22). In the article our aims, goals, and 
techniques were gratefully mentioned. However, the article claimed that our work was the first to 
embed carbon and boron in solid hydrogen. This is not the case. Many groups have preceded 
our work (Weltner, Fajardo, Andrews, Momose, and Shida to list a few), and we have benefited 
greatly from them. 

Further Studies 

The theoretical specific impulses calculated by Carrick (23) and quoted by Fajardo et al (1) 
assumed the ability to produce solely atomic species at a concentration of 5% in solid hydrogen. 
This would result in a significant improvement in the specific impulse of the rocket. We have 
demonstrated that the vaporization of either boron or carbon results in the generation of a 
plethora of molecular species, which is entirely consistent with previous studies of the systems in 
noble gas matrices and mass spectrometric studies. 

In order to accurately access the validity of high energy densified fuels based on impurity doped 
solid hydrogen as the new generation of rocket fuels, calculations need to be made on solid 
hydrogen doped with molecular species such as C 3 or BH 4 , rather than the imagined atomic 
species. 
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The next step in this work would be to improve the quality of the infrared spectra and to confirm 
the spectral assignments. The use of isotopes such as "C, C, B, B, H, and "H (D) is one of 
the most powerful tools available to the matrix isolationist for spectral identification. Further 
analysis of the data will be needed to order to attempt to quantify the relative concentrations of 
the chemical species, in particular the carbon system. The possible presence of C(H 2 ) in the 
carbon system offers the exciting possibility of being able to monitor the interaction of an atomic 
species with hydrogen, whether it is with the bulk hydrogen or a single hydrogen molecule. 

The use of infrared spectroscopy to monitor the generation of species produced by the reaction of 
boron atoms and dimmers with trace impurities such as water and carbon dioxide (not just 
hydrogen) is of immense value as impurities will always be present and they should not be 
ignored. 

We are at present unable to pursue any of this exciting research, as our funding was not renewed. 
We are aware that this research offers few if any immediate applications or financial savings to 
NASA, but it is much more difficult to dismiss the long term applications of such exciting 
research results. 
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Figure (1):- Cross Sectional View of the Cryostat 
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Figure (2):-Typical Experimental Setup (See text for explanation of symbols) 



Figure 2. Typical experiment in progress 

The cryostat (a) is mounted on a custom made optical bench, which enables consistent 
positioning of the cryostat relative to the external sampling port of the Bruker infrared 
spectrometer (b). The movable pumping station (c), consisting of a diffusion pump and a 
backing rotary pump, is connected to the cryostat via a flexible vacuum line, (d) is the computer 
console for running the spectrometer and analyzing spectra, (e) is the hydrogen handling and 
purification system. Liquid helium is transferred from the liquid helium dewar (f) via the 
flexible transfer line (g) to the cryostat, exhaust gas is collected by the flexible recovery line (h), 
which is connected to the permanent laboratory recovery system (i). The Nd:YAG laser (j) with 
steering optics have been temporarily moved aside. The movable pumping station with flexible 
vacuum line (c) and the flexible recovery line (h) are features, which enable the rotation of the 
top section of the cryostat and the inner cold finger under cryogenic and high vacuum conditions 
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Figure ( 3 ):- Infrared Spectrum of Solid Hydrogen Film at 4 K 
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Figure (4):- Infrared Spectrum Produced by Laser Ablating Boron and Trapping in Solid Hydrogen at 4 K 
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Figure (5):- Infrared Spectrum Produced by Laser Ablating Boron and Trapping in Solid Hydrogen at 4K 
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Figure (6):- Infrared Spectrum Produced by Laser Ablating Carbon and Trapping in Solid Hydrogen at 4 K 
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